有哪位matlab高手帮我看下下面这段程序，是利用四个方程解四个未知数

syms k k1 k2 k3
s=150,d=3800,d1=1200,d2=1500,d3=1100,b=122,b1=80,b2=60,b3=70,sl=0.96,h=6,h1=12,h2=15,h3=10,l11=0.01,l12=0.01,l13=0.02,l01=0.15,l02=0.16,l03=0.12,s1=50,s2=60,s3=40
a1=0.436183,a2=-0.120167,a3=0.9372986,a=0.33267
t=1/(1+a*k);
t1=1/(1+a*k1);
t2=1/(1+a*k2);
t3=1/(1+a*k3);
f=(1/(2*pi)^(1/2))*exp(-1/2*k^2);
g=1-1/(2*pi)^(1/2)*exp(-k^2/2)*(a1*t+a2*t^2+a3*t^3);
f1=(1/(2*pi)^(1/2))*exp(-1/2*k1^2);
g1=1-1/(2*pi)^(1/2)*exp(-k1^2/2)*(a1*t1+a2*t1^2+a3*t1^3);
f2=(1/(2*pi)^(1/2))*exp(-1/2*k2^2);
g2=1-1/(2*pi)^(1/2)*exp(-k2^2/2)*(a1*t2+a2*t2^2+a3*t2^3);
f3=(1/(2*pi)^(1/2))*exp(-1/2*k3^2);
g3=1-1/(2*pi)^(1/2)*exp(-k3^2/2)*(a1*t3+a2*t3^2+a3*t3^3);
p=k*g+f-k;
q=1+k*f-g+k^2*g-k^2*g^2-f^2-2*k*f*g;
p1=k1*g1+f1-k1;
q1=1+k1*f1-g1+k1^2*g1-k1^2*g1^2-f1^2-2*k1*f1*g1;
p2=k2*g2+f2-k2;
q2=1+k2*f2-g2+k2^2*g2-k2^2*g2^2-f2^2-2*k2*f2*g2;
p3=k3*g3+f3-k3;
q3=1+k3*f3-g3+k3^2*g3-k3^2*g3^2-f3^2-2*k3

你把我的matlab憋死好几次。

clc;clear
syms k k1 k2 k3
s=150,d=3800,d1=1200,d2=1500,d3=1100,b=122,b1=80,b2=60,b3=70,sl=0.96,h=6,h1=12,h2=15,h3=10,l11=0.01,l12=0.01,l13=0.02,l01=0.15,l02=0.16,l03=0.12,s1=50,s2=60,s3=40
a1=0.436183,a2=-0.120167,a3=0.9372986,a=0.33267
t=1/(1+a*k);
t1=1/(1+a*k1);
t2=1/(1+a*k2);
t3=1/(1+a*k3);
f=(1/(2*pi)^(1/2))*exp(-1/2*k^2);
g=1-1/(2*pi)^(1/2)*exp(-k^2/2)*(a1*t+a2*t^2+a3*t^3);
f1=(1/(2*pi)^(1/2))*exp(-1/2*k1^2);
g1=1-1/(2*pi)^(1/2)*exp(-k1^2/2)*(a1*t1+a2*t1^2+a3*t1^3);
f2=(1/(2*pi)^(1/2))*exp(-1/2*k2^2);
g2=1-1/(2*pi)^(1/2)*exp(-k2^2/2)*(a1*t2+a2*t2^2+a3*t2^3);
f3=(1/(2*pi)^(1/2))*exp(-1/2*k3^2);
g3=1-1/(2*pi)^(1/2)*exp(-k3^2/2)*(a1*t3+a2*t3^2+a3*t3^3);
p=k*g+f-k;
q=1+k*f-g+k^2*g-k^2*g^2-f^2-2*k*f*g;
p1=k1*g1+f1-k1;
q1=1+k1*f1-g1+k1^2*g1-k1^2*g1^2-f1^2-2*k1*f1*g1;
p2=k2*g2+f2-k2;
q2=1+k2*f2-g2+k2^2*g2-k2^2*g2^2-f2^2-2*k2*f2*g2;