4. (1)
syms x positive
assumeAlso(x,'integer')
f1=x^4-5.1*x^3+58.04*x^2-264.384*x+321.408; r1=solve(f1)
(2)
syms x clear
syms a positive
solve(x^2-a*x+a^2)
5. A=sym('[a11,a12,a13;a21,a22,a23;a31,a32,a33]')
或 A=sym('a%d%d',[3,3])
Dt=det(A)
Iv=inv(A)
[Ivs,d]=subexpr(Iv,'d')
6. syms x k
f=x^k;
s=symsum(f,k,0,inf)
subs(s,x,{sym('-1/3'),sym('1/pi'),sym('3')})
7. syms k
syms x positive
f=2/(2*k+1)*((x-1)/(x+1))^(2*k+1);
s=simple(symsum(f,k,0,inf))
8. syms t
y=abs(sin(t))
d=diff(y)
d0=limit(d,t,0,'left')
d1=limit(d,t,pi/2) (or subs(d,t,sym('pi/2'))) 9. syms x clear
syms x
f=exp(-abs(x))*abs(sin(x))
vpa(int(f,x,-5*pi,1.7*pi),64)
10. syms x y