function z = f3(x,y)

global alf

 z = -4*sin(x) -4*x*cos(x) + (x*x+y*y)*sin(x) + alf*(2*x*sin(x) ...
     + (x*x+y*y) *cos(x) );