> restart;
> rivn := c*Udv = Tm*Te*diff(diff(w(t),t),t)+Tm*diff(w(t),t)+w(t);

> Tm:=0.2;Te:=0.02;c:=1;Udv:=1;Ky:=1000;i:=2;S:=0.01;







> Tm := .7;

> ans:=dsolve(rivn);

> plot([t,1-1.06927*exp(-1.521124*t)-0.06924*exp(-23.478875*t),t=0..5],0..3,0..1.2);

> plot([Re(Complex((c*Ky*S/i)/(Tm*Te*(I*w)^3+Tm*(I*w)^2+(I*w)^1))),Im(Complex((c*Ky*S/i)/(Tm*Te*(I*w)^3+Tm*(I*w)^2+(I*w)^1))),w=0..10],-10..10,-10..10);

> X:=w->-Tm*w^2+1;
> Y:=w->-Tm*Te*w^3+w;
> plot([X(w),Y(w),w=0..10],X=-40..3,Y=-10..5);



> plot([w,sqrt(X(w)^2+Y(w)^2),w=0..10]);

> arg:=w->-Pi-arctan(w);
> plot([w,arg(w),w=-20..20],w=-20..20,Phi=0..-2*Pi);
>
>
>


> Xroz:=w->X(w)*S*Ky/i;Yroz:=w->Y(w)*S*Ky/i;


> plot([Xroz(w),Yroz(w),w=0..100]);

> plot([w,sqrt(Xroz(w)^2+Yroz(w)^2),w=0..10]);

> plot([w,arctan(Yroz(w),Xroz(w)),w=0..40]);

> plot([log10(w),20*log10(sqrt(Xroz(w)^2+Yroz(w)^2)),w=0..10]);

> plot([log10(w),arctan(Yroz(w),Xroz(w)),w=0..40]);

> k:=1.7;
> Xzam:=w->1.2*k-1.2*k*Re((Xroz(w)+I*Yroz(w))/(1+Xroz(w)+I*Yroz(w)));
> Yzam:=w->k*Im((Xroz(w)+I*Yroz(w))/(1+Xroz(w)+I*Yroz(w)));
> plot([Xzam(w),Yzam(w),w=0..100],x=-1.2..1.2,y=-1.2..1.2);
>




> plot([w,sqrt(Xzam(w)^2+Yzam(w)^2),w=0..10]);

> plot([w,arctan((Yzam(w)),(Xzam(w))),w=0..40]);

> plot([log10(w),20*log10(sqrt(Xzam(w)^2+Yzam(w)^2)),w=0..10]);

> plot([log10(w),arctan(Yzam(w),Xzam(w)),w=0..40]);

>
> plot([1-1.06927*exp(-1.521124*t)-0.06924*exp(-23.478875*t),1-exp(-t/1.2),1],t=0..6,color=[black,black,black], style=[line,line,point],thickness=[5,1,1]);