PRINCIPAL VALUE + 3D NINTEGRATE+IMPROVED
- To: mathgroup at smc.vnet.net
- Subject: [mg81585] PRINCIPAL VALUE + 3D NINTEGRATE+IMPROVED
- From: djokic at phy.bg.ac.yu
- Date: Fri, 28 Sep 2007 02:14:55 -0400 (EDT)
A big HI to everybody !!!
A problem I am faced with is about numerically solving principal value
of a 3D integral whose integrand contains an infinity number of
singularities on a surface. So, in order not to be explaining that in
detail, here my problem is:
Q[y_] := (52.*y)/ArcTanh[(2.*y)/5.];
g[y_] := FindRoot[Q[x] == y, {x, 5., 2.}];
R[y_] := Re[g[y][[1]][[2]]];
J[x_, y_, z_] := 4.92*(Cos[x] + Cos[y] + Cos[z]) +
3.88*(Cos[(x + y)/2.] + Cos[(y + z)/2.] + Cos[(z + x)/2.]);
j[x_, y_, z_] := 3.88*(Cos[(x - y)/2.] + Cos[(y - z)/2.] + Cos[(z - x)/
2=2E]);
W[x_, y_, z_] := Sqrt[(39.942 + J[x, y, z] - j[x, y, z])*
(39.942 - J[x, y, z] - j[x, y, z])];
=CE=A9[K_, T_, =CF=89_] := ((2.*K^2*R[T]^3)/Pi^3)*
NIntegrate[(((39.942 - j[x, y, z])*(39.942 - j[x, y, z]))/W[x, y,
z])*
(1/(=CF=89^2 - 4.*R[T]^2*W[x, y, z]^2 + 0.000001))*
Coth[(R[T]*W[x, y, z])/(1.44*T)], {x, -Pi, Pi}, {y, -Pi, Pi},
{z, -Pi, Pi}];
s = Table[=CE=A9[K = 2., T, =CF=89 = 50.], {T, 1., 152., 2.}]
When obtaining the function =CE=A9[T] and drawing it depending on parameter
T, the data result achieved in such a way does not look so smooth,
particularly in the case when =CF=89 <200. Well, how to get this function
as smoothed as possible?
All the best,
Dejan Djokic
P=2E S. 0.000001 has ben put in the integration to avoid
singularities on the surface.