function h = rbf(C);
  
  h = exp(-0.5*LD1(C));
