我的程序是：
syms theti thets
bc=1;
M=10;
l=0.2;
h=1/1000000;
   thi=theti.*3.1415926/180;
   sini=sin(thi);
   cosi=cos(thi);
   ths=thets.*3.1415926/180;
   sins=sin(ths);
   coss=cos(ths);
   k=2*3.1415926/bc;
   a=(1+cos(thi+ths)).^2;
   b=k^3./cosi;
   c=exp(-(cosi+coss).^2.*(k^2*h^2));
   det=0;
   out=1;
for m=1:M
    out=out*m;
    e=exp(-(sins-sini).^2.*(k^2*l^2)/(4*m));
    g=(((cosi+coss).*k).^(2*m-2)/out).*(h^(2*m)*l/(2*sqrt(m*3.1415926))).*e;
    det=det+g;
end
   det=det.*a.*b.*c（前面的程序不用怎么看）
到这里为止求出det是theti和thets的函数
以下就不会了：
让theti[1:0.2:89]（角度变化）
当theti 每取一个值时 比如：theti=2时，然后对thets从[0,pi/2]积分
画出theti和det的图形
俺在线期待各位高手的帮忙 qq也行：502513092