clear

clc

close all

%% 第一題

% （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

%% 第三題

% 對(duì)sin（x）進(jìn)行多次求導(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);