add_highs_solver function

Add a HiGHS solver