pp = [0 0 0.015 0.18 0.7 0.96 0.7 0.18 0.015 0 0];
t = -2.5:0.5:2.5;

v = -0.000:-0.001:-1.999;


for ix1 = 1:length(v),
	    disp(ix1);
 for ix2 = 1:length(v),
  p = exp(v(ix1)*t.^2+v(ix2)*t.^4);
  r(ix1,ix2) = norm(p./max(abs(p)) - pp./max(abs(pp)));
 end;
end;