數(shù)值計(jì)算方法丁麗娟課后習(xí)題答案

數(shù)值計(jì)算方法丁麗娟課后習(xí)題答案【篇一：北京理工大學(xué)數(shù)值計(jì)算方法大作業(yè)數(shù)值實(shí)驗(yàn)1】)書p14/4分別將區(qū)間[?10,10]分為100,200,400等份，利用mesh或surf命令畫出二元函數(shù)的三維圖形。z=|??|+??+??+??++??【matlab求解】[x,y]=meshgrid(-10:0.1:10);a=exp(-abs(x));b=cos(x+y);c=1-/(x.A2+y.A2+1);z=a+b+c;mesh(x,y,z);[x,y]=meshgrid(-10:0.05:10);a=exp(-abs(x));b=cos(x+y);c=1-/(x,A2+y,A2+1);z=a+b+c;mesh(x,y,z);[x,y]=meshgrid(-10:0.025:10);a=exp(-abs(x));b=cos(x+y);c=1./(x.A2+y.A2+1);z=a+b+c;mesh(x,y,z);(二)書p7/1?3?2數(shù)值計(jì)算的穩(wěn)定性⑴取=??c語言程序—不穩(wěn)定解+=ln1.2，按公式=??(n=1,2，…)#includestdio.h#includeconio.h#includemath.hvoidmain(){floatm=log(6.0)-log(5.0),n;inti;i=1;printf(y[0]=%-20f,m);while(i20){n=1/i-5*m;printf(y[%d]=%-20f,i,n);m=n;i++;if(i%3==0)printf(\n);}getch();}(ii)c語言程序一穩(wěn)定解???[??+??+????+??按公式=??(??)#includestdio.h#includeconio.h#includemath.hvoidmain(){floatm=(1/1O5?0+1/126?0)/2,n;k=n,n-1,n-2,…)(【篇二：北京理工大學(xué)數(shù)值計(jì)算方法大作業(yè)數(shù)值實(shí)驗(yàn)4】p260/1考紐螺線的形狀像鐘表的發(fā)條，也稱回旋曲線，它在直角坐標(biāo)系中的參數(shù)方程為=??????????????????????????????=????????????????曲線關(guān)于原點(diǎn)對稱，，取a=i/參數(shù)s的變化范圍[-5,5],，容許誤差限'分別是，,和。選取適當(dāng)?shù)墓?jié)點(diǎn)個(gè)數(shù)，利用數(shù)值積分方法計(jì)算曲線上點(diǎn)的坐標(biāo),并畫出曲線的圖形。【matlab求解】??x=zeros(i00,i);y=zeros(i00,i);func1=inline(cos(1/2*(t■八幼)；func2=inline(sin(1/2*(t.A2)));n=i;fori=-5:0.i:5x(n,i)=quad(funci,0,i,ie-6);y(n,i)=quad(func2,0,i,ie-i0);n=n+i;endplot(x,y,b);xlabel(x);ylabel(y);(代碼運(yùn)行圖)2■書p261/3二重積分的計(jì)算（1）i=????????????????????（????+）dy???（??+??）（2）i=??????????dy??+??+??【matlab求解】（1）首先，用matlab自帶的求解二重積分的命令對兩個(gè)二重積分分別求解【第一問解】【篇三：北京理工大學(xué)徐特立學(xué)院數(shù)值分析大作業(yè)上機(jī)實(shí)驗(yàn)】ss=txt＞機(jī)實(shí)驗(yàn)選做教材：數(shù)值計(jì)算方法（2011第一版）.丁麗娟,程杞元.高等教育出版社八以下代碼作者原創(chuàng)八超鏈接：1.2（，，，）2.2（，，，）3.1（，，，，）5.1（，，，）5.3（，，，）第一章：數(shù)值計(jì)算中的誤差2、題目簡介：利用pi/4=1-1/3+1/5-1/7。。。級數(shù)計(jì)算pi的近似值。輸入：誤差值輸出：求和項(xiàng)數(shù)，并輸出pi值工具：c語言運(yùn)行環(huán)境：vc-6.0計(jì)算公式及原理：利用pi/4=1-1/3+1/5-1/7oo。級數(shù)計(jì)算pi的近似值，由數(shù)學(xué)原理可知誤差會小于首次舍棄的項(xiàng),可以編寫循環(huán)實(shí)現(xiàn)。程序運(yùn)行結(jié)果截圖：程序代碼：（c語言）#includestdio.hvoidmain（）{printf（第一章第2題求pi,歡迎使用，請按提示操作。\n）;inti=1,n=0,k=1;doublee,pi,er;printf（請輸入誤差（例如1e-4）:）;scanf（%lf,e）;}printf(請稍候。。。\n);er=e;pi=O;while(er=e){}pi*=4;printf(%d項(xiàng)求和后可以達(dá)到％?10lf精度，這時(shí)pi=%.10lf\n,n,e,pi);getchar();getchar();pi+=k*1.0/i;k=-k;er=1.0/i;i+=2;n++;第二章：解線性方程組的直接方法2、題目簡介：用m