应用背景：求解光在单轴晶体中，各传播方向的光波与折射率的关系
将（48）、（49）代入光的波法线菲涅尔方程
得到
利用Matlab求解符号代数方程求解（50）
>> syms n n0 ne a
y=(n^2-n0^2)*(n^2*(n0^2*(sin(a))^2+ne^2*(cos(a))^2)-n0^2*ne^2)
>> solve(y,n)
ans =
n0
-n0
(n0*ne)/(n0^2*sin(a)^2 + ne^2*cos(a)^2)^(1/2)
-(n0*ne)/(n0^2*sin(a)^2 + ne^2*cos(a)^2)^(1/2)