plot_sol_script.m