clear all, clc
Tarray=[0.01,1,6,100];
for j=1:4
    T=Tarray(j);    
    X(j,:)=[0:1:1000];
    Y(j,:)=exp(-X(j,:)/T);
end

figure()
h = axes('GridLineStyle','--',...
       'FontName','times',...
       'FontAngle','italic',...
       'FontSize',16);
plot(X(1,:),Y(1,:),'b','LineWidth',2), hold on
xlabel('recovery time [months]')
ylabel('distance out of equilibrium [m]')
axis([0 100 0 (Y(j,:)/Y(1,:))])
plot(X(2,:),Y(2,:),'r','LineWidth',2)
plot(X(3,:),Y(3,:),'g','LineWidth',2)
plot(X(4,:),Y(4,:),'k','LineWidth',2)
legend('T=0.01','T=1','T=6','T=100')

saveas(gcf,['parameter T.png'])

