机械优化设计二次插值法(上机)

二次插入法的c语言程序，机械优化设计（第四版）第二次上机内容哦，根据课本程序框图写的，分享给工科学生。

——CXK #define f(x) sin(x)//定义函数，函数可以任意改的。

程序：

#include

#include

#define f(x) sin(x)

void main()

{

double y1,y2,y3,x1,x2,x3,e,h,xr,yr,c1,c2,xp,yp;

printf("请输入x1,x2,x3,e,h\n");

scanf("%lf%lf%lf%lf%lf",&x1,&x2,&x3,&e,&h);

y1=f(x1);

y2=f(x2);

y3=f(x3);

do

{

c1=(y3-y1)/(x3-x1);

c2=((y2-y1)/(x2-x1)-c1)/(x2-x3);

xp=0.5*(x1+x3-c1/c2);

yp=f(xp);

if(fabs((y2-yp)/y2)

break;

if((xp-x2)*h>0)

{

if(y2>=yp)

{x1=x2;

y1=y2;

x2=xp;

y2=yp;}

else

{x3=xp;

y3=yp;}

}

else

{

if(y2>=yp)

{x3=x2;

y3=y2;

x2=xp;

y2=yp;}

else

{x1=xp;

y1=yp;}

}

}while (fabs((y2-yp)/y2)>=e);

if(y2

{xr=x2; yr=y2;

printf("\nxr=%f\nf(xr)=%f",xr,yr);} else

{ xr=xp; yr=yp;

printf("\nxr=%f\nf(xr)=%f",xr,yr);} system("pause") ;

return main();

}

下面是运行程序图：

喜欢的话记得给五星~