function du_dy = economics_cd_ode(omega,u)

% ode for fourier transformed consumption deprivation equation:
% -iw u + u u' - v_0 u'' = 0
% -iw u + u z - v_0 z' = 0

% initialise v_0, the initial diffusion constant (normally time dependant)
v_0 = -2.295677427741977e+13;

% u' = z
% z' = (1/1+v_0)( -iw u + u z )

du_dy = [u(2); (1/(1+v_0))*(-1i*omega*u(1) + u(1)*u(2))];

end