function [mpi] = vypocet_pi_vektorove(N)

x = rand(1,N);
y = rand(1,N);
s = sum(x.^2+y.^2 < 1);
mpi = s / N * 4;