matlab结果中有虚数,Matlab数值积分结果是虚数

小木虫的各位童鞋们，我现在调试matlab程序遇到一个问题。一个很复杂的关于角度theta的数值积分，我调用自带的quadl函数，结果出现虚数。我后来又根据变步长的simpson算法编了一个数值积分程序，结果还是有虚数，而且我的积分程序和quadl算出的结果还有一些误差。积分函数包含了很多三角函数，还有e的指数函数在里面，但是积分过程中并无傅里叶变换等可以带来虚数的过程，不知道虚数是怎么来的。我把积分函数曲线都画出来了，波动比较大，但是作为一个面积，怎么会出现虚数呢？我把我的八个积分函数表达式以及各自的函数曲线贴出来，希望大神们帮我看看，给点意见。论文卡在这里，一直得不到解决。

%========function group of theta=====================

function F1=Ff1(theta)

%======== basic formula=============================

% TO identify whether the rotating axis is inside the circle

if rrat<=0

rm=0.5*(exp((theta-theta0)*tan(phi))-rrat*exp((theta-theta0)*tan(phi)));

R=0.5*(exp((theta-theta0)*tan(phi))+rrat*exp((theta-theta0)*tan(phi)));

else

rm=0.5*(exp((theta-theta0)*tan(phi))+rrat*exp(-(theta-theta0)*tan(phi)));

R=0.5*(exp((theta-theta0)*tan(phi))-rrat*exp(-(theta-theta0)*tan(phi)));

end

a=sin(theta0)/sin(theta)-rm;

F1=2*cos(theta)*(((R^2-a^2)^(1/2))*(0.125*R^2*a-0.25*a^3-0.5*a*rm^2+(R^2-a^2)*rm)...

+(0.5*pi-asin(a/R))*(0.125*R^4+0.5*rm^2*R^2));

end

function F2=Ff2(theta)

%========basic formula==&#