解説Rational Elliptic Function (abs, n=3, x=(0,5)).svg	English: A graph of thr absolute value of the third-order rational elliptic function, R3(ξ,x) over the interval [0,5]. Here, ξ (the selectivity factor) is 1.4. Also shown is the discrimination factor, Ln.
日付	2009年1月29日
原典	投稿者自身による著作物
作者	Inductiveload
許可 (ファイルの再利用)	Public domainPublic domainfalsefalse この作品の著作権者である私は、この作品についての権利を放棄しパブリックドメインとします。これは全世界で適用されます。 一部の国では、これが法的に可能ではない場合があります。その場合は、次のように宣言します。 私は、あらゆる人に対して、法により必要とされている条件を除き、如何なる条件も課すことなく、あらゆる目的のためにこの著作物を使用する権利を与えます。

xp2[xi_] :=
  Module[{g, num, den},
   g = Sqrt[4*xi^2 + (4*xi^2*(xi^2 - 1))^(2/3)];
   num = 2*xi^2*Sqrt[g];
   den = Sqrt[8*xi^2*(xi^2 + 1) + 12*g*xi^2 - g^3] - Sqrt[g^3];
   num/den
   ];
xz2[xi_] := xi^2/xp2[xi];

t[xi_] := Sqrt[1 - 1/xi^2];

(*Use the particular forms for these low-order REFs*)
r1[xi_, x_] := x;
r2[xi_, x_] := ((t[xi] + 1)*x^2 - 1)/((t[xi] - 1)*x^2 + 1);
r3[xi_, x_] := 
  x*((1 - xp2[xi])*(x^2 - xz2[xi]))/((1 - xz2[xi])*(x^2 - 
        xp2[xi]));
r4[xi_, x_] :=
  Module[{num, den},
   num = (1 + t[xi]) (1 + Sqrt[t[xi]])^2*x^4 - 
     2 (1 + t[xi]) (1 + Sqrt[t[xi]])*x^2 + 1;
   den = (1 + t[xi]) (1 - Sqrt[t[xi]])^2*x^4 - 
     2 (1 + t[xi]) (1 - Sqrt[t[xi]])*x^2 + 1;
   num/den
   ];
   
LogPlot[
 xi = 1.4;
 Abs[r3[xi, x]],
 {x, 0, 5},
 PlotRange -> {0.01, 1000}]