1、实验三：1,x=input(输入X的值：);if x=0y=input(请输入工作时长：);if y120a=84*120+(y-120)*84*0.15;disp(您本月工资是： ,num2str(a)elseif y=10-5)x1=a/(b+x0);x0=x1;n=n+1;endnx04、clear allfor n=1:4if n=1f1=1;elseif n=2f2=0;elseif n=3f3=1;elsea=f3-2*f2+f1;b=a-2*f3+f2;c=b-2*a+f3;d=c-2*b+a;H=1,0,1,a,b,c,d;for m=8:4:99a=d-2*c+b;b=a-2
2、*d+c;c=b-2*a+d;d=c-2*b+a;H=H,a,b,c,d;endf100=d-2*c+b;endendmax=max(H);min=min(H);sum=sum(H);disp(最大值是：,num2str(max)disp(最小值是：,num2str(min)disp(各数和是：,num2str(sum)k=0;l=0;p=0;for e=Hif e0k=k+1;elseif e0.5n=n+1;else continue;endendp=n/30000;disp(均值= ,num2str(meanA);disp(标准差=,num2str(stdA);disp(最大值max=,
3、num2str(maxA);disp(最小值min=,num2str(minA);disp( 百分比P=,num2str(p);实验九：2,(1).x=0:pi/100:2*pi;y=sqrt(cos(x.2)+4*sin(2*x).2+1);trapz(x,y)g=inline(sqrt(cos(x.2)+4*sin(2*x).2+1);i=quadl(g,0,2*pi)j=quad(g,0,2*pi)format short2,(2).x=0:0.001:1;y=log(1+x)./(1+x.2);trapz(x,y)g=inline(log(1+x)./(1+x.2);i=quadl(g
4、,0,1)j=quad(g,0,1)format short3, a=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2;b=-4;13;1;11;x1=abl,u=lu(a);x2=u(lb)q,r=qr(a);x3=r(qb)format short5，function f=fz(x);f=3*x+sin(x)-exp(x);6，（1）：f=inline(x3+cos(x)+x*log(x)/exp(x);fminbnd(f,0,1)（2）：function f=fxy(u)x=u(1);y=u(2);f=2*x.3+4*x.*y.2-10*x.*y+y.2;实验十：
5、1, x=sym(6);y=sym(5);z=(x+1)/(sqrt(3+x)-sqrt(y)eval(z)2, syms x y z=x4-y4factor(z)factor(sym(5135) %5135整数分解3, syms x ym=sin(x)*cos(y)-cos(x)*sin(y)n=(4*x2+8*x+3)/(2*x+1)simple(m) %显示过程simplify(m) simplify(n) %最优4, p1=sym(0,1,0;1,0,0;0,0,1)p2=sym(1,0,0;0,1,0;1,0,1)A=sym(a,b,c;d,e,f;g,h,i)B=p1*p2*A C
6、=inv(B) %逆矩阵B.*C-C.*B %验证逆矩阵的正确性triu(B) %上三角阵tril(B) %下三角阵diag(B,1) %主对角线det(B)5, syms x a t ;A=(x*(exp(sin(x)+1)-2*(exp(tan(x)-1)/sin(x)3B=(sqrt(pi)-sqrt(acos(x)/sqrt(x+1)C=(1-cos(2*x)/xD=sym(ax,t3;t*cos(x),log(x)E=det(D)F=(x2-2*x)*exp(-x2-y2-x*y)a=limit(A)b=limit(A,0)c=limit(A,x,0)d=limit(B,x,-1,r
7、ight)e=diff(C) %一阶导数f=diff(C,2) %二阶导数g=diff(E,x) %对x求偏导h=diff(E,t,2) %对t求偏二阶导i=diff(diff(E,x),t) %d(de/dx)/dtj=-diff(F,x)/diff(F,y)%对y求x的一阶偏导数k=diff(diff(F,x),y)%对f求二阶偏导数x=0;y=1;eval(j)eval(k)6.x=sym(x);f=1/(1+x.4+x.8);int(f)f=1/asin(x).2.*sqrt(1-x.2);int(f)f=(x.2+1)/(x.4+1);int(f,x,0,inf)f=exp(x).*(1+exp(x).2;int(f,x,0,log(2)
