#
# uupk01.rb
#
# $Id: uupk01.rb,v 1.1 2000/11/23 08:19:14 keiko Exp $
#

require "narray"
require "numru/dcl"

include NumRu
include Math


n = 30
y = NArray.sfloat(n)

#---------------------------- DATA DEFINITION --------------------------
r = 0.2
a = 3.7
for i in 0..n-1
  r = a*r*(1.0-r)
  y[i] = r
end

y[n/2-1] = 999.0
DCL::gllset('LMISS', true)
rundef = DCL::glrget('RUNDEF')

#----------------------------- GRAPH -----------------------------------
iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
DCL::gropn iws

DCL::grfrm
DCL::sglset('LCLIP', true)

DCL::grswnd(0.0, 1.0, rundef, rundef)
DCL::usspnt(rundef, y)
DCL::uspfit
DCL::grstrf

DCL::usdaxs

DCL::uuslni(2)
DCL::uulin(rundef, y)

DCL::uusmkt(4)
DCL::uumrk(rundef, y)

DCL::grcls