de := diff(theta(t), `$`(t, 2))+c*(diff(theta(t), t))+9.8*sin(theta(t))/L = 0; c := 2; L := 2;
init := theta(0) = 0, (D(theta))(0) = 4;
sol := dsolve({de, init}, theta(t), numeric, output=listprocedure);
Theta := subs(sol, theta(t));
plots[animate](plots[arrow],[[0,0],[L*sin(Theta(t)),-L*cos(Theta(t))],colour=red],t=0..5, axes=none,scaling=constrained);
these lines of codes were posted by robert Israel in one of his posts.