How to solve two generation uility funcion on DSGE model, or if not how to solve it using other model?

Have a look at