distances = {'kl', 'var', 'l1'};
for n=0:9,
  p_joint = rand(2);
  for distid=1:3,
    dist = distances{distid}
    mykl(p_joint, dist);
    saveas(gcf, sprintf('%s-%d', dist, n), 'png');
    close(gcf)
  end
end
