f(x,y):=(mode_declare([x,y],float),
  (X^3+Y^4-0.2*X)*EXP(-X^2-Y^2)+0.3*EXP(-(X-1.225)^2-Y^2));

compile(f);

grindef(f);

plot_f():= PLOT3D('f,-3,3,-3,3);