making solve1.m robust to sparse Jacobian
adding solve_algo=9: using trust_region algorithm on the entire model