#! /bin/sh

wait=$1
file=$2
#----------------------------------------------------------
sm -s <<FIN
if ($wait) {xterm_s} else {laser_sp $file\n lweight 2}
#----------------------------------------------------------

#-----------------------read data
data file.dat
read {x 1 y 2}

#-----------------------create a new array
set q = 0,10,0.1

#-----------------------do some array math and print to a file
set z = x**2 + sqrt(y)
print output.dat <z>

#-----------------------do the same thing using a loop
set dimen(w)=dimen(x)
do i=0,dimen(x)-1,1 {
  
  set w[\$i] = x[\$i]**2 + sqrt(y[\$i])
  echo \$(w[\$i])
}

#-----------------------make plot
erase limits x y 
expand 1.1
box
expand 1

ptype 8 3
points x y

ltype 0
connect x y

expand 1.5
xlabel x
ylabel y
expand 1

#----------------------------------------------------------
if ($wait) {!sleep $wait} else {hardcopy}
FIN
#----------------------------------------------------------