第7章例题源程序..docVIP

例7.1 当λ取何值时，齐次线性方程组： 有非零解。 syms lamda A=[1-lamda,-2,4;2,3-lamda,1;1,1,1-lamda]; D=det(A); factor(D) 例7.2 求下列极限。 （1） （2） （3） （4） syms a m x; f=(x^(1/m)-a^(1/m))/(x-a); limit(f,x,a) %求极限(1) f=(sin(a+x)-sin(a-x))/x; limit(f) %求极限(2) f=x*(sqrt(x^2+1)-x); limit(f,x,inf,'left') %求极限(3) f=(sqrt(x)-sqrt(a)-sqrt(x-a))/sqrt(x*x-a*a); limit(f,x,a,'right') %求极限(4) 例7.3 求下列函数的导数。 (1) ，求y' 。 (2)y=xcos(x)，求y''、y'''。 (3)，求、。 (4)，求、。 (5)z=f(x,y)由方程x2+y2+z2=a2定义，求、。 syms a b t x y z; f=sqrt(1+exp(x)); diff(f) %求(1)。未指定求导变量和阶数，按默认规则处理 f=x*cos(x); diff(f,x,2) %求(2)。求f对x的二阶导数 diff(f,x,3) %求(2)。求f对x的三阶导数 f1=a*cos(t);f2=b*sin(t); diff(f2)/diff(f1) %求(3)。按参数方程求导公式求y对x的导数 %求(3)。求y对x的二阶导数 (diff(f1)*diff(f2,2)-diff(f1,2)*diff(f2))/(diff(f1))^3 f=x*exp(y)/y^2; diff(f,x) %求(4)。z对x的偏导数 diff(f,y) %求(4)。z对y的偏导数 f=x^2+y^2+z^2-a^2; zx=diff(f,x)/diff(f,z) %求(5)。 zy=diff(f,y)/diff(f,z) %求(5)。 例7.4 在曲线y=x3+3x-2上哪一点的切线与直线y=4x-1平行。 x=sym('x'); y=x^3+3*x-2; %定义曲线函数 f=diff(y); %对曲线求导数 g=f-4; solve(g) %求方程f-4=0的根，即求曲线何处的导数为4 例7.5 求下列不定积分。 (1) (2) (3） (4) x=sym('x'); f=(3-x^2)^3; int(f) %求(1) f=sin(x)^2; -1/2*cos(x)*sin(x)+1/2*x syms alpha t; f=exp(alpha*t); int(f) %求(3) f=5*x*t/(1+x^2); int(f,t) %求(4) 例7.6 求下列定积分。 (1) (2) (3) (4) x=sym('x');t=sym('t'); int(abs(1-x),1,2) %求(1) f=1/(1+x^2); int(f,-inf,inf) %求(2) f=x^3/(x-1)^10; I=int(f,2,3) %求(3) double(I) %将上述符号结果转换为数值 int(4*x/t,t,2,sin(x)) %求(4) 例7.7 求椭球的体积。 syms a b c z; f=pi*a*b*(c^2-z^2)/c^2; V=int(f,z,-c,c) 例7.8 求空间曲线c从点(0，0，0)到点(3，3，2)的长度。设曲线c的方程是：