[NEW] Pass solver name as argument for functions solving optimization models

Make all optimization models accept a solver_name argument, set default to glpk. Make it possible to pass in the solver_name when calling the function from command line.