Recent entries

set xrange[-10:10];
set yrange[-10:10];
unset key;
set palette model RGB defined (0 'red', 0.5 'cyan', 1 'yellow');
set palette maxcolors 10;
set cbrange[0:1];
set object rect from 0,0 to 2,5 fc palette frac 1 fs noborder;
plot sin(x) lt palette frac 0.3;
last modified 22 hour(s) ago

plot sin(x) linestyle 1
last modified 2 day(s) ago

A = -5.1399428033
B = 0.6352217088
C = 2.2704694343
D = 22.6206871393

E = 2*pi*B / 360.0
F = 2*pi*C / 360.0
set samples 1200
# Either:
# set angles degrees
# plot A * sin(B*x+C) + D, 'ExtData1' lc rgb '#0060ad' pt 7 ps 1 notitle

# ... or ...
plot A * sin(E*x+F) + D, 'ExtData1' lc rgb '#0060ad' pt 7 ps 1 notitle
last modified 3 day(s) ago

plot sin(x) linestyle 1
last modified 3 day(s) ago

plot sin(x) linestyle 1
last modified 4 day(s) ago

plot ExtData1 using 1:5:6:7 with yerrorbars
last modified 4 day(s) ago

plot (sin(1/x)*x)
last modified 4 day(s) ago

plot ExtData1 using 1:2:5:6 with yerrorbars
set xlabel 'I_HHS in A'
set ylabel 'F_LS in mN'
last modified 4 day(s) ago

set xdata time
set timefmt "%Y-%m-%d %H:%M:%S"
plot '-' using 1:3 with points ls 1 # This is a test
2013-05-15 10:15:21.187 186.209729771549
2013-05-15 10:14:20.933 224.545927174927
2013-05-15 10:13:20.537 188.126047613839
2013-05-15 09:54:15.707 254.789197859055
2013-05-15 09:53:15.580 213.93299371326
2013-05-15 09:52:15.547 280.716197626502
e
last modified 4 day(s) ago

set arrow 1 from pi/2,1 to pi/2,0
set arrow 2 from pi*3/2,-1 to pi*3/2,0
plot[0:2*pi] sin(x)
last modified 5 day(s) ago

plot sin(x) linestyle 1
last modified 5 day(s) ago

#input
set timefmt "%d.%m.%Y"
set datafile separator ";"

#output
unset key
set style data lines
set style fill transparent solid 0.4
set grid

set xdata time
set xlabel 'Jahr'
set xrange ['01.01.1800':'31.12.2009']
set xtics 788400000
set mxtics 5
set format x '%Y'

set ylabel 'Einwohnerzahl'
set ytics 50000
set format y '%.0f'


set term svg size 800,400 font "Arial,10"
set output 'bochum_population.svg'

plot '-' using 1:($2) with filledcurves below x1 lt rgb 'dark-blue' lw 2
30.06.1500;1000
30.06.1650;1400
30.06.1719;1663
30.06.1750;1449
30.06.1775;1456
30.06.1800;1636
30.06.1808;2138
30.06.1812;2102
01.12.1816;2148
01.12.1819;2122
01.12.1825;2524
03.12.1834;3402
03.12.1840;3859
03.12.1843;4282
03.12.1846;4563
03.12.1849;4877
03.12.1852;5833
03.12.1855;6670
03.12.1858;8812
03.12.1861;9840
03.12.1864;11766
03.12.1867;15000
01.12.1871;21192
01.12.1875;28368
01.12.1880;33440
01.12.1885;40767
01.12.1890;47601
02.12.1895;53842
01.12.1900;65554
31.12.1901;68600
31.12.1902;69840
31.12.1903;72490
31.12.1904;115210
01.12.1905;118464
31.12.1906;122465
31.12.1907;126158
31.12.1908;129851
31.12.1909;133544
01.12.1910;136931
31.12.1911;142806
31.12.1912;148314
31.12.1913;155363
01.12.1916;137392
05.12.1917;141499
08.10.1919;142760
31.12.1919;147069
31.12.1920;150630
31.12.1921;156760
31.12.1922;160202
31.12.1923;155174
31.12.1924;156949
16.06.1925;156762
31.12.1925;158160
31.12.1926;213462
31.12.1927;216495
31.12.1928;217552
31.12.1929;321836
31.12.1930;320853
31.12.1931;318049
31.12.1932;316798
16.06.1933;314546
31.12.1933;311401
31.12.1934;313419
31.12.1935;314587
31.12.1936;313470
31.12.1937;312625
31.12.1938;312123
17.05.1939;305485
31.12.1939;318495
31.12.1940;301488
31.12.1941;280000
31.12.1942;257543
31.12.1943;185198
31.12.1944;175162
08.05.1945;160000
31.12.1945;229089
29.10.1946;246477
31.12.1946;249670
31.12.1947;265261
31.12.1948;276420
31.12.1949;285760
13.09.1950;289804
31.12.1950;292583
31.12.1951;301386
31.12.1952;312376
31.12.1953;326997
31.12.1954;335560
31.12.1955;342303
25.09.1956;345614
31.12.1956;350718
31.12.1957;357570
31.12.1958;361678
31.12.1959;364243
31.12.1960;365745
06.06.1961;361382
31.12.1961;367338
31.12.1962;367579
31.12.1963;367552
31.12.1964;365379
31.12.1965;363697
31.12.1966;361465
31.12.1967;356662
31.12.1968;354165
31.12.1969;356447
27.05.1970;343968
31.12.1970;341935
31.12.1971;341767
31.12.1972;340338
31.12.1973;338022
31.12.1974;335867
31.12.1975;414842
31.12.1976;412889
31.12.1977;409242
31.12.1978;405554
31.12.1979;402988
31.12.1980;400757
31.12.1981;398492
31.12.1982;393799
31.12.1983;389064
31.12.1984;384774
31.12.1985;382041
31.12.1986;381216
25.05.1987;386271
31.12.1987;386638
31.12.1988;389087
31.12.1989;393053
31.12.1990;396486
31.12.1991;398578
31.12.1992;400356
31.12.1993;401058
31.12.1994;401129
31.12.1995;400395
31.12.1996;398467
31.12.1997;395837
31.12.1998;393236
31.12.1999;392830
31.12.2000;391147
31.12.2001;390087
31.12.2002;388869
31.12.2003;387283
31.12.2004;388179
31.12.2005;385626
31.12.2006;383743
31.12.2007;381542
31.12.2008;378596
31.12.2009;376319
e
last modified 5 day(s) ago

