solve_model function

Solve_model