R1=1.5958 ;
H1=3;

pi=3.14;

c1=1.6903 ;



R2=1.85901;
H2=10;
H2max=H2;



c2=0.638877 ;









%charakterystyka statyczna modelu nieliniowego
q1=[0:0.00001:1.40];
%h20=h2/2 => q=1,43
h2=((q1.*q1)/(c2*c2));




h20=0.5*H2max;
q10=sqrt(h20)*c2
h10=(c2/c1)^2*h20;


deltah2=0.1*h20;

deltaq=((c2/(2*sqrt(h20)))*deltah2)/2;




h2z=((2*sqrt(h20)/c2))*(q1-q10)+h20;



 %xlabel('q1');
 %ylabel('h2');
%plot(q1,h2,'b');
%hold on;
%plot(q1,h2z,'r');

s=tf('s');

A=pi*R1*R1;
B=c1/(2*sqrt(h10));
C=(pi*R2*R2*h20*h20)/(H2*H2);
D=(c2)/((2*sqrt(h20)));

K=B/((A*s+B)*(C*s+D))


T1=pi*R1*R1*2*sqrt(h10)/c1
T2=2*sqrt(h20)*R2*R2*h20*h20*pi/(c2*H2*H2)
k=2*sqrt(h20)/c2

G=tf(K);

chce zrobic wykres od czasu h2 i h2zlin
plot(tout,h2);
hold on;
plot(tout,h2z);