set datafile separator ";"
set title 'Timing for Unsharp Masking in Spatial Domain'
set grid
set dgrid3d 20,20,20
set view 40,30,0.8,2.0
set xlabel 'KernelSize in px'
set ylabel 'ImageSize in px*px'
set zlabel 'Time in sec'
set autoscale
set surface
set pm3d
set xlabel offset -1,0,1
set ylabel offset 3,0,0
set zlabel offset -2,0,1
set samples 100
set isosamples 100
splot ExtData1 using 1:2:3 title ""
last modified 5 day(s) ago

f(x) = c*exp(-a*(x-b)**2)
fit f(x) ExtData1 via a,b,c

plot ExtData1 u 1:2 ls 1 w lp, f(x) t 'Fit' ls 2
last modified 5 day(s) ago

splot sin(x,y)
last modified 6 day(s) ago

set terminal postscript eps color font 5
set output "25vplitologias.eps"
set key left top
set title 'velocidad de Fase e=0.25 d=0.1'
set ylabel " velfase_{Vp} [m/s]"
set xlabel " theta [rad]"
set grid
#Datos de anisotropia de Thomsen y Velocidades [m/s] ondas P y S y densidad [kg/m3] clay shale
epsi1=0.334
deltaE1=0.818
delta1=0.730
gamma1=0.575
vp01=3928.00
vs01=2055.00
rho1=2590
#Datos de anisotropia de Thomsen y Velocidades [m/s] ondas P y S y densidad [kg/m3] oil shale listo
epsi2=0.200
deltaE2=0.000
delta2=-0.100
gamma2=0.145
vp02=4231.00
vs02=2539.00
rho2=2370
#Datos de anisotropia de Thomsen y Velocidades [m/s] ondas P y S y densidad [kg/m3] Green river shale listo
epsi3=0.040
deltaE3=-0.013
delta3=-0.010
gamma3=0.030
vp03=4167.00
vs03=2432.00
rho3=2310
#Datos de anisotropia de Thomsen y Velocidades [m/s] ondas P y S y densidad [kg/m3] Sandstone listo
epsi4=0.036
deltaE4=-0.037
delta4=-0.039
gamma4=0.030
vp04=4237.00
vs04=3018.00
rho4=2690
#Datos de anisotropia de Thomsen y Velocidades [m/s] ondas P y S y densidad [kg/m3] Toba
epsi5=0.020
deltaE5=-0.003
delta5=-0.030
gamma5=0.105
vp05=4846.00
vs05=1856.00
rho5=2330
#Rango para x
set xrange[0:1.6]
#c para clay shale
c11=(1-(vs01**2)/(vp01**2))
c21=4*(deltaE1/(c11**2))
c31=(4*(c11+epsi1)*epsi1)/(c11**2)
#c para oil shale
c12=(1-(vs02**2)/(vp02**2))
c22=4*(deltaE2/(c12**2))
c32=(4*(c12+epsi2)*epsi2)/(c12**2)
#c para green river shale
c13=(1-(vs03**2)/(vp03**2))
c23=4*(deltaE3/(c13**2))
c33=(4*(c13+epsi3)*epsi3)/(c13**2)
#c para green sandstone
c14=(1-(vs04**2)/(vp04**2))
c24=4*(deltaE4/(c14**2))
c34=(4*(c14+epsi4)*epsi4)/(c14**2)
#c para tobas
c15=(1-(vs05**2)/(vp05**2))
c25=4*(deltaE5/(c15**2))
c35=(4*(c15+epsi5)*epsi5)/(c15**2)




