電気力線プログラム

関数の定義

SetDirectory[ToFileName[Extract["FileName"/.NotebookInformation[EvaluationNotebook[]],{1},FrontEnd`FileName]]];
(* カレントディレクトリをノ-トブックのあるところに移動 *)

<<"FieldLines.m"

h=1.0;
r[x_,y_]:=(Sqrt[(x^2)+(y^2)]);
v[x_,y_]:=(q1/(r[x-h,y]))+(q2/(r[x+h,y]));
reg = 2;

電気力線

q1=-0.5;q2=1.0;
n = 12;

ϕ=0.;
w1=Table[FieldLine[{x,-(∂_x)v[x,y],-h+ 0.1 Cos[θ]},{y,-(∂_y)v[x,y],0.1 Sin[θ]},{t,32}],{θ,π+((2 π)/(2 n)),3 π-0.001`,((2 π)/n)}];

Ef1=Graphics[{Thickness[0.005`],RGBColor[0,0,0],w1},PlotRange->{{0,1.5 reg},{-1.5 reg,1.5 reg}}];
Show[Ef1, PlotRange->{{-1.5 reg,1.5 reg},{-1.5 reg,1.5 reg}}]