# Kommenttirivi: Tama on octave ohjelmani Osub2.m
function [n, z] = rayleightest (t, f)
   n=length(t) ;
   x=2.*pi*f*t ;
z=(1./n)*(sum(cos(x))^2 + sum(sin(x))^2);
endfunction
# ____________________________________________________________
clear ; clc                        # Poistetaan ... Tyhjennetaan ...
t=rand(10,1)*10.D0 ;               # Satunnaisajat 0 < t_i < 10 (n=10)
f=1. ;                             # Frekvenssi on f=1.
[n, z]=rayleightest(t,f) ;         # Rayleigh testin n ja z(f)
printf("n = %f\n",n)               # Tarkistus
printf("z = %f\n",z)               # Tarkitus
m=rayleightest(t,f);               # Koe: Kaytetaan vain yhta muuttujaa 
printf("m = %f\n",m)               # Tulos: antaa VAIN ensimmaisen arvon
printf("class(m) = %s\n",class(m)) # String
