二次插值法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]; ...
二次插值法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(fr
f2)
{ 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;
}
本文档为【二次插值法C语言程序】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。