请教：同一积分表达式采用不同积分函数结果相差巨大

chunshui2003

积分表达式如图所示

分别采用quadgk    quadl    quad   来进行计算，具体参数如下：

1. 
   
2. 
   
3. %%%% 电磁刚度参数数值 %%%%%
   
4. delta_0 = 18*10^-3;        
   
5. Rg = 6.24;     
   
6. L1 = 2.1;       %转子有效长度m
   
7. k_u = 1.102;     %饱和度
   
8. mu_0 = 4*pi*10^(-7);     %空气导磁系数
   
9. f = 50;         
   
10. theta = 30.64/180*pi;      
    
11. phi = acos(0.875);         %功率因数角
    
12. p = 40;                %磁极对数
    
13. Fsm = 19210;           
    
14. Fjm = 24214;         
    
15. Lambda_0 = mu_0./(k_u*delta_0);      %发电机均匀气隙磁导
    
16. sigma = k_u*delta_0;         %饱和度与均匀气隙大小的乘积
    
17. omega_f = 2*pi*f./p;        %发电机同步转速
    
18. 
    

复制代码

积分表达式如下：

1. 
   
2. 
   
3. coeff_1 = Rg.*L1.*Lambda_0./(2*sigma.^2);
   
4. 
   
5. B = @(alpha,t) Fsm.*cos(omega_f.*t-p.*alpha);  
   
6. 
   
7. C = @(alpha,t) Fjm.*cos(omega_f.*t-p.*alpha+theta+phi+pi/2);
   
8. 
   
9. C11 = @(alpha) 1+cos(2.*alpha);
   
10. 
    
11. K11_exp2 = @(alpha,t) coeff_1.*C11(alpha).*(B(alpha,t)+C(alpha,t)).^2;
    
12. 
    
13. K11_quadgk = quadgk(@(alpha) K11_exp2(alpha,0),0,2*pi);
    
14. 
    
15. K11_quadl = quadl(@(alpha) K11_exp2(alpha,0),0,2*pi);
    
16. 
    
17. K11_quad = quad(@(alpha) K11_exp2(alpha,0),0,2*pi);
    
18. 
    

复制代码

计算结果分别为

5.0701e^8    3.1042e^8    3.8713e^8

看了一些关于积分函数的介绍，其中quadgk几乎是一个万能的而且精度较高。
三个函数得到三个不同的结果，请问根据这个积分表达式哪一个积分函数得出的结果更加可信呢？