z(2)=(z(1)+(m)/(z(1)))/2;
i=1;
while
abs(z(i+1)-z(i))>=1.0000e-005
z(i+2)=(z(i+1)+(m)/(z(i+1)))/2;
i=i+1;
end
yf=z(i+1)
以上程序已經經過測試
別忘了：執行程序之前，先給m賦值，否則會報錯的！