如何調用和使用函數在matlab中的遞歸性

Question

我有一個牛頓的運動方式計算y(x+1)=y(x)-f(x)/f'(x)在這個函數我需要y（x）和爲此我使用函數遞歸y（1）& y（2）它的工作，因爲y（1）有公式y(1)=R*T/p，用於保存y（x）我使用zeros()來計算y（x + 1）時使用，但對於x> 2我得到相同的答案，我錯過了什麼？我可以用什麼來代替零（）根據@patrik指導保存和訪問蠑螈（X-1）

function y= newt(x) 
%define beta,gamma,delta,....there 
y(1)=R*T/p; 
answ=zeros(1,20); 
z=0; 
if x==1 
    f=(R*T*y(1)^3)+(beta*y(1)^2)+(gamma*y(1))+delta-(p*y(1)^4);   
    f1=(3*y(1)^2*R*T)+(2*y(1)*beta)+gamma+(4*p*y(1)^3);   
    answ(1) = y(1); 
    fprintf('n=1 v=%f\n',y(1));  
else 
    y=newt(x-1); 
    f=(R*T*y^3)+(beta*y^2)+(gamma*y)+delta-(p*y^4); 
    f1=(3*y^2*R*T)+(2*y*delta)+gamma+(4*p*y^3); 
    z=y-f/f1 
    answ(1,2:x)=z; 
end 
    answ(1) = y(1); 
    answ(1,2:x)=z; 

Answer 1

，我發現我的答案

f=@(v)p*v^4-R*T*v^3-beta*v^2-gamma*v-delta; 
f1=@(v)4*p*v^3-3*R*T*v^2-2*beta*v-gamma;