alpha := vector([cos(theta), sin(theta), sin(k*theta)]);
alpha1 := map(diff, alpha, theta);
alpha11 := map(diff, alpha1, theta);
with(linalg):
cr := multiply(alpha1, alpha11)/simplify(multiply(alpha1, alpha1));
F := vector([0, 0, -m*g]);
c := multiply(F, alpha1/simplify(multiply(alpha1, alpha1)))/m;
motion := diff(f(theta), theta$2) + cr*diff(f(theta), theta)^2 = c;
pdsolve and dsolve not work, how to solve above motion equation with maple?