有朋友解释一下这几行MATLAB代码不？
syms x Q Qm jieguo1 jieguo2;
f=1/(4.74*sqrt(2*pi)).*exp(-1/2*((x-11.73)/4.74).^2);
eq=int(f,x,-inf,Q)-0.77;
Q=solve(eq,Q);
VQ=vpa(Q);
Q=fix(VQ);
disp(Q);
eq1=int(f,x,-inf,Qm)-0.99;
Qm=solve(eq1,Qm);
VQm=vpa(Qm);
Qm=fix(VQm);
disp(Qm);
m=fix(Q);
n=fix(Qm);h=1:1:n+1;f1=1:1:n+1;h(1)=int(f,-inf,0);f1(1)=h(1);for i=2:n+1 h(i)=int(f,-inf,i); f1(i)=h(i)-h(i-1);endjieguo1=0;jieguo2=0;for i=1:n+1 jieguo1=jieguo1+15*(m-i+1)*f1(i);endfor i=1:m+2 jieguo2=jieguo2+15*(m+1-i+1)*f1(i);endfor i=m+1:n+1 jieguo1=jieguo1+50*(i-1-m)*f1(i);endfor i=m+2:n+1 jieguo2=jieguo2+50*(i-1-m-1)*f1(i);end
jieguo1=vpa(jieguo1);
disp(fix(jieguo1))
jieguo2=vpa(jieguo2);
disp(fix(jieguo2))