c数值分析c16-5.pdfVIP

732 Chapter 16. Integration of Ordinary Differential Equations dy[j]=yest[j]; } else { for (k=1;kiest;k++) fx[k+1]=x[iest-k]/xest; for (j=1;j=nv;j++) { Evaluate next diagonal in tableau. v=d[j][1]; d[j][1]=yy=c=yest[j]; h r P C S t e t a e o a p r m for (k=2;k=iest;k++) { : d m p / a y p / i r w b1=fx[k]*v; b s i l w l s g e b=b1-c; e i h p w f o t a . i n (