gtk = graphics_toolkit();

mu_s = [0,1,2,3];
for idx = 1:length(mu_s)
  mu_noise   = 0;
  var_noise  = 1;

  mu_signal  = mu_s(idx);
  var_signal = 1;

  t   = linspace(mu_noise,mu_signal,2500);
  Pfa(:,idx) = 0.5*erfc( (t-mu_noise)/sqrt(2*var_noise) );
  Pd(:,idx)  = 0.5*erfc( (t-mu_signal)/sqrt(2*var_signal) );
 
end

plot(Pfa,Pd);

fname = ['plot_',gtk]
saveas(gcf(), [fname,'.eps']);
saveas(gcf(), [fname,'.fig']);
