Pole magnetyczne wytworzone przez obwód z prądem w kształcie “skrzywionej” pętli kołowej, która nie jest płaska !

rp={R*Cos[t],R*Sin[t],(R/2)*Sin[2*t]}

R=2;plot1=ParametricPlot3D[rp,{t,0,2*Pi},AxesLabel->{"x","y","z"},PlotStyle->{Thick}];R=.;

Show[plot1]

dlpodt=D[rp,t]

r={x,y,z}

v1=Cross[dlpodt,r-rp]

s1=Simplify[Sqrt[(r-rp).(r-rp)]]

s2=Simplify[v1/s1^3,R>0&&t>0&&t<2*Pi]

f1=s2/.{R->u*z}

f2=Simplify[Normal[Series[f1,{u,0,2}]]]

B2=Integrate[f2,{t,0,2*Pi}]

B22=Simplify[B2/.{u->R/z}]

(* Efektywnym polem powierzchni tej "pokrzywionej" ramki jest też pole powierzchni koła o promieniu R *)

Simplify[B22/(Pi*R^2)]

(* Ten wynik należy pomnożyć przez mu0*I/(4*Pi) *)

(* Teraz dla dipola elektrycznego *)

ClearAll[a,r1,r2,r];

rp1={0,0,-L/2}

rp2={0,0,L/2};

r={x,y,z};

E1=k*Q*(r-rp2)/Sqrt[(r-rp2).(r-rp2)]^3-k*Q*(r-rp1)/Sqrt[(r-rp1).(r-rp1)]^3

E2=E1/.{L->w*z}

E3=Simplify[Normal[Series[E2,{w,0,2}]]]

E4=E3/.{w->L/z}

Simplify[E4/(k*Q*L)]