Matlab

Consider the following three-dimensional nonlinear dynamical system de y1 dt dt 3 dt It is known that the solution to (4) is chaotic in time and it settles on a strange attractor. By using the function AB2.m you coded in Exercise 1, compute the numerical solution to (4). To this end, set NSTEPS-1e5, DT- 1e-3, 10STEP 50, y0-11; 2; 3, and write a function function [y,t]-solve ODEsystem() Output y: 3S matrix collecting S time snapshots of the solution to (4). Note that S-floor (NSTEPS/10STEP) +1-2001 t: 1 S vector collecting the time instants at which the solution is saved in the output matrix y. The function solve-ODE-system should also return the following items 1. The graphs of yi(t), (t) and ys(t) versus time in figure (1). 2. A three-dimensional plot of the curve (Y1(t),U2(t), (t)) in figure (2) (use the Matlab command plot3() see the documentation)