gaussint
gaussint 10000 0.0 1.0
gaussint vec 10000 0.0 1.0
vstat vec
create n 0 vec[]
lim n vec; clear; box; plot n vec
create f 0 1000
set F = rnd(f)
set r = F^2
histogram r nr 0 1 0.001
lim r nr; clear; box; plot r nr -x 1
lim r nr; clear; box; plot r nr -x 0
lim r nr; clear; box; plot r nr -x 2
lim f F; clear; box; plot f F
lim f r; clear; box; plot f r
histogram r nr 0 1 0.001 -range dr
lim dr nr; clear; box; plot dr nr -x 1
set r = sqrt(F)
histogram r nr 0 1 0.001 -range dr
lim dr nr; clear; box; plot dr nr -x 1
create f 0 10000
set F = rnd(f)
set r = sqrt(F)
histogram r nr 0 1 0.001 -range dr
histogram r nr 0 1 0.002 -range dr
lim dr nr; clear; box; plot dr nr -x 1
create f 0 1000000
set F = rnd(f)
set r = sqrt(F)
histogram r nr 0 1 0.002 -range dr
lim dr nr; clear; box; plot dr nr -x 1
set theta = rnd(f)*360.0
set x = r*dcos(theta); set set y = r*dsin(theta)
set x = r*dcos(theta); set y = r*dsin(theta)
lim x y; clear; box; plot x y
section default -imtool none; parity 0 1; densify im x y -graph; tv im -1 10
clear -graph
set r = 2000*sqrt(F)
gaussdev
gaussdev z r[] 0.0 50.0
vectors
set L = theta
set D = sqrt(r^2 + z^2)
set B = (z @ r)
set Balt = atan2 (z , r)
vectors
vstat B
vstat Balt
set dB = B - Balt
vstat dB
set RA = L
set DEC = B
csystem
csystem G C RA DEC
vectors
region 0 0 90 ait
box
tv
tv im
section default -imtool none ; parity 0 1
cdensify
cdensify img RA DEC
tv img -1 10
parity 0 1
cur
cgrid -c red
cur
cur
tvcolor rainbow
tv img -1 100
cur
stat img
set limg = log(img)
tv limg -0.1 2.0
pw
input mana.sh
mkgal; showgal
showgal
showgal
head img
echo {3600*2.7691528201E-01}
stat img
keyword img CDELT2 dX
echo $MAX
echo {$MAX / $dX^2}
mkgal; showgal
$Z_GAL = 100
mkgal; showgal
stat img
echo {$MAX / $dX^2}
input mana.sh
input mana.sh
mkgal; showgal
input mana.sh
mkgal; showgal
buffers
stat img
echo {$MAX / $dX^2}
cur
coords
coords img -c $R1 $D1
stat
stat img 670 207 20 20
echo $TOTAL
echo {$TOTAL / $NPIX / $dX^2}
echo {3000 / 5.5}
input mana.sh
mkgal; showgal
stat img 670 207 20 20
echo {$TOTAL / $NPIX / $dX^2}
stat img
echo {$MAX / $dX^2}
echo {3000 / 34}
