I'm currently studing the behaviour of some integral for which I cant find an analytic solution.

As this function depend on a vector on the unit sphere, I first defined it with 3 parameters before making a second one taking eulers angle and calling the first one.

However, when I try to plot my function I don't get my result. Either :

  • I get my first graph, and the worksheet resets

  • Mathematica exits

I'm using on linux (Archlinux 3.13.6-1 x86_64)

Here is my worksheet

w[a_, b_, c_] := 
 NIntegrate[(a x + b y + c) (x^2 + y^2 + 1)^(-2)*
   HeavisideTheta[a x + b y + c], {x, -1, 1}, {y, -1, 1}]
w3D[theta_, phi_] := 
 w[Sin[theta] Cos[phi], Sin[theta] Sin[phi], Cos[theta]]
w[0, 0, 1] + 4 w[0, 1, 0]

Plot[w3D[theta, 0], {theta, 0, Pi}]
Plot[w3D[theta, Pi/4], {theta, 0, Pi}]
Plot3D[w3D[theta, phi], {theta, 0, Pi}, {phi, 0, Pi/2}]

Edit: adding

PrecisionGoal -> 10, AccuracyGoal -> 5

to my Nintegrate solves the problem

  • $\begingroup$ On win7-64bit running 9.0.1 all 3 plots appear (after some minutes). Neither the kernel nor the frontend appear to quit. When you say quit what do you mean? $\endgroup$ – Ymareth Mar 19 '14 at 16:33
  • $\begingroup$ I mean the application vanish and I can't find it using ps -A It's strange since I happened to see my computer run for long minutes during complexe calculation. $\endgroup$ – Amxx Mar 19 '14 at 16:38
  • 2
    $\begingroup$ I can reproduce the kernel crash on OS X 10.9.1, but it's not consistent. Sometimes it crashes sometimes it doesn't. You should report it to WRI. $\endgroup$ – Szabolcs Mar 19 '14 at 16:57
  • 1
    $\begingroup$ Intermittent crashes on 9.01 / Win7 64bit here as well. $\endgroup$ – Yves Klett Mar 19 '14 at 17:07
  • $\begingroup$ Can someone please give me a link to the place where I could report this ? $\endgroup$ – Amxx Mar 19 '14 at 17:40

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.