o
    gx                     @   sH   d dl mZ d dlmZ d dlmZmZmZ d dlmZmZ dd Z	dS )    )symbols)dynamicsymbols)ReferenceFramePointParticle)LagrangesMethod
Lagrangianc                  C   sL  t d\} }t dd\}}td\}}}td}td}|d| |j }	|	|||j  |d||j }
|
|||j  td|	|}|| d  d |_td	|
|}|| | d  d |_|	| | |j f|
| | |j fg}t	|||}t
|| |f||d
}|  |jd }|||ksJ |jd }|||ksJ d S )Nzq1, q2   zm, k, nuNOP1pP1   pP2)	forcelistframer   )r   r   r   r   	locatenewxset_velr   potential_energyr   r   form_lagranges_equationseomdiff)q1q2q1dq2dmknur
   r   r   P2r   r   r   lagl_methodeq1eq2 r%   c/var/www/visachat/venv/lib/python3.10/site-packages/sympy/physics/mechanics/tests/test_lagrange2.pytest_lagrange_2forces	   s.   

r'   N)
sympyr   sympy.physics.mechanicsr   r   r   r   r   r   r'   r%   r%   r%   r&   <module>   s
    