%devct1.m: Procedure to return c derivative of the expected value
%function at a particular (w,c) point in the state space, only for t=T
%in the final horizon case.  CRRA utility

function[x]=devct2(c);
global bq wgrid i; 
dufunc1h=@dufunc1;
x=feval(dufunc1h,c)-bq*feval(dufunc1h,wgrid(i)-c);