求助高手一个关于matlab编程的问题, 求解nd的值,非常感谢给予指导 
a=2\pi*0.41kHz; b=0.02kHz; c=1000kHz;T=1nK; E0=2.0kHz;

Z=\frac{-ib^2}{\sqrt{E^2+a^2}}\left(
                                \begin{array}{cc}
                                  E & a \\
                                  a & -E \\
                                \end{array}
                              \right)

 G0=\frac{1}{(E-E0)(E-E0-c/2)}\left(
                                \begin{array}{cc}
                                  E-E0-\frac{c}{2}+\frac{c}{2}nd & 0 \\
                                  0& -E-E0-\frac{c}{2}+\frac{c}{2}nd \\
                                \end{array}
                              \right)

G=\frac{1}{1/G0-Z}

nd=i\int\frac{dE}{2\pi}\frac{1}{e^{E/T}-1}(G-G')


求nd