function f = f( x, y, a, b,c )
     f=exp(a*(2*x-a)+b*(2*y-b)+2*c*(x-a)*(y-b));
end