D1(x)=(1/2)*c11*((1+c21*(sin(x)**2)*(cos(x)**2)+c31*(sin(x)**4))**(1/2)-1)
vp1(x)=vp01*sqrt(1+epsi1*sin(x)**2+D1(x))


D2(x)=(1/2)*c12*((1+c22*(sin(x)**2)*(cos(x)**2)+c32*(sin(x)**4))**(1/2)-1)
vp2(x)=vp02*sqrt(1+epsi2*sin(x)**2+D2(x))


D3(x)=(1/2)*c13*((1+c23*(sin(x)**2)*(cos(x)**2)+c33*(sin(x)**4))**(1/2)-1)
vp3(x)=vp03*sqrt(1+epsi3*sin(x)**2+D3(x))


D4(x)=(1/2)*c14*((1+c24*(sin(x)**2)*(cos(x)**2)+c34*(sin(x)**4))**(1/2)-1)
vp4(x)=vp04*sqrt(1+epsi4*sin(x)**2+D4(x))


D5(x)=(1/2)*c15*((1+c25*(sin(x)**2)*(cos(x)**2)+c35*(sin(x)**4))**(1/2)-1)
vp5(x)=vp05*sqrt(1+epsi5*sin(x)**2+D5(x))


plot vp1(x) ti 'Ecuacion Thomsen (exacta)1986 Clay shale', vp2(x) ti 'Ecuacion Thomsen (exacta)1986 Oil shale', vp3(x) ti 'Ecuacion Thomsen (exacta)1986 G.R shale', vp4(x) ti 'Ecuacion Thomsen (exacta)1986 Sandstone', vp5(x) ti 'Ecuacion Thomsen (exacta)1986 Toba'
last modified 6 day(s) ago

reset
set term png enhanced truecolor size 1024,1024
unset key
set xrange [2**4:2**8]
set logscale x 2
set grid xtics mxtics ytics
set ytics 1
set multiplot layout 5,2
set title 'default'; set mxtics 8; plot log(x)/log(2)
set title 'xtics, 8'; set xtics 4; set mxtics 8; plot log(x)/log(2)
set title 'xtics, 2'; set xtics 4; set mxtics 2; plot log(x)/log(2)
set title 'xtics, 3'; set xtics 4; set mxtics 3; plot log(x)/log(2)
set title 'xtics, 4'; set xtics 4; set mxtics 4; plot log(x)/log(2)
set title 'xtics, 5'; set xtics 4; set mxtics 5; plot log(x)/log(2)
set title 'xtics, 6'; set xtics 4; set mxtics 6; plot log(x)/log(2)
set title 'xtics, 7'; set xtics 4; set mxtics 7; plot log(x)/log(2)
set title 'xtics, 8'; set xtics 4; set mxtics 8; plot log(x)/log(2)
set title 'xtics, 9'; set xtics 4; set mxtics 9; plot log(x)/log(2)
unset multiplot
last modified 6 day(s) ago

plot 'ExtData1' with linespoints
set xrange [0:0.0005]
plot 'ExtData2' with linespoints
last modified 6 day(s) ago

s(e)=sqrt((1-e)/(1+e))
E(x,e)=2*atan(s(e)*tan(x/2))
e=0.99
M(x,e)=E(x,e)-e*sin(E(x,e))
set grid
w=0

plot [0:1] x ls 1
#(M(pi/2-w,x)-M(0-w,x))/pi ,w=pi/2,(M(pi/2-w,x)-M(0-w,x))/pi, w=pi/4,
(M(pi/2-w,x)-M(0-w,x))/pi, w=pi/3,(M(pi/2-w,x)-M(0-w,x))/pi, w=pi*0.1,(M(pi/2-w,x)-M(0-w,x))/pi

last modified 6 day(s) ago

unset dgrid3d
set hidden3d
set pm3d map
plot ExtData1 using 2:1:3 with image
plot ExtData1 using 2:1:3 with points pt 5 ps 5 palette
last modified 6 day(s) ago

1 - 2 - 3 - 4 - 5 - 6 - 7 - 8 - 9 - 10 - 11 - 12 - 13 - 14 - 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24 - 25 - 26 - 27 - 28 - 29 - 30 - 31 - 32 - 33 - 34 - 35 - 36 - 37 - 38 - 39 - 40 - 41 - 42 - 43 - 44 - 45 - 46 - 47