main.cpp

本人自编的共轭梯度法。为c源码

#include "iostream.h"
const max_count=100;
double x[max_count-1];
double matrix[max_count-1][max_count-1];
double work_outx;
int count;
void main()
{
int i,j,t;
double sum,temp;
while (true)
{
cout<<"输入NewTon插值计算阶数:"<<endl;
cin>>count;
     for(i=0;i<=count;i++)
	 {
          cout<<"X"<<i<<"值:"<<endl;
		  cin>>x[i];
          cout<<"X"<<i<<"值所对应的f(x)值:"<<endl;
		  cin>>matrix[i][0];
	 }
	 cout<<"输入插值后计算的X值:"<<endl;
	 cin>>work_outx;
	 		 for (t=1;t<=count;t++)  //列
			 {
				 for(j=0;j<=count;j++)  //行
				 {
				     if (x[j+t]-x[j]!=0) 
					 {
					 matrix[j][t]=(matrix[j+1][t-1]-matrix[j][t-1])/(x[j+t]-x[j]);
					 }
				 }
			 }
			  for (i=0;i<=count;i++)
				 cout<<"插值系数"<<i+1<<":"<<matrix[0][i]<<endl;
			  sum=matrix[0][0];
			for (j=1;j<=count;j++)
			{
			  temp=1;
			  for (i=0;i<=j-1;i++)
			  {
				  temp=temp*(work_outx-x[i]);
			  }
			  temp=temp*matrix[0][j];
			  sum=sum+temp;
			}
			cout<<"所求值: "<<sum<<endl;
}
}