hi, how can i accomplish this densityplot?
- To: mathgroup at smc.vnet.net
- Subject: [mg107392] hi, how can i accomplish this densityplot?
- From: çåå <hiphopphysics at yahoo.com.tw>
- Date: Thu, 11 Feb 2010 06:53:44 -0500 (EST)
f[x_, y_] :=
NIntegrate[
r0*Exp[-(r0 +
2*(x^2 + y^2 + r0^2 +
2*r0*((x^2 + y^2)^0.5)*(Cos[theta0]*(y/(x^2 + y^2)^0.5) -
Sin[theta0]*(x/(x^2 + y^2)^0.5)))^0.5)], {r0, 0,
Infinity}, {theta0, 0, Pi}]
DensityPlot[f[x_, y_], {x, -2, 2}, {y, -2, 2}]
i entered the code above but mathematica return the error message:
NIntegrate::inumr: The integrand E^(-r0-2 (r0^2+<<1>>+<<1>>+2 r0
Power[<<2>>] Plus[<<2>>])^0.5) r0 has evaluated to non-numerical
values for all sampling points in the region with boundaries {{\
[Infinity],0.},{0,3.14159}}. >>
General::stop: Further output of NIntegrate::inumr will be suppressed
during this calculation.
what's wrong with my code? please help me! thanks a lot!