科学计算与MATLAB语言-第六单元.pdf

第六单元

1求定积分的近似值

f=@(x)sqrt(cos(x.^2)+1);

%下面三种方法选择一种即可

%基于全局自适应计分方法

I1=integral(f,0,2*pi)

I1=6.0401

%辛普森方法

I2=quad(f,0,2*pi)

I2=6.0401

%Guss-Lobatto

I3=quadl(f,0,2*pi)

I3=6.0401

%梯形积分法

clear

f=@(x)sqrt(cos(x.^2)+1);

a=0;

b=2*pi;

n=1000;

h=(b-a)/n;

x=linspace(a,b,n+1);

y=f(x);

I_trapz=h*(sum(y)-0.5*(y(1)+y(end)))

I_trapz=6.0401

2在给定的初值x0=1,y0=1,z0=1下，求方程组的数值解

%方法一，定义独立的函数，下面这段代码需要一函数的形式存在单独的函数文件中

functionF=F0502(X)

x=X(1);

y=X(2);

z=X(3);

F(1)=sin(x)+y^2+log(z)-7;