czas = 500000;



Tp = 0.3*10^5;
To = 0.02*10^5;
delta_x = 0.49;
k = delta_x/dfmw;

L=[k];
M=[Tp,1];

sim('Lab3_pelen_domek1', czas);

s = tf('s');
test = (k/(Tp*s+1))*exp(-s*To);


hold on

plot(time,Test2);
%plot(time,Twew1)