begin1 = 150;
end1 = 160;
begin2 = .03;
end2 = .05;
n1 = 60;
n2 = 60;
for j = 1:n2
   x2(j) = begin2 +(j-1)*(end2-begin2)/n2;
   for i = 1:n1
      x1(i) = begin1 +(i-1)*(end1-begin1)/n1;
      u(i,j) = lsfct2([x1(i) x2(j)]);
   end
end
%mesh(x1,x2,u)
contour(x1,x2,u',[8 9 9.5 10 10.5 11 12])