Um Pouco de Teoria

5.1- Princípio de Fermat

5.2- Equações de Lagrange - Dinâmica

Programação

5.2 Equações de Lagrande - Dinâmica

kill(all)$

rx:0$   ry:y(t)$   rz:0$

vx:diff(rx,t)$

vy:diff(ry,t)$

vz:diff(rz,t)$

T:m*vx^2/2+ m*vy^2/2 + m*vz^2/2;

U:m*g*ry;

L:T-U;

/* este bloco deve ser calculado para cada coordenada generalizada */

q:y(t)$

t1:diff(L,q)$

t2:diff(L,diff(q,t))$

t3:diff(t2,t)$

Eq:t1-t3=0; 

solve(Eq, diff(q,t,2));