源代码见文末
部分源代码：
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') % 微分方程组合上面的方程组道理一样