close all
clear all;
pkg load parallel

function [ mPV ]= get_ests(th0)
    % rand("state","reset") % this  helps ...
    for itr=1:100 % run over trials and estimate each time
        t(itr)=randn();
    end
    mPV= mean(t);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

th_ar = 1:14;

mPV_ar= pararrayfun(nproc,@(th) get_ests(th), th_ar)

