这是我的一段程序,但是运行之后发现圆相当粗,为什么？


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');