elseif n==4;

t=(0:4)*pi/2+pi/2; else

t=[f(p);f(p)+pi/p]; end

0.80.60.40.20-0.2-0.4-0.6-0.8-0.8-0.6-0.4-0.200.20.40.60.81

实验十 Matlab的符号计算

21?x?2xy?y1. 展开??,它有多少项？多项式的最大系数是多少？

20最高幂次是多少？ 解：>> syms x y; >> s=(1+x+2*x*y^2+y)^20;

>> expand(s)

2. 分解因式xn?1, n=20,50,342,678。 >> syms x n; >> s=x^20-1; >> factor(s) 23. 对有理式：?4x?x?4?3x?2?x?x2?x?1?x2进行通分，因式分解，拆分，约

分。

解：通分：>> syms x; 因式分解：>> factor(s) 拆分：>> expand(s) 约分：>> collect(s) 4.

y?xarctgx，求其高阶导数y(100)及其在0点的值。

解：>> syms x; >> y=x*atan(x); >> diff(y,100) >> subs(ans,'0')

5. 计算外径为5，内径为1的圆环的体积。 解：有问题

6. 计算曲线Sin(x),0?x?2?绕x轴旋转一周所围部分的体积。

解：>> y=sym('pi*sin(x)^2'); >> int(y,0,2*pi)

7. 求解下面的不定积分和定积分：

??x?2?0x?3x?53322??x?x?1??2x?6dx dx

xcosxdx?asinx25?3sin(x)

?2ax?ax422dx23a

??/2d?0???12?13cos??/23?x?2?dx

?x?2??323解：>> y=sym('(x^3+3*x^2-5)/((x^2-2*x-6)*(x^3+x+1))'); >> int(y,x)

>> y=sym('a^x*sin(x)*cos(x^2)'); >> int(y,x)

>> y=sym('1/(5+3*sin(x))'); >> int(y,0,2*pi)

>> y=sym('sqrt(x^2-a^2)/x^4'); >> int(y,'a','2*a')

>> y=sym('1/(12+13*cos(x))'); >> int(y,0,pi/2) /5)/2) >> clear

>> y=sym('(x-2)^(2/3)/((x-2)^(2/3)+3)'); >> int(y,3,pi/2)

8. 求直线ax+by=1和双曲线x2?y2?1的交点。

解：>> [X,Y]=solve('a*x+b*y=1','x^2-y^2=1','x,y')

9. 解以下方程

5x?3x?5?0

42

?x2?y3?xy ??x?y?xy??1解：>> solve('5*x^4+3*x^2-5=0')

>> [X,Y]=solve('x^2+y^3-x*y=0','x+y+x*y-1=0')

10.求解以下微分方程