1、首先看一个数值计算例子

%定义矩阵A
A = [ 8     1     6;
      3     5     7;
      4     9     2; ];
%获取矩阵A的特征值和特征向量
[eVec, eVal] = eig(A) 

2、符号计算例子

（1）定义一个函数f：

%定义函数
syms theta
f   = sin(theta)^2
%计算微分
df  = diff(f, theta)

（2）定义另一个符号计算例子

一个点移动的动能

syms t x(t) mass 
velocity =  diff( x(t), t);

KE       = (1/2)*mass*velocity^2
%使用符号V代替dx/dt----diff(x(t),t)，subs函数
syms v
KE_again = subs(KE, diff(x(t),t), v)

对KE进行速度微分【动能对时间的导数为功率；动能对速度的微分为动量（质量乘以速度=mv）】

%动量
der_of_KE = diff(KE_again, v)