Pmec=90.3;
t1=5.66;
omg0=151.63;
Ji=Pmec*t1/omg0^2
kv=0.003262;
mf=(Ji/t1-kv)*omg0
tau=Ji/kv
t=0:20;
omg1=(omg0+mf/kv)*exp(-t/tau)-mf/kv;
plot(t,omg1)
grid