code prepared for course

sec3_predefined_utilities.m

@@ -58,18 +58,17 @@ F        = integrator('F', 'cvodes', ode, opts);
 figure; grid on; hold on;
 
 xi = [2; 0];
-trajectory_vals = [xi];
+traj_vals = [xi];
 for i = 1:floor(T / T_sample)
     if (T_sample * i < T_star)
-        trajectory = F('x0', trajectory_vals(:, end), 'p', 1);
+        trajectory = F('x0', traj_vals(:, end), 'p', -1);
     else
-        trajectory = F('x0', trajectory_vals(:, end), 'p', -1);
+        trajectory = F('x0', traj_vals(:, end), 'p', 1);
     end
-    trajectory_vals = [trajectory_vals, full(trajectory.xf)];
+    traj_vals = [traj_vals, full(trajectory.xf)];
 end
 %get the trajectory samples
-plot(trajectory_vals(1, :), trajectory_vals(2, :));
-
+plot(traj_vals(1, :), traj_vals(2, :));
 
 %% If else instruction
 % Saturate function