求大神帮助啊
用MATLAB画出激光腔的稳定图（D为横坐标d为纵坐标）-2<tr(b)<2.传输矩阵及数据在下面给出。a=cosd(22);
syms d D
r=0.1;
l=0.01;
n=1.76;
b=[1/n l/2*n;0 n*n]*[1 d;0 1]*[1 0;-2/r*a 1]*[1 D;0 1]*[1 0;-2/r*a 1]*[1 d;0 1]*[n l/2*n*n;0 1/n*n];