function VOI_Gaussianprojects;

% vector of correlations
rhoV=[-0.99:0.09:0.99];
tauvec=[0.1 0.5 1];

   for i=1:length(rhoV),
    rho=rhoV(i);
for j=1:length(tauvec),
    tau=tauvec(j);
 
    
    % IMPERFECT    
    sumIMP2(i,j)=(2/sqrt(2*pi))*sqrt(((1+tau^2)*(1+rho^2)-2*rho^2)/((1+tau^2)^2-rho^2));

     
    end
       % PERFECT
    sumPERF2(i)=(1+abs(rho))/sqrt(2*pi);

end

figure(1)
clf;
plot(rhoV,sumPERF2,'kx');
hold;
plot(rhoV,sumIMP2(:,1),'-.k')
plot(rhoV,sumIMP2(:,2),'o-k')
plot(rhoV,sumIMP2(:,3),'-k')

sl=legend('Perfect partial','Total Imperfect, \tau=0.1','Total Imperfect, \tau=0.5','Total Imperfect, \tau=1','Location','SouthEast');
set(sl,'FontSize',14);
sy=ylabel('VOI');
sx=xlabel('Correlation parameter');
set(sx,'FontSize',14);
set(sy,'FontSize',14);
axis([-1 1 0 0.9])

