dcauchy(x,l,s) = 1.0 / (pi * s * (1 + ((x-l)*1.0/s)**2))
pcauchy(x,l,s) = 0.5 + atan((x-l)*1.0/s) / pi
qcauchy(p,l,s) = l - s * 1.0/tan(pi * p)
set samples 1001
set terminal postscript enhanced color solid lw 2 "Times-Roman" 20
set output
set xtics 1
set ytics 0.1
#f(x,y,z) = dcauchy(x,y,z)
#set key 3.6,0.65
f(x,y,z) = pcauchy(x,y,z)
set key 3.6,0.2
plot [-4:4] \
f(x, 0, 1/3) title "{/Times-Italic x}_0 = 0, {/Symbol g} = 1/3", \
f(x, 0, 1) title "{/Times-Italic x}_0 = 0, {/Symbol g} = 1.0", \
f(x, 0, 1.5) title "{/Times-Italic x}_0 = 0, {/Symbol g} = 1.5", \
f(x, 0, 2) title "{/Times-Italic x}_0 = 0, {/Symbol g} = 2.0"
本作品为自由软件,您可以依据自由软件基金会发行的GNU通用公共许可证第2版或任意后续版本的条款,传播和/或修改本作品。本作品发表时预期有用,但对此无任何保证,亦无隐含的可以销售或适合特定目的的保证。详情请见GNU通用公共许可证第2版和第3版。http://www.gnu.org/licenses/gpl.htmlGPLGNU General Public Licensetruetrue