a2=(Ft2-F2)./(Q(2)*m); t2=trapz(ua2,1./a2)/3.6;
ua3=uamax(2):uamax(3);
n3=ua3*Ig0(k)*Ig5(3)/(0.377*0.367);
Ttq3=-19.313+295.27*(n3/1000)-165.44*(n3/1000).^2+40.874*(n3/1000).^3-3.8445*(n3/1000).^4;
Ft3=Ttq3*Ig0(k)*Ig5(3)*0.85/0.367;
F3=0.013*3880*9.8+2.77*ua3.^2/21.15; a3=(Ft3-F3)./(Q(3)*m); t3=trapz(ua3,1./a3)/3.6;
ua4=uamax(3):70;
n4=ua4*Ig0(k)*Ig5(4)/(0.377*0.367);
Ttq4=-19.313+295.27*(n4/1000)-165.44*(n4/1000).^2+40.874*(n4/1000).^3-3.8445*(n4/1000).^4;
Ft4=Ttq4*Ig0(k)*Ig5(4)*0.85/0.367;
F4=0.013*3880*9.8+2.77*ua4.^2/21.15; a4=(Ft4-F4)./(Q(4)*m); t4=trapz(ua4,1./a4)/3.6; t(k)=t2+t3+t4; end
f=polyfit(Qz,t,2);
QQ=Qz(1):0.001:Qz(5); b=polyval(f,QQ) plot(QQ,b,Qz,t,'*')
gtext('5.17') gtext('5.43') gtext(