function plotVasicek(a,b,sigma,r0,T,delta)

N = T/delta;

for u = 1:N
    g(u) = Vasicek(u*delta,a,b,sigma,r0);
    v(u) = u*delta;
end

plot(v,g,0)

endfunction