三维磁晶各向异性能（草稿）

(1) y=a*(sin(pi/180*x))^2+b*(sin(pi/180*x))^4     易001
    y=a*(cos(pi/180*x))^2+b*(cos(pi/180*x))^4  易010 100
得到a，b，即为r1，r2





(2) x(u,v)=(r1*sin(u)^2+r2*sin(u)^4)*sin(u)*cos(v)     易001
y(u,v)=(r1*sin(u)^2+r2*sin(u)^4)*sin(u)*sin(v)
z(u,v)=(r1*sin(u)^2+r2*sin(u)^4)*cos(u)



double r1 = 
double r2 = 



(2) x(u,v)=(r1*cos(u)^2+r2*cos(u)^4)*sin(u)*cos(v)    易010 100
y(u,v)=(r1*cos(u)^2+r2*cos(u)^4)*sin(u)*sin(v)
z(u,v)=(r1*cos(u)^2+r2*cos(u)^4)*cos(u)