如何找到每個初始點

Question

function p = newton_hw(p0,tol,Nmax) 

%NEWTON'S METHOD: Enter f(x), f'(x), x0, tol, Nmax 
f = @(x) x*cos(x)-((sin(x))^2); 
fp= @(x) -x*sin(x)+ cos(x)-2*sin(x)*cos(x); 
p = p0 - (f(p0)/fp(p0)); 
y1=f(p); 
fprintf('y1=%f',y1) 
i = 1; 
while (abs(p - p0) >= tol) 

    p0 = p; 
    p = p0 - f(p0)/fp(p0); 
    i = i + 1; 
    if (i >= Nmax) 
     fprintf('Fail after %d iterations\n',Nmax); 
     break 
    end 

    y=f(p); 
    fprintf('a=%f,y=%f,\n',p,y); 
end 

end 

的牛頓迭代這是我的問題：

如何循環訪問爲每個p0 = 0,.1,.2,...,49,5。

Answer 1

使用迭代用戶的步驟可以做到這樣：

for i = 0:0.1:5 

而且任何索引可以做同樣的方式：x = [0:2:50]。

如果你的功能才能正常工作（我想是這樣），我們可以這樣說：

k = 1; 
for i = 0:0.1:5 
res(k) = newton_hw(i,0.001,1000); 
k = k+1; 
end 

但是，我們也可以在一行做到這一點 - 風格：

res = arrayfun(@(x) newton_hw(x, 0.001, 1000), I)