lisp抛物线插值_抛物线插值

#include&ltstdio.h&gt

#include&ltmath.h&gt

#define m 10

void main()

{

int i,n;

float x[m],y[m],xgiven,yresult,yresult10,yresult20,yresult30,yresult11,yresult21,yresult31;

printf("\n&lt抛物线插值法&gt\n");

printf("请输入已知点个数(少于%d个):\n",m);

scanf("%d",&n);

printf("\n请输入已知坐标点:\n");

for(i=0;i&ltn;i++)

{

printf("\n点(%d)=",i);

scanf("%f%f",&x[i],&y[i]);

}

printf("\n请输入求解坐标点:");

scanf("%f",&xgiven);

if(xgiven&ltx[0]||xgiven&gtx[n-1])

printf("越界\n");

else

{

for(i=1;i&ltn-1;i++)

{

if(xgiven&lt=x[i])

{

if(i=1)

break;

else if(fabs(xgiven-x[i-1])&lt=fabs(x[i]-xgiven))

{

i=i-1;

break;

}

}

}

if(x[n-2]&ltxgiven&&xgiven&lt=x[n-1])

i=n-2;

yresult10=y[i-1]*(xgiven-x[i])*(xgiven-x[i+1]);

yresult11=(x[i-1]-x[i])*(x[i-1]-x[i+1]);

yresult20=y[i]*(xgiven-x[i-1])*(xgiven-x[i+1]);

yresult21=(x[i]-x[i-1])*(x[i]-x[i+1]);

yresult30=y[i+1]*(xgiven-x[i-1])*(xgiven-x[i]);

yresult31=(x[i+1]-x[i-1])*(x[i+1]-x[i]);

yresult=yresult10/yresult11+yresult20/yresult21+yresult30/yresult31;

printf("xgiven=%f,yresult=%f\n",xgiven,yresult);

}

}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值