求值

			#include "StdAfx.h"
			#include <stdio.h>
			#define MAX 10
			void  single(){
				double x1,x2,y1,y2,f,x;
				printf("请再输入一组x,y(先输入x再输入y)");
				scanf("%lf%lf",&x1,&y1);
				printf("请再输入一组x,y");
				scanf("%lf%lf",&x2,&y2);
				printf("请再输入您需要求解的自变量x");
				scanf("%lf",&x);
				f=y2*(x-x1)/(x2-x1)+y1*(x-x2)/(x1-x2);
				printf("该函数值为%lf",f);


			}
			void dou(){
				double x1,x2,x3,y1,y2,y3,f,x;
				printf("请再输入一组x,y(先输入x再输入y)");
				scanf("%lf%lf",&x1,&y1);
				printf("请再输入一组x,y");
				scanf("%lf%lf",&x2,&y2);
				printf("请再输入一组x,y");
				scanf("%lf%lf",&x3,&y3);
				printf("请再输入您需要求解的自变量x");
				scanf("%lf",&x);
				f=y1*(x-x2)*(x-x3)/((x1-x2)*(x1-x3))+y2*(x-x1)*(x-x3)/((x2-x1)*(x2-x3))+y3*(x-x2)*(x-x1)/((x3-x2)*(x3-x1));
				printf("该函数值为%lf",f);


			}

			void nev(){
			int n,i,j;
			double a[MAX],b[MAX],c[MAX][MAX],x;
			printf("请输入您将输入的数据的组数(注意:不能超过10组)");
			scanf("%d",&n);
			while(n>10||n<0){
				printf("请重新输入您将输入的数据的组数(注意:不能超过10组)");
				scanf("%d",&n);
			}
			for(i=0;i<n;i++){
				printf("请输入x y");
				scanf("%lf%lf",&a[i],&b[i]);
			}
			printf("请输入您需要求得的x");
			scanf("%lf",&x);
			for(i=0;i<n;i++){
				for(j=0;j<n;j++){
					c[i][j]=0;
				}
				c[i][0]=b[i];
			}
			for(i=1;i<n;i++){
				for(j=1;j<n;j++){
					c[j][i]=c[j-1][i-1]+(c[j][i-1]-c[j-1][i-1])/(a[j]-a[j-1])*(x-a[j-1]);
				}
			}
			for(i=0;i<n;i++){
				for(j=0;j<=i;j++){
					printf("%lf  ",c[i][j]);
				}
				printf("\n");
			}
			}
			void main(){
				 single();//线性插值法
				 dou();//抛物插值法
				 nev();//Neville算法
			}