Как вызвать и использовать функцию рекурсивности в matlab

Question

У меня есть упражнение с Ньютоном способ расчета y(x+1)=y(x)-f(x)/f'(x) в этой функции мне нужен y (x) и для этого я использую функцию рекурсивности для y (1) и y(2), она работает, потому что y (1) имеет формулу y(1)=R*T/p, для сохранения у (х) я использую zeros()использовать при вычислении y(x+1), но для x>2 я получаю тот же ответ, что мне не хватает?? что я могу использовать вместо нулей () для сохранения и доступа к newt(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

Согласно руководству @patrik, я нахожу свой ответ

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;