Added calls to integration array Y and derivative DYDX in main loop
Added k-grid