例 1:计算 sin(x)/(x^2+4x+3) 的四阶导数
f=sin(x)/(x^2+4*x+3); y=diff(f,x,4)
例 2:计算下面的多项式方程的根
s^6+9*s^5+135/4*s^4+135/2*s^3+1215/16*s^2+729/16*s+64=0
p=[1,9,135/4,135/2,1215/16,729/16,64];
roots(p)%数值解
p1=poly2sym(p);
solve(p1)%解析解
%对解析解和数值解的理解
%y=2x
%y=0.5x——解析解
%x=1时,y=0.5——数值解
例 3:考虑如下 Hilbert 矩阵
H=sym(hilb(20));
det(H)
例 4:计算 Van der Pol 微分方程 (不太明白)
mu=1000;
f=@(t,x)[x(2);-mu*(x(1)^2-1)*x(2)-x(1)];
[t,x]=ode15s(f,[0,3000],[-1,1]);
plot(t,x)
例 5:最优化问题
&#x