受限于设备的原因，计算时间较长

代码如下，请电脑好的帮忙计算一下结果
syms n2 Thta;
eta_p=1;%量子效率
nt=0.8;%掺杂浓度
Nho=3*4.55*6.02e23*nt/59370;%格位浓度
syms r z;
lr=5;%晶体长度
sigma_abs=1.2e-20;%吸收截面
arf=Nho*sigma_abs;
eta_arf=1-exp(-arf*lr);%吸收效率
h=6.626e-34;%普朗克常数
c=3e10;n=1.82;Pin=500;lamda_p=1908e-7;lamda_l=2090e-7;%光速单位cm折射率，输入功率，泵浦波长，激光波长
v_p=c/lamda_p;v_l=c/lamda_l;%泵浦和激光的频率
R=eta_p*eta_arf*Pin*(1-n2/Nho)/(h*v_p);%考虑基态损耗下的抽运速率
M2=3;%泵浦光的光束质量
wp0=200e-4;%束腰半径
z0=lr/2;%假设束腰在晶体中心
wp=wp0*sqrt(1+((z-z0)*M2*lamda_p/(pi*wp0^2))^2);
rp=2*arf*exp(-arf*z)*exp(-2*r^2/wp^2)./(pi*wp^2*(1-exp(-arf*lr)));%泵浦光归一化分布
wl=200e-4;
lc=20;%物理腔长
lcc=lc+(n-1)*lr;%腔的光学长度
j=2*exp(-2*r^2/wl^2)/(pi*wl^2*lcc);%激光的归一化分布
detaa=0.01;%腔的固有损耗
T=0.5;%输出镜透过率
deta=detaa+T;
tao_c=2*lcc/(c*deta);%腔寿命
tao2=8e-3;%上能级寿命
k=1.5e-18;sigma_em=1e-20;%发射截面
fu=0.0994;fl=0.017;%上下能级的玻尔兹曼集居数
f1=R*rp-k*n2^2-n2/tao2-sigma_em*c/n*((fu+fl)*n2-fl*Nho)*Thta*j;
N2=solve(f1==0,n2)
J=int(int((fu*N2-fl*(Nho-N2))*Thta*j,z,0,lr),r,0,0.5);
f2=sigma_em*c/n*2*pi*J-Thta/tao_c;
Th=solve(f2==0,Thta)

设备不行