Please do this in Matlab Thanks

Consider the Van der Pol oscillator

$y^{\text{'}\text{'}}-(1-y^2)y^{\text{'}}+y=F(t),0t100$

$y(0)=0,y^{\text{'}}(0)=2$

You can first let $y_1=y$ and split the IVP $(1)$ into two ODE's

$y_1^{\text{'}}=y_2,$

$y_2^{\text{'}}=(1-y_1^2)y_2-y_1+F(t)$

$y_1(0)=0,y_2(0)=2$

Part A: $F(t)=0.$

Part B: $F(t)=$Asin$(t).$

Please do the following:

I. Solve the IVP $(2)$ with $=\{0,0\mathrm{.}1,0\mathrm{.}5,1,5,10\},A=1\mathrm{.}5,=\frac{}{5},$ a your choice of

step size $h$ and generate solution plots $(y$ vs $t$ plot$)$ using the Euler's Method

II$.($Only for Part A$)$ Make $y_1$ vs $y_2(y_1$ is the $x-$axis; $y_2$ is the $y-$axis$)$ plot for each case

of your $.$

III. On top of the plot in part II$,$ draw the direction field plot.