myfun.m

自己遍的牛顿迭代法求解二元二次方程的程序 迭代步数可以自己设定

function y = myfun(x)
t=35;% x值 15，20，25，30，35
y(1)=x(1);
y(2)=4/pi*(1/3*(2.5-(1/(0.011998*t+x(1)))^2)*sqrt(1-(1/(0.011998*t+x(1)))^2)+(1/(0.011998*t+x(1)))*asin(1/(0.011998*t+x(1))))

y=[y(1),y(2)];