求助！哪位高手帮我看哈这个程序！！！

somnus

解一个6元1次方程组，
U=f[(x-d)cosr+ysinr](x^2+y^2+e^2-a^2)-ex[(x-d)^2+y^2+f^2-c^2]
V=f[(x-d)sinr-ycosr](x^2+y^2+e^2-a^2)+ex[(x-d)^2+y^2+f^2-c^2]
W=2efsinr[x(x-d)+y^2-dycotr]
U+V=W
取正弦函数sin(pix)上的6个点（0,0）；（1/6,1/2）； [1/4，sin(pi/4)];   [1/3,sin(pi/3)];   (1/2,1)  ;  (1 ,0) 代入上式得到6元1次方程组，求解未知数  a ,c , d, e, f, r,谢谢了！
function p=mybar(u1,u2,u3,u4,u5,u6)
u1=f;
u2=d;
u3=r;
u4=e;
u5=a;
u6=c;
p(1)=u1*((0-u2)*cos(u3)+0*sin(u3))*(0^2+0^2+u4^2-u5^2)-u4*0*((0-u2)^2+0^2+u1^2-u6^2)+u1*((0-u2)*sin(u3)-0*cos(u3))*(0^2+0^2+u4^2-u5^2)+
u4*0*((0-u2)^2+0^2+u1^2-u6^2)-2*u4*u1*sin(u3)*(0*(0-u2)+0^2-u2*0*cot(u3);

p(2)=u1*((1/6-u2)*cos(u3)+1/2*sin(u3))*(1/6^2+1/2^2+u4^2-u5^2)-u4*1/6*((1/6-u2)^2+1/2^2+u1^2-u6^2)+u1*((1/6-u2)*sin(u3)-1/2*cos(u3))*(1/6^2+1/2^2+u4^2-u5^2)+
u4*1/2*((1/2-u2)^2+1/6^2+u1^2-u6^2)-2*u4*u1*sin(u3)*(1/6*(1/6-u2)+1/2^2-u2*1/2*cot(u3);

p(3)=u1*((1/4-u2)*cos(u3)+sin(pi/4)*sin(u3))*(1/4^2+sin(pi/4)^2+u4^2-u5^2)-u4*1/4*((1/4-u2)^2+sin(pi/4)^2+u1^2-u6^2)+u1*((1/4-u2)*sin(u3)-sin(pi/4)*cos(u3))*(1/4^2+sin(pi/4)^2+u4^2-u5^2)+
u4*1/4*((1/4-u2)^2+sin(pi/4)^2+u1^2-u6^2)-2*u4*u1*sin(u3)*(1/4*(1/4-u2)+sin(pi/4)^2-u2*sin(pi/4)*cot(u3);

p(4)=u1*((1/3-u2)*cos(u3)+sin(pi/3)*sin(u3))*(1/3^2+sin(pi/3)^2+u4^2-u5^2)-u4*1/3*((1/3-u2)^2+sin(pi/3)^2+u1^2-u6^2)+u1*((1/3-u2)*sin(u3)-sin(pi/3)*cos(u3))*(1/3^2+sin(pi/3)^2+u4^2-u5^2)+
u4*1/3*((1/3-u2)^2+sin(pi/3)^2+u1^2-u6^2)-2*u4*f*sin(u3)*(1/3*(1/3-u2)+sin(pi/3)^2-u2*sin(pi/3)*cot(u3);

p(5)=u1*((1/2-u2)*cos(u3)+1*sin(u3))*(1/2^2+1^2+u4^2-u5^2)-u4*1/2*((1/2-u2)^2+1^2+u1^2-u6^2)+u1*((1/2-u2)*sin(u3)-1*cos(u3))*(1/2^2+1^2+u4^2-u5^2)+
u4*1/2*((1/2-u2)^2+1^2+u1^2-u6^2)-2*u4*u1*sin(u3)*(1/2*(1/2-u2)+1^2-u2*1*cot(u3);

p(6)=u1*((1-u2)*cos(u3)+0*sin(u3))*(1^2+0^2+u4^2-u5^2)-u4*1*((1-u2)^2+0^2+u1^2-u6^2)+u1*((1-u2)*sin(u3)-0*cos(u3))*(1^2+0^2+u4^2-u5^2)+
u4*1*((1-u2)^2+0^2+u1^2-u6^2)-2*u4*u1*sin(u3)*(1*(1-u2)+0^2-u2*0*cot(u3);
p=[p(1),p(2),p(3),p(4),p(5),p(6)];