matlab数学实验复习题(有标准答案)
复习题
1、写出3
2、i ｎv(A)表示Ａ的逆矩阵；
3、在命令窗口健入
clc,４、在命令窗口健入clea
5、在命令窗口健入６、x=-1:０.2:17、det(A)表示计算Ａ的行列式的值；8、三种插值方法：拉格朗日多项式插值,分段线性插值，三次样条插值。
9、若A=123456789⎡⎤
⎢⎥⎢⎥⎢⎥⎣⎦,则fliplr （Ａ）＝321654987⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦
Ａ-３=210123456--⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦A ．＾２=149162536496481⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦
tril(A)=100450789⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦ tri ｕ（Ａ，－1)=123456089⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦ｄiag(A ）＝100050009⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦
Ａ(:,2),＝2
58A(3,:)=369
10、nor ｍcd ｆ（1，1，2)=0.5%正态分布mu=1,s ｉgm ａ＝２,x ＝1处的概率
ｅ45(@ｆ,［ａ,b ］,x0),中参数的涵义是@fun 是求解方程的函数M 文
件,[a,b ］是输入向量即自变量的范围a 为初值,ｘ0为函数的初值，t 为输出指定的[ａ,b],x 为函数值
１5、写出下列命令的功能：te ｘt （1,2,‘ｙ=s ｉｎ(x)’
hold on 16fun ｃｔion 开头;
17
,4)
３，４）
21、设x 是一向量，则）的功能是作出将Ｘ十等分的直方图 ２2、interp １(［1,2,3]，[3,４,5]，2.5)
Ans=4.５
23、建立一阶微分方程组⎩
⎨⎧+='-='y x t y y x t x 34)(3)(2
的函数M 文件。
(做不出来） 二、写出运行结果:
１、>＞ｅy ｅ（3，4)=1000
01000010
2、>>s ｉze([1,2,3])=1;3
3、设b=ro ｕnd （unifrnd(-５，5,１，4)),则=3 ５ 2 -5 >>[x,m]=ｍin(b)；x ＝－5；m=4
，［x,n ］=sort(b ）
－５ 2 3 ５
4 3 １ 2
ｍea ｎ(ｂ）=１.2５，m ｅdian(b)＝2.5,range(ｂ）=10
４、向量ｂ如上题,则
＞>ａn ｙ(b),all(ｂ<2),ａｌl(b<6）
Ans ＝１ 0 １
５、＞>[5 6;7 8]>[7 8；5 6]=00
11
6、若1234B ⎡⎤=⎢⎥⎣⎦
,则 7、>>ｄｉag(d ｉag （B ）)=10
04
8、＞>[4:-2:１].＊[-1,６]=-4 12
9、>>acos(0.5),a ｔan(1）
ans ＝
1.６5９8
ans=。
matlab复习题标准答案
matlab复习题标准答案# MATLAB复习题标准答案一、选择题1. MATLAB中的矩阵默认是按什么顺序存储的？- A. 行主序- B. 列主序- C. 随机序- D. 行主序2. 下列哪个函数用于计算矩阵的行列式？- A. `det`- B. `diag`- C. `trace`- A. `det`3. MATLAB中如何创建一个5x5的单位矩阵？- A. `eye(5,5)`- B. `zeros(5,5)`- C. `ones(5,5)`- A. `eye(5,5)`4. 下列哪个命令用于生成0到1之间的随机数？- A. `rand()`- B. `randi()`- C. `randn()`- A. `rand()`5. MATLAB中如何计算向量的范数？- A. `norm(vector)`- B. `abs(vector)`- C. `sum(vector)`- A. `norm(vector)`二、填空题1. MATLAB中的`%`符号用于______。
- 答案：注释2. 要创建一个从10到1的递减向量，可以使用______。
- 答案：`1:-1:10`3. MATLAB中，`zeros(3,4)`会生成一个______。
- 答案：3行4列的全零矩阵4. 函数`max(A)`会返回矩阵A中的______。
- 答案：最大值5. `plot(x,y)`函数在MATLAB中用于______。
- 答案：绘制x和y的数据点的图形三、简答题1. 请简述MATLAB中矩阵的基本操作有哪些？- 答案：矩阵的基本操作包括矩阵的创建、转置、求逆、矩阵乘法、矩阵加法、矩阵减法等。
2. MATLAB中如何实现循环结构？- 答案：MATLAB中实现循环结构主要有两种方式，一种是使用`for`循环，另一种是使用`while`循环。
3. 请解释MATLAB中的函数文件和脚本文件的区别。
- 答案：函数文件通常包含返回值和参数，可以被其他脚本或函数调用；而脚本文件主要用于执行一系列命令，不包含返回值。
（完整版）MATLAB）课后实验答案[1]
![（完整版）MATLAB）课后实验答案[1]](https://360docs-1393072327.cos.ap-beijing.myqcloud.com/306ff29b690203d8ce2f0066f5335a8102d266cb.webp)
（完整版）MATLAB）课后实验答案[1]实验⼀ MATLAB 运算基础1. 先求下列表达式的值，然后显⽰MATLAB ⼯作空间的使⽤情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =，其中2120.455i x +??=?- (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--L (4) 2242011122123t t z t t t t t ?≤=-≤，其中t =0:0.5:2.5 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建⽴⼀个字符串向量，删除其中的⼤写字母。
解：(1) 结果：(2). 建⽴⼀个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验⼆ MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S=?，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对⾓阵，试通过数值计算验证2 2E R RS A O S +??=。
解: M ⽂件如下；5. 下⾯是⼀个线性⽅程组：1231112340.951110.673450.52111456x x x ??=???(1) 求⽅程的解。
(2) 将⽅程右边向量元素b 3改为0.53再求解，并⽐较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解： M ⽂件如下：123d4e56g9实验三选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ?+-<≠-?=-+≤<≠≠??--?且且及其他⽤if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
解：M ⽂件如下：2. 输⼊⼀个百分制成绩，要求输出成绩等级A、B、C、D、E。
matlab20道试题及解答
试题1．“数学黑洞”：任意一个4位自然数，将组成该数的各位数字重新排列，形成一个最大数和一个最小数，之后两数相减，其差仍为一个自然数。
重复进行上述运算，最终会出现一个神秘的数，请编程输出这个神秘的数。
clear;a=input('请输入一个四位正整数：');str_a=num2str(a); %将a转化为一个字符串b_min=str2double(sort(str_a)); %形成最小数b_max=str2double(sort(str_a,'descend')); %形成最大数b=b_max-b_min; %求最大数与最小数之差while (b~=a)a=b;str_a=num2str(a); %将a转化为一个字符串b_min=str2double(sort(str_a)); %形成最小数b_max=str2double(sort(str_a,'descend')); %形成最大数b=b_max-b_min; %求最大数与最小数之差endb试题2．将数字1、2、3、4、5、6填入一个2行3列的表格中，要使得每一列右边的数字比左边的数字大，每一行下面的数字比上面的数字大。
请编写程序求出按此要求可有几种填写方法。
a(1)=1;a(6)=6;count=0; %用来计数b=perms('2345'); %产生2345的全排列[m,n]=size(b);for i=1:mtemp=b(i,:);a(2)=str2double(temp(1));a(3)=str2double(temp(2));a(4)=str2double(temp(3));a(5)=str2double(temp(4));if ((a(4)>a(2))&&(a(4)>a(3))&&(a(5)>a(3)))count=count+1;c=reshape(a,2,3); %将a向量转化为2*3矩阵输出disp(c);endenddisp(['共有',num2str(count),'种填写方法']); %输出填写方法的种数试题3．编写成绩排序程序。
数学实验(matlab版)过程考试试卷及答案完整版
试绘出三种产品产量与季度的三维垂直方向条形图(分组式). >> x=[8,8,9;11,7,8;12,6,9;10,6,10]; bar3(x,'group') 图形如下：
2/2
第一题：编程计算下面问题, x 值由键≥ 1 y = x 2 , −1 ≤ x < 1 2 x − 1, x < −1
>> x=input('输入 x：'); if x>=1 y=x^2+1; end if x<-1 y=x^2-1; end if x>=-1&x<1 y=x^2; end y 输入 x：5 y= 26 第二题：某人做一种材料的伸缩实验,t 为温度(℃),L 为长度(mm),实验数据见下表 t 20 25 30 35 40 L 81 82.3 84 86.8 89
f = x 4 − xy + y 2 ,求
>> syms x y
∂f ∂ 3 f , ∂x ∂y 3
1/1
f=x^4-x*y+y^2; dx=diff(f,x,1) dy3=diff(f,y,3) dx = 4*x^3-y dy3 = 0 第四题：某厂生产三种产品,某年四季度的产量如下 A 产品产量 笫一季度 笫二季度 笫三季度 笫四季度 8 11 12 10 B 产品产量 8 7 6 6 C 产品产量 9 8 9 10
用二阶拟合法,求 L 与 t 的表达式.要求:1.编程;2.写出 L 与 t 的关系式. >> t=[20,25,30,35,40]; L=[81,82.3,84,86.5,89]; k=polyfit(t,L,2) k= 0.0091 -0.1446 80.2114 L=0.0091 t^2 —0.1446t+ 80.2114 第三题：求微分与积分(编程)
MATLAB考试题目及解答
1.圆盘上有如图所示的二十个数，请找出哪4个相邻数之和最大，并指出它们的位置和数值。
（10分）2011841361015217319716811149125解答：%1.圆盘上有如图所示的二十个数，请找出哪4个相邻数之和最大，并指出它们的位置和数值。
（10分）A=[1 18 4 13 6 10 15 2 17 3 19 7 16 8 11 14 9 12 5 20];% 程序位置规定：从1开始顺时针方向计数；NumA=size(A); Num=NumA(1,2); sum(1)=A(1); for i=1:(Num-3)sum(i)=A(i)+A(i+1)+A(i+2)+A(i+3); endmaxresult=max(sum(:));%找出4个相邻数之和最大值 maxresult %4个相邻数之和最大值 Position=i %四个数起始位置 FourNumber=A(1,i:(i+3)) %四个数的值及顺序 运行结果：maxresult =50Position =17FourNumber =9 12 5 202.甲、乙、丙三人上街买糖果。
三人都买好后，甲对乙、丙说，我可以按你们现有的糖果数再送你们每人一份。
甲送给乙、丙后，乙也按甲、丙现有的糖果数，送给甲、丙每人各一份糖果。
丙也如此送了甲、乙各一份。
互相赠送后，每人恰好各有64颗糖果。
问甲、乙、丙原来各买了多少糖果？（10分）解答：%由代数关系构造矩阵 a=[1 -1 -1;0 2 0;0 0 2]; b=a([2 1 3],:); b=b(:,[2 1 3]); c=a([2 3 1],:); c=c(:,[2 3 1]); d=64*ones(3,1); result=a\(b\(c\d))运行结果：result=104 56 323.求n S a aa aaa aaa a =++++ 的值。
a 的值为1~9之间的一个整数，n S 中每一项aaa a 共有n 位。
MATLAb与数学实验 第五章习题解答
1 3 2 4 1 3 5 2 4 3 5
j=
1 1 2 2 3 3
3 4 4 5 5
s=
1 3 1 3 2 1 3 2 1 2 1 A1 =
10200 01020 30102 03010 00301 (2) n=nnz(A)
n=
13 ans =
1 2 -1 3 -2 4 1 -3 5 2
-4 3 -5
>> nx=nzmax(A)
nx =
35
>> [i,j,s]=find(A)
i=
1 2 1 3 2 4 1 3 5 2 4 3 5
j=
1 2 3 3 4 4 5 5 5 6 6 7 7
s=
1 2 -1 3 -2 4 1 -3 5 2 -4 3 -5
A1 =
1 0 -1 0 1 0 0
0 2 0 -2 0 2 0
C=
11 11 11
>> B=[3 5 7;0 1 0]
B=
357 010
>> D=[0 0 0;0 0 0]
D=
000 000
2.随机生成：（1）一个含有五个元素的列向量. （2）一个数值在 0~100 之间的三行四列的矩阵.
答 (1) rand(5,1)
ans =
0.9501
0.2311 0.6068 0.4860 0.8913
0 0 3 0 -3 0 3
0 3 0 1 0
0 0 0
4
0
4
0
0 0 3 0 1
0 0 0 0 5 0 5
答（1） n=nnz(A)
n=
11
>> nonzeros(A)
MATLAB数学实验第二版课后练习题含答案
MATLAB数学实验第二版课后练习题含答案课后练习题MATLAB数学实验第二版的课后练习题如下：第一章课后练习题1.编写MATLAB程序，计算并输出下列公式的结果：y = \\frac{1}{\\sqrt{2\\pi\\sigma^2}} e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}其中，x, $\\mu$, $\\sigma$ 分别由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');mu=input('请输入 mu 的值:');sigma=input('请输入 sigma 的值:');y=1/sqrt(2*pi*sigma^2) *exp(-(x-mu)^2/ (2*sigma^2));fprintf('y = %.2f\', y);2.编写MATLAB程序，求解下列方程的解：4x + y = 11\\\\x + 2y = 7答案如下：A= [4,1;1,2];B= [11;7];X=inv(A) *B;fprintf('x = %.2f, y = %.2f\', X(1), X(2));第二章课后练习题1.编写MATLAB程序，计算下列多项式的值：P(x) = x^4 - 2x^3 + 3x^2 - x + 1其中，x 由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');y=x^4-2*x^3+3*x^2-x+1;fprintf('P(%.2f) = %.2f\', x, y);2.编写MATLAB程序，绘制下列函数的图像：f(x) = \\begin{cases} x + 1, & x < 0 \\\\ x^2, & 0 \\leq x < 1 \\\\ 2x - 1, & x \\geq 1 \\end{cases}答案如下：x=-2:0.01:2;y1=x+1;y2=x.^2.* ((x>=0) & (x<1));y3=2*x-1;plot(x,y1,x,y2,x,y3);legend('y1 = x + 1','y2 = x^2','y3 = 2x - 1');总结本文提供了《MATLAB数学实验第二版》的部分课后练习题及其答案。
北京科技大学-MATLAB数学实验期末考试题目及答案
第1张1.figure(1)x=-pi:pi/50:pi;y=cos(2*(x.^2));plot(x,y,'r-.')figure(2)x=-pi:pi;y=cos(2*(x.^2));xi=-pi:pi/8:pi;yi1=interp1(x,y,xi,'*nearest');yi2=interp1(x,y,xi,'*linear');yi3=interp1(x,y,xi,'*spline');yi4=interp1(x,y,xi,'*cubic');plot(x,y,'ro',xi,yi1,'--',xi,yi2,'-',xi,yi3,'k.-',xi,yi4,'m:')2.function xdot=exf(t,x)u=t^3;xdot=[0 1;-5 0]*x+[0 3]'*u;clf,t0=0;tf=10;x0t=[1;2];[t,x]=ode23('exf',[t0,tf],x0t) y1=x(:,1),y2=x(:,2);plot(t,y1,'-',t,y2,'o')第2张1.第一张dsolve('D2y+4*y=t*cos(t)','y(0)=1,Dy(0)=1','t')第二张function xdot=exf(t,x)u=t*cos(t);xdot=[0 1;-4 0]*x+[0 1]'*u;第三张t0=0;tf=10;x0t=[1;1];[t,x]=ode23('exf',[t0,tf],x0t)y1=x(:,1)y2=cos(2*t) + (2*sin(2*t))./9 - cos(2*t).*(sin(3*t)/36 +sin(t)./4 - (t.*cos(3*t))./12 - (t.*cos(t))./4) +sin(2*t).*(cos(3*t)./36 + cos(t)./4 + (t.*sin(3*t))./12 + (t.*sin(t))./4)plot(x,y1,'b-',x,y2,'r-.')legend('y的图像','y的一阶导图像')2.x=0:2*pi/9:2*pi;y=(sin(2*x)).^2p3=polyfit(x,y,3);p5=polyfit(x,y,5);p7=polyfit(x,y,7);disp('三阶拟合函数'),f3=poly2str(p3,'x')disp('五阶拟合函数'),f5=poly2str(p5,'x')disp('七阶拟合函数'),f7=poly2str(p7,'x')x1=0:pi/50:2*pi;y3=polyval(p3,x1);y5=polyval(p5,x1);y7=polyval(p7,x1);plot(x,y,'rp',x1,y3,'--',x1,y5,'k-',x1,y7,'-');legend('拟合点','三阶拟合','五阶拟合','七阶拟合')第3张1.t1=-pi:pi/4:pi;x1=cos(t1);y1=sin(t1);z1=t1;subplot(1,2,1),plot3(x1,y1,z1)t2=-4:0.25:4[x2,y2]=meshgrid(t2);z2=sin(x2+y2);subplot(1,2,2),surf(x2,y2,z2),title('sin(x+y)')2.function xdot=exf(t,x)u=0;xdot=[0 1;-17/4 -1]*x+[0 0]'*u;clf,t0=0;tf=8;x0t=[-1;2];[t,x]=ode23('exf',[t0,tf],x0t)y1=x(:,1)y2=x(:,2)plot(t,y1,'b-',t,y2,'r-.')legend('y的图像','y的一阶导图像')第4张1、figure(1)[x,y]=meshgrid(-2:.5:2);z=x.^2+y.^2+3;[x1,y1]=meshgrid(-2:.1:2);z1=x1.^2+y1.^2+3;subplot(1,2,1),mesh(x,y,z),title('数据点') subplot(1,2,2),mesh(x1,y1,z1),title('函数图象') [xi,yi]=meshgrid(-2:.125:2);zi1=interp2(x,y,z,xi,yi,'*nearest');zi2=interp2(x,y,z,xi,yi,'*linear');zi3=interp2(x,y,z,xi,yi,'*spline');zi4=interp2(x,y,z,xi,yi,'*cubic');figure(2)subplot(221),mesh(xi,yi,zi1),title('最近点插值') subplot(222),mesh(xi,yi,zi2),title('线性插值') subplot(223),mesh(xi,yi,zi3),title('样条插值') subplot(224),mesh(xi,yi,zi4),title('立方插值')2、function xdot=exf(t,x)u=t;xdot=[0 1;-2,2]*x+[0 1]'*u;t0=0;tf=4;x0t=[1;pi];[t,x]=ode23('exf',[t0,tf],x0t); y1=x(:,1);y2=x(:,2);plot(t,y1,'r-',t,y2,'b-')legend('y','y一阶导数')第5张1.x=-5:0.01:0;y=3.*exp(sin(x)).*cos(x)+x.^(-2); plot(x,y,'r.')hold onx=0:0.01:5;y=3.*exp(sin(x)).*cos(x)+x.^(-2); plot(x,y,'g.')2.function xdot=exf(t,x)u=t.*tan(t);xdot=[0 1;-4 0]*x+[0 1]'*u; clf,t0=0;tf=10;x0t=[1;2];[t,x]=ode23('exf',[t0,tf],x0t)y1=x(:,1),y2=x(:,2);plot(t,y1,'-',t,y2,'o')第6张1.t=0:pi/50:2*pi;x=2*cos(t);y=2*sin(t);subplot(1,1,1),plot(x,y,'b-')2.function xdot=exf(t,x)u=1-t^2;xdot=[0 1;0 t]*x+[0 1]'*u;clft0=0;tf=2;x0t=[0.1;0.2];[t,x]=ode23('exf',[t0,tf],x0t) y1=x(:,1),y2=x(:,2);plot(t,y1,'r.',t,y2,'b-')第7张1.t1=-8:0.1:8;[x1,y1]=meshgrid(t1);z1=x1.^2/16-y1.^2/9; subplot(1,2,1),mesh(x1,y1,z1)t2=0:pi/50:2.*pi;x2=3*cos(t2);y2=2*sin(t2);subplot(1,2,2),plot(x2,y2)2.for n=1:40p0(n)=prod(100:-1:100-n+1)/100^n;p1(n)=1-p0(n);endn=1:40;plot(n,p0,n,p1,'-')xlabel('人数'),ylabel('概率')legend('房间各不同的概率','至少两人相同的概率') axis([0 40 -0.1 1.1]),grid on。
MATLAB数学实验 练习题(附答案)教学教材
%2、求100~999之间的水仙花数 clear all； clc; for n=100:999 n1=floor(n/100); %取出百位数字n1 n2=mod(floor(n/10),10); %取出十位数字n2 n3=mod(n,10) ; %取出个位数字n3 if n1^3+n2^3+n3^3==n
MATLAB数学实验 练习题(附答 案)
4.作yoz面上的曲线z=sin(y)+1绕z轴旋转所得到的图形。
5.作xoy面上的曲线(x-a)2+y2=r2绕z轴旋转所得到的图形 (0<r<a),取r=1,a=2. 6. 作x2/2-y2/3=2z图形. 7. 作x2/2+字符，若输入一个大写字母，则输出其对应的小写字母，若输入一个 小写字母，则输出其对应的大写字母;若为其它字符，则原样输出。 clear all,clc; ch=input('请输入一个字符:','s'); if ch>='A' && ch<='Z'%判断条件 ch=char(ch+32);%大写字母与小写字母ASCII值相差32 else ch=char(ch); end ch
此课件下载可自行编辑修改，仅供参考！ 感谢您的支持，我们努力做得更好！谢谢
%6. 作x2/2-y2/3=2z图形. clear,clf,clc; [X,Y,Z]=meshgrid(-1:0.1:1); V=X.^2./2+Y.^2./2-2*Z; fv=isosurface(X,Y,Z,V,0);%显示V=0等值面 p=patch(fv); set(p,'FaceColor','blue','EdgeColor','none');%设置图形的相关属性 view(3) hold on camlight
数学实验(matlab)样题及参考解答
东华大学高等数学实验试题A考试时间：90分钟（附参考解答）班级 学号 姓名 得分 上机考试说明：1. 开考前可将准备程序拷到硬盘, 开考后不允许用移动盘，也不允许上网；2. 领座考生试卷不同，开卷，可利用自己备用的书和其他资料，但不允许讨论，也不允许借用其他考生的书和资料。
3. 解答(指令行，答案等)全部用笔写在考卷上。
一、 计算题（70分）要求：写出M 函数（如果需要的话）、MATLAB 指令和计算结果。
1． 解线性方程组⎪⎪⎩⎪⎪⎨⎧-=+=+--=-+=-+14235231543421431321x x x x x x x x x x x 并求系数矩阵的行列式。
指令行：A=[5 1 –1 0;1 0 3 –1;-1 –1 0 5;0 0 2 4];b=[1;2;3;-1]; x=A\b,d=det(A) 结果：x 1=1.4, x 2= -5.9, x 3=0.1, x 4= -0.3. 行列式=70.2． 设 f(x,y) = 4 sin (x 3y)，求 3,22==∂∂∂y x y x f 。
指令行：syms x y; f=diff(4*sin(x^3*y),x); f=diff(f,y); f=subs(f,x,2); f=subs(f,y,3)结果：1063.63． 求方程 3x 4+4x 3-20x+5 = 0 的所有解。
指令行：roots([3 4 0 –20 5])结果：-1.5003 - 1.5470i, -1.5003 + 1.5470i, 1.4134, 0.25394． 使用两种方法求积分dx e x 210221-⎰π的近似值。
方法一：指令行：syms x; s=int(1/sqrt(2*pi)*exp(-x^2/2),0,1); vpa(s,5)结果：0.34135方法二：指令行：x=0:0.01:1; y=1/sqrt(2*pi)*exp(-x.^2/2);trapz(x,y)结果：0.3413方法三：M 函数ex4fun.mfunction f=ex4fun(x)f=1/sqrt(2*pi)*exp(-x.^2/2);指令行：s=quadl(@ex4fun,0,1)结果：0.34135． 求函数 f(x,y) = 3x 2+10y 2+3xy-3x +2y 在原点附近的一个极小值点和极小值。