example5_6.java

清华大学2002年出版的《科学与工程数值计算算法Java》配套源码

import java.applet.*;
import java.awt.*;
public class Example5_6 extends Applet
{
 public double eepqs(double x,double h,int n,double[] y,double t)
 {int i,j,k,l;
  double[] b=new double[n];
  double f;
  b[0]=y[0];
  for(i=1;i<n;i++)
  {b[i]=y[i];
   for(j=0;j<i;j++)b[i]=(i-j)*h/(b[i]-b[j]);
  }
  f=b[n-1];
  for(i=n-2;i>=0;i--)f=b[i]+(t-x-i*h)/f;
  return f;
 }

 public void paint(Graphics g)
 {double[] y=new double[11];
  double f,h;int i;
  f=-1;
  for(i=0;i<11;i++){y[i]=1.0/(1+25*f*f);f+=0.2;}
  f=eepqs(-1,0.2,11,y,-0.25);h=1.0/(1+25*0.25*0.25);
  g.drawString("t=-0.85  f(t)="+f+"  实际值="+h,10,20);
  f=eepqs(-1,0.2,11,y,0.75);h=1.0/(1+25*0.75*0.75);
  g.drawString("t=0.25   f(t)="+f+"  实际值="+h,10,40);
  }
}