一
x=10+round(rand(1,10)*90); x1=diff(x) x2=diff(x,2) x3=diff(x,3)
二
x=1;i=1; a='det([x x^2 x^3;1 2*x 3*x*x;0 2 6*x])'; %原行列式,题目求x=1,2,3时矩阵的数值导数 f=inline(a); while x<=3.01 g(i)=f(x); i=i+1; x=x+0.01; %因为差分的使用使维数减少 end g; dx=diff(g)/0.01; f1=dx(1) f2=dx(101) f3=dx(length(g)-1)
三
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)
四
function f=test9_4_1(x,y) f=1./sqrt(x.^2+y.^2); I1=integral2(@test9_4_1,0,1,0,1)
function f=test9_4_2(x,y) f=y.*sin(x)+3.*cos(y)-1; I2=integral2(@test9_4_2,0,2*pi,0,pi)
五
x=0.3:0.2:1.5; Y=[0.3895 0.6598 0.9147 1.1611 1.3971 1.6211 1.8325]; trapz(x,Y)%梯形公式
x=0.3:0.2:1.5; Y=[0.3895 0.6598 0.9147 1.1611 1.3971 1.6211 1.8325]; s=sum(Y(1:end-1).*(diff(x)))%矩形公式