#! /bin/bash
# Last edited on 2009-11-30 00:12:42 by stolfi

# Plots the windowing function {omega(B,E,t)}


epsfile="/tmp/$$.eps"

gnuplot <<EOF
set term postscript eps color solid "TimesRoman" 24
set size 1.0,0.75
set output "${epsfile}"
set nokey

B=20
E=70
W=10

set xzeroaxis
set xtics 10
set grid xtics
set xrange[(B-2*W):(E+2*W)]

set yrange [-0.2:+1.2]
set ytics 0.50 format "%4.2f"
set mytics 5

sigma(z)=(abs(z) > 1 ? 0/0 : sin(0.5*pi*sin(0.5*pi*z)))
omega(z)=(z<=B-W ? 0 : (z < B+W ? 0.5*(1+sigma((z-B)/W)) : (z <= E-W ? 1 : (z < E + W ? 0.5*(1+sigma((E-z)/W)) : 0 ))))

plot omega(x) notitle with lines lw 3.0 lt rgb '#227700'
quit
EOF

cat ${epsfile}
gv ${epsfile}
rm -f ${epsfile}