  We study FRW cosmology for scalar tensor theory where two scalar functions
nonminimally coupled to the geometry and matter Lagrangian. In a framework to
study stability and attractor solutions of the model in the phase space, we
simultaneously solve the dynamical system and best fit the stability and model
parameters with the observational data. The approach imposes restrictions on
the model constraints while providing information about the universe dynamics.
The model predict current accelerating universe, with a phantom crossing in
near future.
