function dx = du(C, sigma_param)
    dx = C^(-sigma_param);  % First derivative of CRRA utility
end
