```from rsf.proj import * ## If you change the order of lines here, do not forget to modify ../rsfplot.tex # Graph, x linear example Flow('sinc',None,'math n1=101 d1=0.2 o1=-10 output="sin(x1)/x1"') Result('sinc','graph title=Sinc Function') # Graph, x squared example Flow('xsq',None,'math n1=10 output="x1*x1"') Flow('sqline','xsq','cmplx \$SOURCE') Result('sqline',''' graph symbol="+" title="Line, x\\^2" symbolsz=12''') Flow('msinc',None, '''math n1=101 d1=0.2 o1=-10 n2=3 output="sin((x2+1)*x1)/((x2+1)*x1)"''') Result('msincs','msinc','graph title="Multiple Sincs, n3 > 1"') Result('mpsincs','msinc',''' graph title="With symbol=" symbol=+*o symbolsz=8''') Flow('mline',None,'math n1=10 n2=3 output="(x2+1)*x1"') Result('dplines','mline',''' dots title="Lines with bars" gaineach=n radius=0.15''') # Box example Plot('label1',None,''' box y0=5.6 x0=6.4 xt=-1 yt=1 label="Line \F10 a"''') Plot('msinc','msinc','graph title="Sincs with a label (overlay)"') Result('mbmincs','msinc label1','Overlay') # Surface and contours Flow('sinsq',None,''' math n1=201 n2=201 o1=-1.0 o2=-1.0 d1=0.01 d2=0.01 output="sin(3.14*x1)*sin(3.14*x2)" label1=x label2=y''') Result('sqsurf','sinsq','grey color=j title=Surface scalebar=y') Result('csqsurf','sinsq','contour title=Isolines allpos=n') # Surface and contours, overlay Plot('psinsq','sinsq','''grey color=j title=Overlay scalebar=y min1=-1 max1=1 min2=-1 max2=1''') Plot('pcsinsq','sinsq','contour title=Overlay allpos=n scalebar=y') Result('csurf','psinsq pcsinsq','Overlay') # Wiggle example Flow('signal',None,'''spike nsp1 n1=301 d1=0.004 o1=0.0 k1=101 l1=101 title= label1=Time unit1=s | ricker1 frequency=15 | spray axis=2 d=0.01 n=51 o=0 label=Distance unit=km''') Flow('vel','signal','math output=2') Flow('gath',['signal','vel'],'inmo velocity=\${SOURCES[1]}') Plot('wgath','gath','''wiggle poly=y min2=0 max2=0.5 plotcol=7 yreverse=y wherexlabel=top gridcol=7''') Plot('rgath','gath','grey transp=n min2=0 max2=0.5') Result('sgath','rgath wgath','Overlay') # Cubeplot example Flow('ganom',None,'''math n1=41 n2=41 n3=41 o1=-1 o2=-1 o3=-1 d1=0.05 d2=0.05 d3=0.05 output="exp(-2*(x1^2+x2^2+x3^2))"''') Result('fganom','ganom','''byte gainpanel=all bar=gbar.rsf allpos=y | cubeplot color=j scalebar=y bar=gbar.rsf title="Slices, flat" frame1=20 frame2=20 frame3=20''') # Cubeplot, non-flat example Result('ganom3','ganom','''byte gainpanel=all bar=gbar.rsf allpos=y | cubeplot color=j scalebar=y bar=gbar.rsf title="Cube" frame1=20 frame2=20 frame3=20 flat=n''') # SideBySide example Flow('ganom1',None,'''math n1=41 n2=41 n3=41 o1=-1 o2=-1 o3=-1 d1=0.05 d2=0.05 d3=0.05 output="exp(-2*(x1^2+x2^2+x3^2))"''') Plot('ganom','ganom','''byte gainpanel=all bar=gbar.rsf allpos=y | cubeplot color=j scalebar=y bar=gbar.rsf title="Surface" frame1=20 frame2=20 frame3=20''') Plot('ganom1','ganom1','''contour3 color=j scalebar=y title=Isolines nc=30 frame1=20 frame2=20 frame3=20''') Result('sganom','ganom ganom1','SideBySideAniso') Result('ganom2','ganom ganom1','SideBySideIso') Result('ganom-7','ganom1 ganom','Overlay') Flow('data1',None,'spike n1=30 k1=15 label1=" " ') Flow('data2','data1','causint') Result('dots','data1 data2',''' cat axis=2 \${SOURCES[1:2]} | dots title="Dirac Delta and Step functions" gaineach=n radius=0.15 yreverse=1 overlap=.8 constsep=1 connect=0''') Result('dplines2','data1 data2', ''' cat axis=2 \${SOURCES[1:2]} | graph ''') End()```