问一下,大侠帮帮忙啊~~
这是我的一段程序,但是运行之后发现圆相当粗,为什么?if xm=='Y'
p=0:360;
x=input('x方向全应力---- (MPa)');
y=input('y方向全应力---- (MPa)');
t=input('xy平面切应力---- (MPa)');
a=input('截面应力角---- (度)');
m=0.5*(x+y)+0.5*(x-y)*cos(2*a)+t*sin(2*a);
n=-0.5*(x-y)*sin(2*a)+t*cos(2*a);
a2=atan((2*t)/(x-y));
M=0.5*(x+y)+(((0.5*(x-y))^2+a^2)^0.5)*cos(p);
N=(((0.5*(x-y))^2+a^2)^0.5)*sin(p);
disp(n);
disp(a2);
plot(M,N,'.r');