Question: How to solve motion on curve?

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?

Please Wait...