
load theta1.mat
solutions_Final_p2c(t, theta1, 0.6, 0.8, 0.1);