(format t "Woot~%")
  (new-random-seheult 1.4 1.6 1 1)
  (let (tempval)
    (handler-case
        (progn
          (setf tempval (sh-model x y sd mu tau)))
      
      (floating-point-overflow ()
        (setf tempval 10)))
    (if (> tempval 10)
        (setf tempval 10))
    (with-open-file (out "randomdata.dat"
                         :direction :output
                         :if-exists :append)
      (format out "~%~%~a ~a ~a" x y tempval))
    tempval)
  (format-gnuplot "set title 'Seheult Data'")
  (format-gnuplot "unset key")
  (set-labels "x-ri" "y-ri" "LR")
  (new-generate-thumbs 30 40 40 30)
  (format-gnuplot "splot 'randomdata.dat' index 0, 
'randomdata.dat' index 1"))