// This function returns the value of a Euclidian Distance minimum
// distance classifier at x, given the mean vectors.
 
function [r]= ed_disc (x, U1, U2)

   g1 = sqrt ((x-U1)' * (x-U1));
   g2 = sqrt ((x-U2)' * (x-U2));
 
   r = g1 - g2;