Dear Fairmat,

I tried to write a plug-in and have some questions. I saw that you have Rectangle, AdaptLobatto, AdaptLobatto, and Romberg integration. However, I would like to use adaptive Simpson quadrature integration. If you have this integration, please give me some examples. I have the functions in Matlab that want to write in fairmat plug-in as

Matlab
F=@(x)(x.^(l-1)).*((1-x).^(n-l-1)).*((1-x.*o).^(-h)).*exp(x.*p);
quad(F,0,1);

This integration will have a problem with singularity. Moreover, I tried adaptive Gauss-Kronrod quadrature (quadgk and adaptiveSimpson), adaptive Lobatto quadrature (quadl and adaptiveLobatto), Gauss-Legendre quadrature (quadgr), some type of Gaussian quadrature (quadg), cumulative integral of adaptive Simpson quadrature (cumquad) in Matlab and they are not as good as adaptive Simpson quadrature integration. I wonder your integration will work same as Matlab or not. If you have any suggestion, please give me some examples in adaptive Simpson quadrature. Or if you don't have adaptive Simpson quadrature, could you show me in AdaptLobatto.

Best Regards,
Khampol