LAB$06$ex$1.$m

clear all; $\%$this deletes all variables

omega $=3$;$c=2$; omega $=3\mathrm{.}6$;

param $=[$omega$0,$ c$,$omega$]$;

$t0=0$;$y0=0$;$v0=0$;$Y0=[y0$;$v0]$;$tf=40$;

options $=$ odeset$(\text{'}$AbsTol$\text{'},1e-10,$ 'relTol', $1e-10$;

$,$ options, param $)$;

$y=Y($:$,1)$;$,v=Y($:$,2)$;

figure $(1)$

plot$($t$,$y$,\text{'}$b$-\text{'})$; ylabel$(\text{'}$y$\text{'})$; grid on;

$t1=13$;$i=$ find $(t>t1)$;

$C=\frac{\text{m}\text{a}\text{x}(Y(i,1))-\text{m}\text{i}\text{n}(Y(i,1))}{2}$;

disp$([\text{'}$computed amplitude of forced oscillation $=\text{'},$ num$2$str$($C$)])$;

Ctheory omega$02$