function d = lvdemo_dist(z,pt)
  % basically the 2-norm of each element of z from pt
  
  d = (real(z - pt)).^2 + (imag(z - pt)).^2;