源代碼見文末
部分源代碼：
1 % 符號變量 兩種表達方式 2 a=sym('a'); 3 class(a); 4
5 syms b; 6 b; 7
8 % 符號常量 9 c=sym('3'); 10 c; 11 % 符號表達式 三種表達方式 12 f1='3*x+4'; % char類型，但是可以參與運算。所以這種形式是可以直接用的，比較方便 13
14 syms x 15 f2=3*x+4; % sym類型 16
17 f3=sym('3*y+4'); % sym類型 18
19 % 符號四則運算 20 fadd1 = f1 + f2; 21 fmu1 = f1*f2; 22
23 % 符號表達式的化解 24 syms x y 25 s = (x^2+y^2)^2 + (x^2-y^2)^2; 26 simplify(s); 27
28 % 符號表達式和數值的轉換 29 eval(c); 30 syms; % 上面已經涉及到 31
32 % 因式分解、展開與合並同類項 33
34 syms a b x y 35 f1=a^3-b^3; 36 factor(f1); 37
38 f2=(3*x^2+8*y^2)*(-x^2+3*y); 39 expand(f2) 40
41 f3 = 2*x+3*x+4*x+y; 42 collect(f3) 43
44 % 符號矩陣 45 a1=[x x+y;y x+y]; 46 a1 47 transpose(a1) % 轉置 48
49 % 符號函數值的求解 50 syms x 51 f1 = x^3 -9; 52 subs(f1,3) 53
54 % 符號極限、符號微分、符號積分 55 syms x 56 y=sin(x); 57 limit(y,0) 58
59 y1 = x^2; 60 diff(y1) % 求導、微分 61 diff(y1,2) % 二階導數 62
63 y3=(3-x^2)^3; 64 int(y3) % 不定積分 65
66 y4=abs(1-x); 67 int(y4,1,2) % 定積分1-2范圍內 68
69 %%符號方程求解 70 %代數方程 71 clear 72 syms x 73 solve(x+x*exp(x)-10); 74
75 solve('x+x*exp(x)-10') % 和前面的符號表達式方法聯系起來了 76 eval(solve('x+x*exp(x)-10')) % 這是數值的表達結果 77
78 %方程組 79 [x,y] = solve('x+y=10','x-y=2','x,y') % 輸出x和y的值 80
81 % 符號的微分方程 82 dsolve('Dy-(x^2+y^2)/x^2/2','x') % 微分方程組合上面的方程組道理一樣