!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
define real ndat
define real tsdev
let tsdev 0.02
!
for idir 1 2 
    pause "Type C to continue"
    ndat = zfocus%dat'idir'%n
    let zfocus%dat'idir'%x[i] = 4*(i-0.5*(ndat+1))/ndat
    exa zfocus%dat'idir'%x  ! [i] = 2*(i-0.5*ndat)
    let zfocus%dat'idir'%y = 1/(1+((zfocus%dat'idir'%x-1.5+idir)/3)^2)                
    for idat 1 to ndat
        say 'noise(tsdev)'
        let zfocus%dat'idir'%y[idat] =  zfocus%dat'idir'%y[idat]+noise(tsdev)
    next
    let zfocus%fun'idir'%flag = 0.
    focus L
next
PAUSE "Continue to see Parabola fitting"
focus P
return
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
