Question

我在这个函数中用牛顿方式计算y(x+1)=y(x)-f(x)/f'(x)我需要y（x），为此我使用y（1）＆amp;的函数Recursivity。 y（2）它正在工作，因为y（1）具有公式y(1)=R*T/p，对于保存y（x）我在计算y（x + 1）时使用zeros()但是对于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;

如何在matlab中调用和使用函数递归

1 个答案: