d[l_, g_ , list_] := 
 Table[{list[[k]][[1]], l*Sin[list[[k]][[2]]]/Tau[l, g]^2}, {k, 1, 
   Length[list]}];
Show[{ListPlot[d[1.45, 9.81, lösung[0.25, 0, 1.45, 0.01]], 
   PlotLegends -> {"Beschleunigung der Stabspitze"}], 
  Plot[9.81, {x, 0, 0.8}, PlotLegends -> {"Erdbeschleunigung"}, 
   PlotStyle -> Pink]}, 
 PlotLabel -> 
  "Fallbeschleunigung einer Punktmasse und Stabspitze im Vergleich", 
 AxesLabel -> {"t in s", "a in m/s^2"}, ImageSize -> Large]