Calculating Multiple Integrals on Mathematica (Electric fields)
- To: mathgroup at smc.vnet.net
- Subject: [mg113049] Calculating Multiple Integrals on Mathematica (Electric fields)
- From: Hugo Almeida <squall.griffinheart at gmail.com>
- Date: Tue, 12 Oct 2010 04:25:48 -0400 (EDT)
So i'm trying to calculate a the electric field generated by 2 plaques
in the Y component. And then the median of every contribuitons from
each element of the plaques.
I have the formulas and etc, and to my best knowledge they are fine.
The problem i'm hitting is that the calculation of the median takes a
lot of time and i get this error about inumr.
ax := 0.09
az := 0.03
d := 0.052
V := 3000
E0 := 8.85*10^(-12)
Ke := 8.9*10^9
Ey1[x_, y_, 0] :=
Ke*E0*(V/d)*
Integrate[((y +
d/2)/(((x - x0)^2 + (y + d/2)^2 + z0^2)^(3/2))), {x0, -ax/2,
ax/2}, {z0, -az/2, az/2}]
Ey2[x_, y_, 0] :=
Ke*E0*(V/d)*
Integrate[-((y -
d/2)/(((x - x0)^2 + (y - d/2)^2 + z0^2)^(3/2))), {x0, -ax/2,
ax/2}, {z0, -az/2, az/2}]
Ey[x_, y_, 0] := Ey1[x, y, 0] + Ey2[x, y, 0]
Eymedio[y_] := (1/ax)*Integrate[Ey[x, y, 0], {x, -ax/2, ax/2}]
Ey1[0, 0, 0]
Ey2[0, 0, 0]
Ey[0, 0, 0]
Eymedio[0]