A question: Mechanical system with white noise external forces

[Peter230655]

I have a mechanical system described as d/dt (q) = f(q, t, F), where q …generalized coordinates, F…external forces to be modeled as white noise (Where I used sympy.physics.mechanics to generate these equations of motion)

What I did in order to apply SDEINT, I linearized about the forces F: B = f(q, t, F).diff(F).subs({F: 0.}) F1 = f(q, t, F).subs({F: 0})

Then I calculated as Result = itoint(F1, B, x0, tspan)

Is this the right approach, or what would be the correct one?

Thanks for any help!

[Peter230655]

Just a reminder, is my idea correct? Thanks a lot for your help!

[Peter230655]

Has sdeint been 'abandoned'? Would be a shame, it is so nice and useful!