求助：牛顿迭代法计算程序

cqrj

处女作， 可惜不成功  ，运行时说“  θ=5;  ”参数有问题   望老师们开导  谢谢了 ！！



clear;
x=8;
y=0;
k=0;
p=31.512678;
A=105;
cotw=0.836649;
θ=5;
for i=1:6
θ=θ-((x*x+k*x*y+p*A*cotw)*cos(θ*3.14159/180)-(x*y+k*y*y+p*A*cotw*k)*sin(θ*3.14159/180)-p*p*θ*sin(θ*3.14159/180)-p*p*k*θ*cos(x*3.14159/180)-(A*x+p*x*cotw+k*A*y+k*p*y*cotw))/
-(x*x+k*x*y+p*A*cotw)*sin(θ*3.14159/180)-(x*y+k*y*y+p*A*cotw*k)*sin(θ*3.14159/180)-p*p*(sin(θ*3.14159/180)+θ*cos(θ*3.14159/180))-p*p*k*(cos(θ*3.14159/180)-θ*sin(θ*3.14159/180))   
end