二次插值法C语言程序

二次插值法C语言程序 #include"stdio.h" #include"math.h" #include"conio.h" void main() { float*area(float a1,float p,float a[3]); float f(float x); float ar,fr; float a1=10,p=0.01,e=0.000001; float pa[3]; area(a1,p,pa); a1=pa[0]; float a2=pa[1]; float a3=pa[2]; float f1=f(a1); float f2=f(a2); float f3=f(a3); do { ar=((a3*a3-a2*a2)*f1+(a1*a1-a3*a3)*f2+(a2*a2-a1*a1)*f3); ar=ar/2/((a3-a2)*f1+(a1-a3)*f2+(a2-a1)*f3); fr=f(ar); if(ar>a2) { if(fr>f2) { a3=ar; f3=fr; } else if(frf2) { a1=ar; f1=fr; } else if(fr=f1) { if(fabs(f2-f1)a3) { temp=a1; a1=a3; a3=temp; } a[0]=a1;a[1]=a2;a[2]=a3; return a; } float f(float x) { float y=pow(fabs(x-1),1.5)+pow(fabs(x-1),2.7); return y; }