function f = f(k)        

global alpha delta

f = k.^alpha + (1 - delta)*k;

