差商法.cpp

数值分析的差商求解法,是用C编写的,希望大家互相参考!!1

#include<stdlib.h>
#include<stdio.h>


/*void main()
{
	int i,j,k,count,N;
	float T[10][10];
	cout<<"请输入T[i][j]的组数，不得超过20组:";//要求用户输入数据组数
    cin>>count;
    for(k=0;k<count;k++)
	{
		printf("请输入第%d组T[%d][%d]=",k+1,k,k);
        scanf("%f",&T[k][k]);
        printf("请输入第%d组T[%d][%d]=",k+1,k,k+1);
        cin>>T[k][k+1];
	}
	cout<<"请输入阶数N:";
	cin>>N;
	for(i=1;i<N;i++)
	{for(j=2;j<N;j++)
	T[i][j]=(T[i][j-1]-T[i-1][j-1])/(T[i][0]-T[i-1][0]);
	
	 cout<<T[i][j];
	}
	
}*/
int n=2;
void main()
{
	int i,n=2;
	//scanf("%d",&n);
	float x[n],y[n];
	for(i=0;i<=n;i++)
	{scanf("%f",&x[i]);
	scanf("%f",&y[i]);
	}
chashang(x,y);
}
 void chashang(float x[],float y[])
{
	int i,j;
	int n=2;
	float T[n+1][n+1],l,g,z;
	for(i=0;i<=n;i++)
		T[i][0]=y[i];
	for(j=1;j<=n;j++)
		for(i=j;i<=n;i++)
			T[i][j]=(T[i][j-1]-T[i-1][j-1])/(x[i]-x[i-j]);
		for(i=0;i<=n;i++)
		{
			for(j=0;i<=n;j++)
				if(i>=j) printf("%f\t",T[i][j]);
				printf("\n");
		}
		/*pirntf("z=?");
		scanf("%f",&z);
		l=0;
		for(i=0;i<=n;i++)
		{
			g=a[i][i];
			for(j=0;j<=i-l;j++)
				g=g*(z-x[j]);*/


}