TR(1)=10 ;
TR(2)=17.76 ;
TR(3)=25.52 ;
TR(4)=33.28 ;
TR(5)=41.04 ;
TR(6)=48.8 ;
TR(7)=56.56 ;
TR(8)=64.32 ;
TR(9)=72.08 ;
TR(10)=79.84 ;
TR(11)=87.6 ;
TD(1)= 5.5908700e+03 ;
TD(2)= 5.2324300e+03 ;
TD(3)= 4.8857600e+03 ;
TD(4)= 4.5508700e+03 ;
TD(5)= 4.2277500e+03 ;
TD(6)= 3.9164100e+03 ;
TD(7)= 3.6168300e+03 ;
TD(8)= 3.3290300e+03 ;
TD(9)= 3.0530100e+03 ;
TD(10)=2.7887500e+03 ;
TD(11)=2.5362700e+03 ;
Tk(1)= 6.1434300e+11 ;
Tk(2)= 5.3482100e+11 ;
Tk(3)= 4.6326700e+11 ;
Tk(4)= 3.9913100e+11 ;
Tk(5)= 3.4188300e+11 ;
Tk(6)= 2.9101100e+11 ;
Tk(7)= 2.4602700e+11 ;
Tk(8)= 2.0645700e+11 ;
Tk(9)= 1.7185100e+11 ;
Tk(10)=1.4177600e+11 ;
Tk(11)=1.1582000e+11 ;
x_Tgang=(TR-10)/77.6;
faiT_r_out = -2.504*x_Tgang.^6 +6.2275*x_Tgang.^5 -5.6202*x_Tgang.^4 +2.1963*x_Tgang.^3 +0.7004*x_Tgang.^2;
faiT_r_in = 0.5357*x_Tgang.^6 -2.2395*x_Tgang.^5 +3.0871*x_Tgang.^4 -1.7684*x_Tgang.^3 +1.385*x_Tgang.^2;
faiT_r_out_dot =(-2.504*6*x_Tgang.^5 +6.2275*5*x_Tgang.^4-5.6202*4*x_Tgang.^3 +2.1963*3*x_Tgang.^2 +0.7004*2*x_Tgang.^1)/77.6;
faiT_r_in_dot =(0.5357*6*x_Tgang.^5 -2.2395*5*x_Tgang.^4+3.0871*4*x_Tgang.^3 -1.7684*3*x_Tgang.^2 +1.385*2*x_Tgang.^1)/77.6;
faiT_r_out_dot2=(-2.504*6*5*x_Tgang.^4+6.2275*5*4*x_Tgang.^3-5.6202*4*3*x_Tgang.^2+2.1963*3*2*x_Tgang.^1+0.7004*2)/(77.6^2);
faiT_r_in_dot2 =(0.5357*6*5*x_Tgang.^4-2.2395*5*4*x_Tgang.^3+3.0871*4*3*x_Tgang.^2-1.7684*3*2*x_Tgang.^1+1.385*2)/(77.6^2);
Mn_in= trapz(TR,TD.*(faiT_r_in.^2)) +240000+56780;
Mn_out=trapz(TR,TD.*(faiT_r_out.^2))+240000+56780;
Kn_out=trapz(TR,Tk.*(faiT_r_out_dot2.^2));
Kn_in=trapz(TR,Tk.*(faiT_r_in_dot2.^2));
(Kn_in/Mn_in)^0.5/2/pi
(Kn_out/Mn_out)^0.5/2/pi
I would be appreciated it if you could reply. Thank you
Best regards
Yingxin