Ringkasan

Mathematica 11.0 code

rs[t_] := {-(E^(-2500 t^2)/(150 Sqrt[\[Pi]])) + t/3 - 1/3 t Erf[50 t],
    0, 0};
c = 1;
q = -1;
\[Beta][t_] := {(1 - Erf[50 t])/3, 0, 0};
d\[Beta][t_] := Chop@Evaluate@D[\[Beta][t1], t1] /. t1 -> t;

Plot[(\[Beta][t])[[1]], {t, -0.1, 0.1}, PlotRange -> All]
Plot[(d\[Beta][t])[[1]], {t, -0.1, 0.1}, PlotRange -> All]
\[Gamma][t_] := 1/Sqrt[1 - Norm[\[Beta][t]]^2];
n[r_, t_] := (r - rs[t])/Norm[r - rs[t]];
tr[r_, t_] := t - Norm[r - rs[t]]/c;
Ef[r_, t_] := (q (n[r, t] - \[Beta][t])/(\[Gamma][
      t]^2 (1 - n[r, t].\[Beta][t])^3 Norm[r - rs[t]]^2) + 
   q n[r, t] \[Cross]((n[r, t] - \[Beta][t])\[Cross]d\[Beta][t])/(
    c (1 - n[r, t].\[Beta][t])^3 Norm[r - rs[t]]) )
\[Tau] = 0.75;
StreamPlot[
 Ef[{x, y, 0}, tr[{x, y, 0}, \[Tau]]][[1 ;; 2]], {x, -1, 1}, {y, -1, 
  1}, VectorPoints -> 20]
Bf[r_, t_] := (q (
    c (\[Beta][t]\[Cross]n[r, t]))/(\[Gamma][
      t]^2 (1 - n[r, t].\[Beta][t])^3 Norm[r - rs[t]]^2) + 
   q n[r, t] \[Cross](n[r, 
        t]\[Cross]((n[r, t] - \[Beta][t])\[Cross]d\[Beta][t]))/((1 - 
       n[r, t].\[Beta][t])^3 Norm[r - rs[t]]) )
S[r_, t_] := Ef[r, t]\[Cross]Bf[r, t]

p1 = Table[
   StreamDensityPlot[{Ef[{x, y, 0}, tr[{x, y, 0}, \[Tau]]][[1 ;; 2]], 
     Log@Norm@Ef[{x, y, 0}, tr[{x, y, 0}, \[Tau]]][[1 ;; 2]]}, {x, -2,
      1}, {y, -1, 1}, ColorFunction -> "SunsetColors", 
    StreamPoints -> Medium, StreamStyle -> {White}, PlotRange -> All, 
    Frame -> False, MaxRecursion -> 5, AspectRatio -> 2/3]
   , {\[Tau], -2.5, 1.5, 0.05}];

ListAnimate[p1]

Perlesenan

Dengan ini, saya selaku pemegang hak cipta karya ini, menerbitkannya di bawah lesen yang berikut:

Fail ini disediakan di bawah Dedikasi Domain Awam Sejagat CC0 1.0 Creative Commons.

Pihak yang mengaitkan karya ini dengan surat ikatan ini telah mendedikasikan karya ini kepada domain awam dengan melepaskan segala haknya kepada karya ini di seluruh dunia di bawah undang-undang hak cipta, termasuk semua hak berkaitan dan berjiran, sejauh mana yang dibenarkan oleh undang-undang. Anda boleh menyalin, mengubah suai, mengedarkan dan mempersembahkan karya ini, termasuk untuk tujuan komersil, semuanya tanpa meminta kebenaran. http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse

Fail ini, yang pada asalnya dihantar ke https://twitter.com/j_bertolotti/status/1056921689643249666, telah disemak pada 30 October 2018 oleh penyemak Ronhjones, yang mengesahkan bahawa ia telah didapati di sana di bawah lesen yang dinyatakan pada tarikh itu.