4d integral

--===============2013263128058904528==
Content-Type: multipart/alternative; boundary="000000000000607532062b89de88"

--000000000000607532062b89de88
Content-Type: text/plain; charset="UTF-8"

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.

https://feynarts.de/cuba/
Non parallel version doesn't crash, 2.1 I think (Cuhre method)

--000000000000607532062b89de88
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

<div dir=3D"auto"><div>Just for the record, when calling colloc, you specif=
y number of points within the interval, for example 60, like this,</div><di=
v dir=3D"auto">[t,~,~,q]=3Dcolloc(60)<span style=3D"color:rgb(0,0,0);font-f=
amily:arial,sans-serif;font-size:14.6667px">;</span><br><br>I&#39;m using C=
uba C integrator (n dimensions) within octave via oct. If some octave dev c=
ould be interested, I could provide sample code. It would be great having n=
dim integrator (cubature) in octave.</div><div dir=3D"auto"><br></div><div =
dir=3D"auto"><a href=3D"https://feynarts.de/cuba/">https://feynarts.de/cuba=
/</a></div><div dir=3D"auto">Non parallel version doesn&#39;t crash, 2.1 I =
think (<span style=3D"background-color:rgb(238,238,255);color:rgb(0,0,0);fo=
nt-family:cmr,times,serif;font-size:14px">Cuhre method)</span></div></div>

--000000000000607532062b89de88--


--===============2013263128058904528==
Content-Type: text/plain; charset="utf-8"
MIME-Version: 1.0
Content-Transfer-Encoding: base64
Content-Disposition: inline

Ci0tLS0tLS0tLS0KV2UgYXJlIHRyYW5zaXRpb25pbmcgdG8gYSB3ZWIgYmFzZWQgZm9ydW0KZm9y
IGNvbW11bml0eSBoZWxwIGRpc2N1c3Npb25zIGF0Cmh0dHBzOi8vb2N0YXZlLmRpc2NvdXJzZS5n
cm91cC9jL2hlbHAK

--===============2013263128058904528==--