image thumbnail

ball in a cone - Lagrange mechanics

version (330 KB) by Lucas Tassilo Scharbrodt
animate the motion of a mass point constrained to move on the wall of a circle cone


Updated 11 Jul 2018

View License

a point mass moves under the influence of gravity on the wall of a circle cone. Equations of motion for the two DOF's r and phi are obtained from the lagrangian L and solved numerically for a certain initial condition:

tspan = [0 T]; % time span for simulation
[r(t=0) r'(t=0) phi(t=0) phi'(t=0) ] initial conditions
y20 = [1.3 0 0 w ]; % w - angular frequency
f = @(l,y2) [y2(2); -g*cos(a) + y2(1)*(y2(4)^2)*((sin(a))^2)-k*(y2(2)^2+y2(4)^2)^0.5;y2(4);(-2*y2(2)*y2(4))/(y2(1))] ;

[l,y2]=ode45(f,tspan,y20); % call ode45 solver

the zip-file contains a mp4-video of the animation (created using matlabs WriteVideo() function)

Cite As

Lucas Tassilo Scharbrodt (2022). ball in a cone - Lagrange mechanics (, MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2017b
Compatible with any release
Platform Compatibility
Windows macOS Linux
Tags Add Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!