% This is the function for Problem 5 of HW 1:

function f=fun4_hw1_p5(x,y);

global a g

f= -a*y + g;