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