function pf = f(k,l)        

global alpha delta A

pf = A.*k.^alpha.*l.^(1-alpha) + (1 - delta).*k;