Hello, I am new with Mathematica.
I want to define a fonction
cem[mr_] := ro[te,mr] / (Norm[ro[te,mr]]^3)
but the parameter te
is calculating by means of
te == tr - Norm[ro[te,mr]]
and it doesn't work. Could you know how to calculate cem[mr_]
with an implicite calculus of te
please ?
Here is the code :
me[t_] := {t*(1 - Erf[t]) - 1/(Sqrt[Pi]*E^t^2), 0, 0}
mr := {x, y, z}
ro[te_, mr_] := mr - me[te]
z = 0;
tr = -5;
(*
Clear[ten]; RecurrenceTable[{ten[n+1] == N[tr-Norm[ro[ten[n], mr]],100], ten[0] == tr}, ten, {n, 0, 12}]
Clear[tel]; tel[ten_]:= N[tr-Norm[ro[ten , mr]],100]; te = NestList[tel,tr,12][[13]]
*)
te == tr - Norm[ro[te, mr]];
cem[mr_] := ro[te, mr]/(Norm[ro[te, mr]]^3)
StreamPlot[cem[mr], {x, -3, 3}, {y, -3, 3}]
Thank you.
Armen SARIAN