Add parallel construct to inner loop of solve_and_plot()