第一題

function y=mat1(x)    %建立函數(shù)文件mat1.m

y=[exp(x),log(x),sin(x),cos(x)];

%在命令窗口調(diào)用上述函數(shù)文件：

y=mat1(1+i)

%第二題程序一

function [a,b,N,M]=shiyanwu2(m,n,t)

A=[m*cos(t*pi/180),-m,-sin(t*pi/180),0;m*sin(t*pi/180),0,cos(t*pi/180),0;0,n,-sin(t*pi/180),0;0,0,-cos(t*pi/180),1];

B=[0,9.8*m,0,9.8*n];

C=inv(A)*B';

a=C(1);

b=C(2);

N=C(3);

M=C(4);

%在命令窗口調(diào)用該函數(shù)文件:

m1=input('m1=');

m2=input('m2=');

theta=input('theta=');

[a1,a2,N1,N2]=shiyanwu2(m1,m2,theta)

%第二題程序二

function X=mat2(m1,m2,t)

g=9.8;

A=[m1*cos(t*pi/180),-m1,-sin(t*pi/180),0;m1*sin(t*pi/180),0,cos(t*pi/180),0;0,m2,-sin(t*pi/180),0;0,0,-cos(t*pi/180),1];

B=[0;m1*g;0;m2*g];

X=inv(A)*B;

%在命令窗口調(diào)用該函數(shù)文件：

X=mat2(1,1,60)

%第三題

function flag=mat3(x)

flag=1;

for i=2:sqrt(x)

if rem(x,i)==0

    flag=0;

    break;

end

end

%在命令窗口調(diào)用該函數(shù)文件：

for i=10:99

j=10*rem(i,10)+fix(i/10);

if mat3(i)&mat3(j)

   disp(i)

end

end

%第四題

function y=fx(x)

y=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);

%在命令窗口調(diào)用該函數(shù)文件：

y=fx(2)

a=[1,2;3,4];

y=fx(a)

%第五題

%(1)

function f1=mat5(n)

f1=n+10*log(n*n+5);

%在命令窗口中調(diào)用該函數(shù)文件：

y=mat5(40)/(mat5(30)+mat5(20))

%(2)方法一

function f2=mat6(n)

f2=0;

for i=1:n

    f2=f2+i*(i+1);

end

%在命令窗口中調(diào)用該函數(shù)文件如：

y=mat6(40)/(mat6(30)+mat6(20))

%(2)方法二

function f2=mat7(n)

i=1:n;

m=i.*(i+1);

f2=sum(m);

end

%在命令窗口中調(diào)用該函數(shù)文件如：

y=mat7(40)/(mat7(30)+mat7(20))