Fo = 1.0; Omega = 1.0; f[t_] := Fo Sin[Omega t]Plot[f[t], {t, 0, 2 π}, AxesLabel -> {"t", "f[t]"}]

k = 1. ; m = 1. ; b = 0.3; Fo = 1. ; Omega = 1. ; ti = 0. ; tf = 35. ; delta = 0.2; x[ti] = 0. ; v[ti] = 0. ; Do[F[t] = -k x[t] - b v[t] + Fo Sin[Omega t]; 
      v[t + delta] = v[t] + ((F[t] delta)/m); 
      x[t + delta] = x[t] + v[t + delta] delta, {t, ti, tf, delta}]

xdata = Table[{t, x[t]}, {t, ti, tf, delta}]; ListPlot[xdata, AxesLabel -> {"t", "x[t]"}, 
      PlotLabel -> "Driving Freq = 1.0 "]

(*    Make a list of the resonance data obtained from running 
the previous program.  Make a graph of the "resonance curve".                            *)
       
ResonanceData = { {0.2, 1.0}, {0.5, 1.4}, {0.75, 2.0}, {1.0, 3.3},
                  {1.25,1.6}, {1.5, 0.8}, {1.9, 0.4} };
ListPlot[ResonanceData,AxesLabel->{"Omega","Amplitude"}],
      PlotStyle->PointSize[0.02],PlotRange->{{0,2},{0,4}}]