function [r]=randn()
    r=rand(1,"normal");
endfunction

function [w]=wiener(dt,n)
    w(1)=0; 
    for i=1:n
        dw=sqrt(dt)*randn();
        w(i+1)=w(i)+dw;
    end;
endfunction

// ---------------------------------
// I. Wienerov proces

dt=0.001;
n=1000;
cas=(0:dt:n*dt);

figure(1);
a=gca();
a.background=color("white");
a=gcf();
a.background=color("white");

plot(cas,wiener(dt,n));

// -------------------------
// II. Brownov pohyb
// DOPISTE ZVYSOK SKRIPTU  - chceme niekolko realizacii a strednu hodnotu

//cas=(0:dt:n*dt);
//mi=3;sigma=0.5;
//
//figure(2);
//a=gca();
//a.background=color("white");
//a=gcf();
//a.background=color("white");
//
//plot(cas,mi*cas+sigma*wiener(dt,n));


// -------------------------
// III. Geometricky Brownov pohyb
// DOPISTE ZVYSOK SKRIPTU  - chceme niekolko realizacii a strednu hodnotu

// plot(cas,x0*exp(mi*cas+sigma*wiener(dt,n)),'r');


