[Scilab-users] ode with algebraic constraint

[Jens Simon Strom]

*Hallo,**
**To solve the equation of motion for a mass point m in a constant force 
field f=[fx; fy; fz] one can use ode(...) with the function**
**
**function dz=EoM(t, z, m,f)//z=[x; y; z;  vx; vy; vz] (6 x 1)**
**   dz(1:3)=z(4:6)**
**   dz(4:6)=f/m**
**endfunction**
**
**If the motion shall be constrained to a sphere of constant radius R - 
is there a way to implement this constraint into the function. **
**
**Of course spherical coordinates may be a better choice, however I am 
looking for a cartesion solution.**
**
**Kind regards **
**Jens**
*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180810/7c05dd40/attachment.htm>