%% 第一題 % （1） x = -2:0.001:2; y = x.*sin(1./x); plot(x,y); xlabel('x') ylabel('y') title('plot畫的圖像'); figure fplot(@(z)z*sin(1./z),[-2 2]); xlabel('x') ylabel('y') title('fplot畫的圖像'); %(2) syms m; f=m*sin(1/m); disp('xsin(1/x)的極限為：'); L=limit(f, m, 0) %% 第二題 A = [4 2 -1 3 -1 2 12 3 0]; B = [2;10;8]; disp('方程組的解為：'); x = inv(A)*B %% 第三題 % 對sin（x）進行多次求導(dǎo)，并畫圖查看 syms x1 y = sin(x1); y1 = diff(y,x1); y2 = diff(y1,x1); y3 = diff(y2,x1); y4 = diff(y3,x1); y5 = diff(y4,x1); y6 = diff(y5,x1); y7 = diff(y6,x1); y8 = diff(y7,x1); figure subplot(411) ezplot(y,[-4*pi :0.01:4*pi]); subplot(412) ezplot(y1,[-4*pi :0.01:4*pi]); subplot(413) ezplot(y2,[-4*pi :0.01:4*pi]); subplot(414) ezplot(y3,[-4*pi :0.01:4*pi]); figure subplot(411) ezplot(y4,[-4*pi :0.01:4*pi]); subplot(412) ezplot(y6,[-4*pi :0.01:4*pi]); subplot(413) ezplot(y7,[-4*pi :0.01:4*pi]); subplot(414) ezplot(y8,[-4*pi :0.01:4*pi]); %% 第五題 A = [1 0 0 -1 0 0 0 1 0 0 -1 0 0 0 1 0 0 -1]; intcon = [1 2 3 ]; B = [200;-100;150;-80;100;-50]; f = [200 150 180]; [x fval exitflag] = intlinprog(f,intcon,A,B); disp('產(chǎn)品A B C分別產(chǎn)量為：'); disp(x); disp('總利潤為：'); disp(fval);