十大解法c语言Word格式文档下载.docx
- 文档编号:19448309
- 上传时间:2023-01-06
- 格式:DOCX
- 页数:12
- 大小:17.97KB
十大解法c语言Word格式文档下载.docx
《十大解法c语言Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《十大解法c语言Word格式文档下载.docx(12页珍藏版)》请在冰豆网上搜索。
i,n;
x[20],y[20],xx,yy;
printf("
Input
n:
"
);
scanf("
%d"
&
n);
if(n>
=20)
{printf("
Error!
The
value
of
n
must
in
(0,20)."
getch();
1;
if(n<
=0)
x[%d]:
i);
%f"
x[i]);
\n"
y[%d]:
y[i]);
xx:
xx);
yy=lagrange(x,y,xx,n);
x=%f,y=%f\n"
xx,yy);
2.牛顿插值多项式,用于离散数据的拟合
void
difference(float
*y,int
n)
*f;
k,i;
f=(float
for(k=1;
k<
=n;
k++)
f[0]=y[k];
k;
f[i+1]=(f[i]-y[i])/(x[k]-x[i]);
y[k]=f[k];
return;
difference(x,(float
*)y,n);
yy=y[20];
for(i=n-1;
i>
=0;
i--)
yy=yy*(xx-x[i])+y[i];
NewtonInter(%f)=%f"
3.高斯列主元消去法,求解其次线性方程组
code#include<
math.h>
#define
N
20
n,i,j,k;
mi,tmp,mx;
a[N][N],b[N],x[N];
\nInput
N)
input
should
in(0,N)!
=0)
Now
a(i,j),i,j=0...%d:
n-1);
n;
a[i][j]);
b(i),i,j=0...%d:
b[i]);
n-2;
for(j=i+1,mi=i,mx=fabs(a[i][j]);
n-1;
if(fabs(a[j][i])>
mx)
mi=j;
mx=fabs(a[j][i]);
if(i<
mi)
tmp=b[i];
b[i]=b[mi];
b[mi]=tmp;
for(j=i;
tmp=a[i][j];
a[i][j]=a[mi][j];
a[mi][j]=tmp;
for(j=i+1;
tmp=-a[j][i]/a[i][i];
b[j]+=b[i]*tmp;
for(k=i;
a[j][k]+=a[i][k]*tmp;
x[n-1]=b[n-1]/a[n-1][n-1];
for(i=n-2;
i--)
x[i]=b[i];
x[i]-=a[i][j]*x[j];
x[i]/=a[i][i];
Answer:
\n
x[%d]=%f\n"
i,x[i]);
0;
#include<
NUMBER
Esc
0x1b
Enter
0x0d
A[NUMBER][NUMBER+1]
ark;
flag,n;
exchange(int
r,int
k);
max(int
message();
{
x[NUMBER];
r,k,i,j;
char
celect;
clrscr();
\n\nUse
Gauss."
\n\n1.Jie
please
press
Enter."
\n\n2.Exit
Esc."
celect=getch();
if(celect==Esc)
exit(0);
\n\n
n="
\n\nInput
matrix
A
and
B:
for(i=1;
a%d1--a%d%d
b%d:
i,i,n,i);
for(j=1;
=n+1;
j++)
A[i][j]);
k++)
ark=max(k);
if(ark==0)
\n\nIt'
s
wrong!
else
if(flag!
=k)
exchange(flag,k);
for(i=k+1;
for(j=k+1;
A[i][j]=A[i][j]-A[k][j]*A[i][k]/A[k][k];
x[n]=A[n][n+1]/A[n][n];
for(
k=n-1;
k>
=1;
k--)
me=0;
me=me+A[k][j]*x[j];
x[k]=(A[k][n+1]-me)/A[k][k];
\n\nx%d=%f"
k)
i;
A[0][i]=A[r][i];
A[r][i]=A[k][i];
A[k][i]=A[0][i];
temp=0;
for(i=k;
if(fabs(A[i][k])>
temp)
temp=fabs(A[i][k]);
flag=i;
temp;
message()
Go
on
Exit
Esc!
switch(getch())
case
Enter:
main();
Esc:
default:
error!
4.龙贝格求积公式,求解定积分
f(x)
(sin(x)/x)
MAX
20
a
2
b
4
e
0.00001
LBG(float
p,float
q,int
sum=0,h=(q-p)/n;
for
(i=1;
sum+=f(p+i*h);
sum+=(f(p)+f(q))/2;
return(h*sum);
n=N,m=0;
T[MAX+1][2];
T[0][1]=LBG(a,b,n);
n*=2;
for(m=1;
m<
MAX;
m++)
m;
T[i][0]=T[i][1];
=m;
T[i][1]=T[i-1][1]+(T[i-1][1]-T[i-1][0])/(pow(2,2*m)-1);
if((T[m-1][1]<
T[m][1]+e)&
&
(T[m-1][1]>
T[m][1]-e))
Answer=%f\n"
T[m][1]);
;
code5.牛顿迭代公式,求方程的近似解
100
PS
1e-5
TA
Newton(float
(*f)(float),float(*f1)(float),float
x0
)
x1,d=0;
k=0;
do
x1=
x0-f(x0)/f1(x0);
if((k++>
N)||(fabs(f1(x1))<
PS))
\nFailed!
exit();
d=(fabs(x1)<
1?
x1-x0:
(x1-x0)/x1);
x0=x1;
x(%d)=%f\n"
k,x0);
while((fabs(d))>
PS&
fabs(f(x1))>
TA)
x1;
f(float
x)
x*x*x+x*x-3*x-3;
f1(float
3.0*x*x+2*x-3;
f(float);
f1(float);
x0,y0;
x0:
x0);
x(0)=%f\n"
x0);
y0=Newton(f,f1,x0);
\nThe
root
is
x=%f\n"
y0);
6.
牛顿-科特斯求积公式,求定积分
NC(a,h,n,r,f)
(*a)[];
h;
n,f;
*r;
nn,i;
ds;
1000||n<
2)
if
(f)
Faild!
Check
1<
n<
1000!
n);
return(-1);
if(n==2)
*r=0.5*((*a)[0]+(*a)[1])*(h);
return(0);
(n-4==0)
*r=0;
*r=*r+0.375*(h)*((*a)[n-4]+3*(*a)[n-3]+3*(*a)[n-2]+(*a)[n-1]);
if(n/2-(n-1)/2<
nn=n;
else
nn=n-3;
ds=(*a)[0]-(*a)[nn-1];
for(i=2;
=nn;
i=i+2)
ds=ds+4*(*a)[i-1]+2*(*a)[i];
*r=ds*(h)/3;
nn)
h,r;
n,ntf,f;
a[16];
the
x[i](16):
=15;
a[i]);
h=0.2;
f=0;
ntf=NC(a,h,n,&
r,f);
if(ntf==0)
\nR=%f\n"
r);
Wrong!
Return
code=%d\n"
ntf);
7.雅克比迭代,求解方程近似解
0.00001
i,j,k;
t;
a[N][N],b[N][N],c[N],g[N],x[N],h[N];
dim
0<
N!
a[i,j],i,j=0…%d:
c[i],i=0…%d:
c[i]);
b[i][j]=-a[i][j]/a[i][i];
g[i]=c[i]/a[i][i];
h[j]=g[j];
for(k=0;
if(j==k)
continue;
h[j]+=b[j][k]*x[k];
t=0;
if(t<
fabs(h[j]-x[j]))
t=fabs(h[j]-x[j]);
x[j]=h[j];
e)
x_i=\n"
i++)
after
%d
repeat
return\n"
MAX);
8.秦九昭算法
qin(float
a[],int
n,float
r=0;
for(i=n;
r=r*x+a[i];
r;
a[50],x,r=0;
n,i;
frequency:
while(n<
1);
value:
x);
r=qin(a,n,x);
9.幂法
3
x[n]={0,0,1};
a[n][n]={{2,3,2},{10,3,4},{3,6,1}};
y[n];
xm,oxm;
oxm=0;
N;
for(j=0
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 解法 语言