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