!set n=$teller
bewerking=bewerking1.proc
wims_read_parm=$empty
nivo_title=Bereken de Mediaan,gemiddelde en nog wat meer...    
counter=0    
data$n=<table border="1">
gemiddelde$n=$empty
som=0
!if $graad=0
    aantal_vragen=5
!else
    aantal_vragen=2
!endif            
!if $soort=1
    klaar=1
    counter=$maxcounter
    !for p=1 to $maxcounter
	meting=!randint 1,$afrondingsfactor
	meting$p=$[1000*(1+ $meting/$afrondingsfactor)]
	som=$[$som + $(meting$p)]
	wims_read_parm=!append item $(meting$p) to $wims_read_parm
	data$n=!append line <td>meting&nbsp;$p</td><td>$(meting$p)&nbsp;seconden</td><tr> to $(data$n)
    !next p
!else
    !if $klaar=1
	!for p=1 to $maxcounter
	    som=$[$som + $(meting$p)]
	    wims_read_parm=!append item $(meting$p) to $wims_read_parm
	    data$n=!append line <td>meting&nbsp;$p</td><td>$(meting$p)&nbsp;seconden</td><tr> to $(data$n)
	!next p
    !else
	 !exit
    !endif
!endif

data$n=$(data$n)</table>
# gemiddelde
GOED1$n=$[$som/$maxcounter]

# standaard wafwijking
# is s geen sigma
# !readproc slib/stat/deviation $wims_read_parm 
# GOED2$n=!item 2 of $slib_out
    
NW=0
!for p=1 to $maxcounter
    w=!item $p of $wims_read_parm
    W=$[($w-$(GOED1$n))^2]
    NW=$[$W + $NW]
!next p
    
GOED2$n=$[sqrt($NW/$maxcounter)]


# eerste & derde kwartiel
    
wims_read_parm=!sort numeric list of $wims_read_parm
!if $[gcd($maxcounter,2)]=2
    even=1
    # mediaan
    # !readproc slib/stat/median $wims_read_parm 
    # GOED3$n=$slib_out
    m1=!item $[0.5*$maxcounter] of $wims_read_parm
    m2=!item $[0.5*$maxcounter+1] of $wims_read_parm
    GOED3$n=$[($m2+$m1)/2]
    
    !if $[gcd((0.5*$maxcounter),2)]=2
        w1=$[0.25*$maxcounter]
	w11=$[0.25*$maxcounter+1]
	w3=$[0.75*$maxcounter]
	w33=$[0.75*$maxcounter+1]
	q1=!item $w1 of $wims_read_parm
	q11=!item $w11 of $wims_read_parm
	q3=!item $w3 of $wims_read_parm
	q33=!item $w33 of $wims_read_parm
    !else
	w1=$[0.25*$maxcounter+0.5]
	w3=$[0.75*$maxcounter+0.5]
	w33=$w3
	q1=!item $w1 of $wims_read_parm
	q11=$q1
    	q3=!item $w3 of $wims_read_parm
	q33=$q3    
    !endif
!else
    even=0
    # mediaan
    # !readproc slib/stat/median $wims_read_parm 
    # GOED3$n=$slib_out
    GOED3$n=!item $[($maxcounter+1)/2] of $wims_read_parm
    $[gcd((0.5*($maxcounter-1)),2)]
    !if $[gcd((0.5*($maxcounter-1)),2)]=2
	w1=$[0.25*$maxcounter - 0.25]
	w11=$[0.25*$maxcounter + 0.75]
	w3=$[0.75*$maxcounter + 0.25]
	w33=$[0.75*$maxcounter + 1.25]
	q1=!item $w1 of $wims_read_parm
	q11=!item $w11 of $wims_read_parm
	q3=!item $w3 of $wims_read_parm
	q33=!item $w33 of $wims_read_parm
    !else
	w1=$[0.25*$maxcounter+0.25]
	w3=$[0.75*$maxcounter+0.75]
	w33=$w3
	q1=!item $w1 of $wims_read_parm
	q11=$q1
	q3=!item $w3 of $wims_read_parm	
	q33=$q3
    !endif
!endif

    
GOED4$n=$[($q1+$q11)/2]
GOED5$n=$[($q3+$q33)/2]
goed1$n=Het gemiddelde is dus $(GOED1$n)
goed2$n=De standaarddeviatie $m_sigma is $(GOED2$n)
goed3$n=De mediaan $m_mu is $(GOED3$n)
goed4$n=Het eerste kwartiel is $(GOED4$n)
goed5$n=Het derde kwartiel is $(GOED5$n)
