''
'' cosc611HW7.sim
'' simc cosc611HW7.sim
'' 
'' Forgot that pulling a random number has to be the same location as 
'' computing a variance value since there can not be two random value
'' from pulling from a list when computing a mu and S^2.
'' it has to be the same number from the list for both computations
'' old script screwed up itself. 
'' sim script.
''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
PREAMBLE
define i as integer variable
define z as integer variable
define tmp as integer variable
define total as integer variable
define mu as real variable
define varvalue as real variable
define muvalue as real 1-dim array
define arrlist as integer 1-dim array
define muarr as real 1-dim array
define vararr as real 1-dim array
define minarr as real 1-dim array
define maxarr as real 1-dim array
define mufinal as real 1-dim array
END

MAIN
'' file output
open 3 for output, name="data.txt"
use 3 for output

reserve muvalue(*) as 6
reserve arrlist(*) as 100
reserve muarr(*) as 100
reserve vararr(*) as 100
reserve minarr(*) as 100
reserve maxarr(*) as 100
reserve mufinal(*) as 100

'' array values
arrlist(1)=21
arrlist(2)=23
arrlist(3)=25
arrlist(4)=26
arrlist(5)=28
arrlist(6)=29
arrlist(7)=31
arrlist(8)=33
arrlist(9)=35
arrlist(10)=41
arrlist(11)=21
arrlist(12)=23
arrlist(13)=25
arrlist(14)=27
arrlist(15)=28
arrlist(16)=29
arrlist(17)=31
arrlist(18)=33
arrlist(19)=36
arrlist(20)=41
arrlist(21)=22
arrlist(22)=24
arrlist(23)=25
arrlist(24)=27
arrlist(25)=28
arrlist(26)=29
arrlist(27)=31
arrlist(28)=33
arrlist(29)=36
arrlist(30)=42
arrlist(31)=22
arrlist(32)=24
arrlist(33)=25
arrlist(34)=27
arrlist(35)=28
arrlist(36)=29
arrlist(37)=31
arrlist(38)=34
arrlist(39)=36
arrlist(40)=42
arrlist(41)=22
arrlist(42)=24
arrlist(43)=26
arrlist(44)=27
arrlist(45)=28
arrlist(46)=29
arrlist(47)=31
arrlist(48)=34
arrlist(49)=37
arrlist(50)=45
arrlist(51)=22
arrlist(52)=24
arrlist(53)=26
arrlist(54)=27
arrlist(55)=28
arrlist(56)=29
arrlist(57)=31
arrlist(58)=34
arrlist(59)=38
arrlist(60)=46
arrlist(61)=23
arrlist(62)=24
arrlist(63)=26
arrlist(64)=27
arrlist(65)=28
arrlist(66)=30
arrlist(67)=32
arrlist(68)=34
arrlist(69)=39
arrlist(70)=47
arrlist(71)=23
arrlist(72)=25
arrlist(73)=26
arrlist(74)=27
arrlist(75)=28
arrlist(76)=30
arrlist(77)=33
arrlist(78)=34
arrlist(79)=39
arrlist(80)=52
arrlist(81)=23
arrlist(82)=25
arrlist(83)=26
arrlist(84)=27
arrlist(85)=28
arrlist(86)=30
arrlist(87)=33
arrlist(88)=35
arrlist(89)=40
arrlist(90)=54
arrlist(91)=23
arrlist(92)=25
arrlist(93)=26
arrlist(94)=28
arrlist(95)=29
arrlist(96)=30
arrlist(97)=33
arrlist(98)=35
arrlist(99)=40
arrlist(100)=58

'' Calculate MU value and calculate S^2 = Variance
''muvalue = 0
for z = 1 to 100
do
    for i = 1 to 6
    do
      let muvalue(i) = arrlist(uniform.f(1, 100, dim.f(seed.v)))
      compute mu as average of muvalue(i)
      compute varvalue as variance of muvalue(i)
    loop
    let muarr(z) = mu
    let vararr(z) = varvalue
loop

'' Calculate Min interval
for i = 1 to 100 
do
   minarr(i) = muarr(i) - 2.015*sqrt.f(vararr(i)/5)
loop

'' Calculate Max interval
for i = 1 to 100 
do
   maxarr(i) = muarr(i) + 2.015*sqrt.f(vararr(i)/5)
loop

'' determine 90% data is less than mu
'' 1=true; 0=false
for i = 1 to 100
do
   if 30.64 > minarr(i) and 30.64 < maxarr(i)
       mufinal(i) = 1
        else
       mufinal(i) = 0        
    endif
loop
          
print 1 line thus
REP   RandSample   X(6)       S(6)      Cal-Min/Max-Interval   MUWithinInterval
for i = 1 to 100
print 1 line with i,  6  , muarr(i),  vararr(i), minarr(i),maxarr(i), mufinal(i) thus
***    *          **.****   **.****      **.****    **.****    *

'' end main
end