求助！电机的Te-s曲线绘制

clc;clear;
Pn=10000;m=3;p=2;f=50;Un=380;Te=zeros(1,100);
R1=1.33;X1=2.43;R2=1.12;X2=4.4;Rm=7;Xm=90;Z1=R1+i*X1;Zm=Rm+i*Xm;
c=1+X1/Xm;
Ns=3000/p;
smax=c*R2/sqrt(R1^2+(X1+c*X2)^2);
Tmax=m*Un^2/((2*pi*Ns/60)*(2*c*(R1+sqrt(R1^2+(X1+c*X2)^2))));
Tst=m*Un^2*R2/((2*Pi*Ns/60)*((R1+c*R2)^2+(X1+c*X2)^2));
Ist=Un/Xm+Un/sqrt((R1+c*R2)^2+(X1+c*X2)^2);
n1=[0,200,400,600,810,1110,1210,1260,1310,1360,1400,1450,1470];
s=(1500-n1)/1500;
L=length(s);
for i=1:L
Te(i)=(m*Un^2*(R2/s(i)))/(((R1+c*R2/s(i))^2+(X1+c*X2)^2)*(2*pi*Ns/60));
end
figure(1);
plot(n1,Te,'b-',n1,Te,'r*');
xlabel('转速n（r/min）');
ylabel('转矩Te(N.m)');
grid on;
fprintf('\n');
fprintf('临界转差率：smax=%5.3f\n',smax);
fprintf('最大转矩：Tmax=%5.2fN.m\n',Tmax);
fprintf('启动电流：Ist=%5.2fA\n',Ist);
fprintf('启动转矩：Tst=%5.3fN.m\n',Tst);