exec("scilab/question-1-1.sci");

function [x] = inverse(u)
  x = cdfnor("X",zeros(u),ones(u),u,ones(u)-u);
endfunction

N=10000;
X=[vector_vdc(2,N);vector_vdc(3,N)];
G_VDC=inverse(X);

Covar(G_VDC(1,:),G_VDC(2,:));
histplot(100,G_VDC(1,:));
xbasc();
isoview(0,1,0,1);
plot2d(G_VDC(1,:),G_VDC(2,:),0);


X=[rand(1,N,"unif");rand(1,N,"unif")];
G_RAND=inverse(X);
Covar(G_RAND(1,:),G_RAND(2,:));
histplot(100,G_RAND(1,:));
xbasc();
isoview(0,1,0,1);
plot2d(G_RAND(1,:),G_RAND(2,:),0);