Just for the record, when calling colloc, you specify number of points within the interval, for example 60, like this,
[t,~,~,q]=colloc(60);
I'm using Cuba C integrator (n dimensions) within octave via oct. If some octave dev could be interested, I could provide sample code. It would be great having ndim integrator (cubature) in octave.
Non parallel version doesn't crash, 2.1 I think (Cuhre method)