matlab实验9答案详解.docVIP

实验九 数值微积分与方程数值求解 一、实验目的 掌握求数值导数和数值积分的方法 掌握代数方程数值求解方法 掌握常微分方程数值求解方法 二、实验内容 1.程序代码 x=1; i=1; f=inline(det([x x^2 x^3;1 2*x 3*x^2;0 2 6*x])); while x=3.01 g(i)=f(x); i=i+1; x=x+0.01; end g; t=1:0.01:3.01; dx=diff(g)/0.01; f1=dx(1) f2=dx(101) f3=dx(length(g)-1) 运行结果 f1 = 6.0602 f2 = 24.1202 f3 = 54.1802 程序代码 f=inline(sqrt(cos(t.^2)+4*sin(2*t).^2+1)); I1=quad(f,0,2*pi) g=inline(log(1+x)./(1+x.^2)); I2=quad(g,0,2*pi) 运行结果 I1 = 10.4285 I2 = 0.9997 程序代码 A=[6 5 -2 5;9 -1 4 -1;3 4 2 -2;3 -9 0 2]; b=[-4 13 1 11]; x=A\b y=inv(A)*b [L,U]=lu(A)