數(shù)值計(jì)算方法實(shí)驗(yàn)

《數(shù)值計(jì)算方法實(shí)驗(yàn)》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、1.已知√（100）=10,√(121)=11,√(144)=12,用拉格朗日插值法求√（115）。C語言代碼如下/*Note:YourchoiceisCIDE*/#include"stdio.h"#defineN3voidmain(){doublex[N]={100,121,144};doubley[N]={10,11,12};doubles,t,y1=0,x1=115.0;intj,k;for(k=0;k

2、t;}printf("每一項(xiàng)的插值是：%lf",t);y1=y1+t*y[k];}printf("y的值是：%f",y1);}然后得到即√（115）≈10.7227562.用迭代法求x^3-x-1=0的一正根。迭代方程：x[k+1]=(xk+1)^(1/3)C語言代碼如下：/*Note:YourchoiceisCIDE*/#include"stdio.h"#include"math.h"voidmain(){doublex,x1;x=1.5;do{x1=x;x=pow((x1+1),1.0/3);}while(fabs(x-x1)>1e-9);print

3、f("x=%4.5f",x);}即x=1.324723.用牛頓迭代法求((2x-4)x+3)x-6=0，x0=1.5C語言代碼如下/*Note:YourchoiceisCIDE*/#include"stdio.h"#include"math.h"voidmain(){floatx1,x0,f,f1;x1=1.5;do{x0=x1;f=((2*x0-4)*x0+3)*x0-6;f1=(6*x0-8)*x0+3;x1=x0-f/f1;}while(fabs(x1-x0)>=1e-5);printf("x=%8.7f",x1);}輸出結(jié)果即x的取值為24.用二分法

4、求f(x)=x^3-3x^2+3x-2的根?？梢运愠鰂’=3x^2-6x+3=3(x-1)^2即f(x)單調(diào)，即有一個(gè)根C語言代碼如下/*Note:YourchoiceisCIDE*/#include"stdio.h"#include"math.h"#definef(x)(x*x*x-3*x*x+3*x-2)voidmain(){doublea=-10,b=10,c,eps=1e-5;inti=0;if(f(a)*f(b)>0){printf("請(qǐng)重新輸入a和b的值");}while((b-a)>eps){i++;c=(a+b)/2;if(f(c)==0)b

5、reak;elseif(f(a)*f(c)<0)b=c;elsea=c;}printf("二分次數(shù)為%d",i);printf("root=%f",c);}輸出結(jié)果它的根準(zhǔn)確為2相對(duì)跟為1.9999985.用高斯消去法解決方程組10x1-7x2+0x3=7的值。5x1-x2+5x3=6-3x1+2x2+6x3=4C語言編程如下：/*Note:YourchoiceisCIDE*/#include"stdio.h"#defineN3voidmain(){floata[N][N];floatb[N];floats,akk,aik;intk,i,j;prin

6、tf("請(qǐng)輸入原始增廣矩陣(按行輸入):");for(i=0;i

7、or(i=k+1;i=0;i--){s=0.0;for(j=2;j>i;j--)s=s+a[i][j]*b[j];b[i]=b[i]-s;}for(i=0;i

8、{printf("x[%d]=%ft",i,b[i