#
# uspk03.rb
#
# $Id: uspk03.rb,v 1.1 2000/11/22 14:36:02 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

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

DCL::grfrm
DCL::grswnd(0.0, 1.0, rundef, rundef)

DCL::usgrph(rundef, y)

DCL::grcls

#-----------------------------------------------------------------------