clear;
format long;

n=2;  % 2 paramters
m=20;  % 20 ensemble members

r=zeros(m,n);
x=zeros(m,n);
y=zeros(1,n);
mu=[0.95 10];  % reference values from Jon
col=[0.3; 5.0]; % variances
row=[0.3 5.0];   % variances
sigma=col*row;
%sigma=[1 1.5; 1.5 3];
%r=mvnrnd(mu,sigma,m);
%plot(r(:,1),r(:,2),'+g');
%hold on;
x=lhsnorm(mu,sigma,m); % call routine
plot(x(:,1),x(:,2),'or');
hold on;
xlabel('qvsurf');
ylabel('evaplimit')
save temp.dat x -ascii  % create ascii file of the results
