set term png
set output "color_plot.png"
set size square

set xrange [-3:3]
set yrange [-3:3]


set palette model HSV functions gray,1,1 # HSV color space
#set palette model RGB rgbformulae 7,5,15 # traditional pm3d (black-blue-red-yellow)
#set palette model RGB rgbformulae 30,31,32 # color printable on gray (black-blue-violet-yellow-white)
#set palette model RGB functions gray,0,0 # gray: 0以上1以下の実数
set key outside
unset colorbox

max_n=4
plot n=0, min=15,\
"mod_data.txt" u ( min < $3 ? $1 : 1/0 ):($2):(0.1) notitle w circles lw 2 palette frac n/4.0,\
n=n+1, max=min, min=5,\
"" u (( min < $3 && $3 <= max ) ? $1 : 1/0 ):($2):(0.04) notitle w circles lw 2 lc palette frac n/4.0,\
n=n+1, max=min, min=-5,\
"" u (( min < $3 && $3 <= max ) ? $1 : 1/0 ):($2):(0.04) notitle w circles lw 2 palette frac n/4.0,\
n=n+1, max=min, min=-15,\
"" u (( $3 <= max ) ? $1 : 1/0 ):($2):(0.1) notitle w circles lw 2 palette frac n/4.0

f(x, a,b)=x/a/a*exp(-x*x/2.0/a/a)*b
a=1
b=300
fit f(x, a, b) "${OUTPUT}" u 1:2 via a,b
plot[0:5] "foo.txt" with boxes, f(x, a, b) w l lw 3 lt 3

gsl-randist 0 1000 rayleigh 1 > hoge.txt

$ which gsl-randist
/usr/local/bin/gsl-randist

[coffee@somewhere:~]gsl-randist
Usage: gsl-randist seed n DIST param1 param2 ...
Generates n samples from the distribution DIST with parameters param1,
param2, etc. Valid distributions are,

beta
binomial
bivariate-gaussian
cauchy
chisq
dir-2d
dir-3d
dir-nd
erlang
exponential
exppow
fdist
flat
gamma
gaussian-tail
gaussian
geometric
gumbel1
gumbel2
hypergeometric
laplace
landau
levy
levy-skew
logarithmic
logistic
lognormal
negative-binomial
pareto
pascal
poisson
rayleigh-tail
rayleigh
tdist
ugaussian-tail
ugaussian
weibull

gsl-randist 0 1000 rayleigh 1 > hoge.txt

cat hoge.txt | gsl-histogram 0 5 100 > foo.txt

plot [0:5]"foo.txt" u 1:3 w boxes

$ gsl-randist 0 1000 rayleigh
Error: arguments should be sigma = scale parameter

set term png truecolor

plot "3.txt"

set xtics rotate by -60
plot "3.txt"

paste data1.txt data2.txt > mod_data.txt

plot "<paste data1.txt data2.txt"

set pm3d map
set xtics 1
set ytics 1
splot "hoge.txt"

set contour base
set cntrparam level incremental 0, 20, 100
unset surface
set table "contour.dat"
splot "hoge.txt"
unset table

plot "contour.dat" w l lw 2 lc -1

set xrange [0:4]
set yrange [0:4]
set xtics 1
set ytics 1
#set isosample 250, 250
set table 'image.txt'
splot "hoge.txt"
unset table

set term png
set output "image.png"
set pm3d map
splot "image.txt"



reset
set contour base
set cntrparam level incremental 0, 20, 100
unset surface
set table 'contour.txt'
splot "hoge.txt"
unset table

set term png
set output "contour.png"
set xtics 1
set ytics 1
plot "contour.txt" w l lw 2 lc -1



reset
set term x11
set xtics 1
set ytics 1
set xrange [0:4]
set yrange [0:4]
unset key
set palette rgbformulae 33,13,10
set term png
set output "op_image.png"
p 'image.txt' with image, 'contour.txt' w l lt -1 lw 2

set term png
set output "op_image2.png"
set pm3d explicit map
splot "hoge.txt" with pm3d, "contour.txt" w l lw 2 lt -1

set terminal gif animate

ptions are 'nocrop font "arial,12" fontscale 2.0 animate delay 10 loop 0 nooptimize size 1080,720 '

plot 1

set parametric
plot [-5:5]1,t

f(x)=1.0/sqrt(2*pi)*exp(-x*x/2.0)
sigma = 1.5set lmargin 5
set rmargin 4
set tmargin 2
set bmargin 2
set size ratio 2.0/7.0
set term pdf size 7, 2
set output "gauss_twoside.pdf"
set multiplot layout 1,2
plot [-4:4][0:0.45] (-sigma < x && x < sigma ? f(x) : 1/0) with filledcurve x1 lc 3 notitle , f(x) w l lw 3 lc -1 notitle
plot [-4:4][0:0.45] (x < -sigma ? f(x) : 1/0) with filledcurve x1 lc 3 notitle , f(x) w l lw 3 lc -1 notitle, (sigma < x ? f(x) : 1/0) with filledcurve x1 lc 3 notitle
unset multiplot

set label 2 center at graph 0.5,1.1 "R.A. [degree]"

set key spacing 1.5

[x_data] [y_data] [contour_data]

set xlabel "R.A. [degree]"

set ylabel "Dec. [degree]"

set key below
#set key outside
set key spacing 0.8
#set key box

set xrange [360:0]
set yrange[-90:90]

set ztics 0,0,0

set term pdf
set output "hoge.pdf"

set contour
unset surface
set view 0,0
set cntrparam levels incremental 0,25,90
splot "data_mod.txt" w l lw 2

set pm3d corners2color c1

set key below

set key outside

set autoscale [軸]

set multiplot layout 2,1
plot "hoge1.txt" u 1:2, "hoge2.txt" u 1:2
plot "hoge1.txt" u 1:3, "hoge2.txt" u 1:3
unset multiplot

set multiplot layout 2,1
plot "hoge1.txt" u 1:2
replot "hoge2.txt" u 1:2
plot "hoge1.txt" u 1:3
replot "hoge2.txt" u 1:3
unset multiplot