niufront43.java

数值分析之计算方法与实验4：插值法 by java

class Niufront43{
  static double [][] Chafen(double a[][]){
    int n=a.length;
    for(int j=2;j<=n+1;j++)
      for(int i=0;i<(n-j+1);i++)
         a[i][j]=a[i+1][j-1]-a[i][j-1];
    return a;     
  }
  static void N3(double t,double a[][]){
  	double y=a[0][1]+a[0][2]*t+a[0][3]*t*(t-1)/2+a[0][4]*t*(t-1)*(t-2)/6;
  	System.out.println(y);
  }	
  	
  public static void main(String[]args){
    double a[][]={{20.0,1.30103,0.0,0.0,0.0,0.0},{21.0,1.32222,0.0,0.0,0.0},{22.0,1.34242,0.0,0.0},{23.0,1.36173,0.0},{24.0,1.38021}};
    a=Chafen(a);
    for(int i=0;i<a.length;i++){
      for(int j=0;j<a.length+1-i;j++)
         System.out.print(a[i][j]+"   ");
       System.out.println();
    }
    N3(1.4,a);
    
  }
}