parallel numerical integration (quad, quadl)

Hey I stumble upon a mini project that is to parallelize a mathematica-serial numerical integration program to matlab. The reason is mathematica is using NIntegrate and this takes 8 hours to run on a core-i7 computer.

I found some good tutorial about using spmd to estimate integration using area-under-the curve algorithm. However, I need high degree of accuracy provided by quadl.

Would you mind point in the right direction to work out parallel quadl/int/quad in matlab? Many thanks!