I am having trouble getting this to work in mat lab and plot

this is what i was initally trying but I'm not very adept with matlab

if you could comment out what you code is doing to help me understand mat lab commands I would appreciate it

i could get the forcing function as a piece wise to plot but that's about all that worked for me and I think the values of an(i) was working...

t=[-pi:.1:pi] pw = piecewise( (-pi

clc clear syms t T=2*pi; w=(2*pi)/T; n = 10; step = .1;

for i = 1 : n-1 an(i) =(2/T)*( ( int( (t+(pi/2))*cos(i*w*t),t,(-pi/2),0 ))*( int( (-t+(pi/2))*cos(i*w*t),t,0,(pi/2)) )); n(i)=i; end %t3=[-pi:step:pi] %for ii = 1 : n-1 % C(ii) = cos(w*n*t) %n(ii)=ii %end

[15 points Consider a damped mass-spring system undergoing forced oscillations governed by the ODE 3. y +0.05y 25y (t), where the driving force is a non-sinusoidal periodic function given by -