上机实验报告
班级:机械工程及自动化02
姓名:鄢文光
学号:40902010221
日期:2012年11月13日
题目1-1:>> syms x;limit(cos(sqrt(x))^(pi/x),x,0,'right')
ans =
exp(-1/2*pi)
题目1-2:>> syms x;limit((3*sin(x)+(x^2)*cos(1/x))/((1+cos(x))*log(1+x)),x,0)
ans =
3/2
题目1-3:>> syms x;limit((sqrt(4*x^2+x-1)+x+1)/sqrt(x^2+sin(x)),x,inf,'left')
ans =
3
题目2-1:>> syms x;f=(tan(sqrt(x+sqrt(x+sqrt(2*x)))))^2;diff(f)
ans =
tan((x+(x+2^(1/2)*x^(1/2))^(1/2))^(1/2))*(1+tan((x+(x+2^(1/2)*x^(1/2))^(1/2))^(1/2))^2)/( x+(x+2^(1/2)*x^(1/2))^(1/2))^(1/2)*(1+1/2/(x+2^(1/2)*x^(1/2))^(1/2)*(1+1/2*2^(1/2)/x^(1/ 2)))
题目2-2:>> syms x;f=(cos(x^2))*((sin(1/x)^3));diff(f)
ans =
-2*sin(x^2)*x*sin(1/x)^3-3*cos(x^2)*sin(1/x)^2*cos(1/x)/x^2
题目3-1:>> syms x;f=sqrt(sin(x)-((sin(x))^3));int(f,0,pi)
ans =
4/3
题目3-2:>> syms x;f=(1/x)*sqrt((x+1)/(x-1)),int(f)
f =
1/x*((1+x)/(x-1))^(1/2)
ans =
((1+x)/(x-1))^(1/2)*(x-1)/((1+x)*(x-1))^(1/2)*(-atan(1/(x^2-1)^(1/2))+log(x+(x^2-1)^(1/2) ))
题目4:>> dsolve('(x^2)*Dy+x*y=y^2','y(1)=1','x')
ans =
2*x/(1+x^2)
题目5:
function ydot=fun1(t,y)
ydot=[y(2);y(3);-4*y(1)-3*y(2)-11*y(3)+3];
%建立函数fun1.m
>> y0=[1,1,1];[t,y]=ode45('fun1',[0:5:100],y0)
t =