%dufunc1.m: Procedure to evaluate the derivative of the utility function for lifecycle saving problem. function[x]=dufunc1(c); global gama; x=(1/c^gama);