![]() h takes the role of a small time step dt. Note that f is time independent, therefore has only 1 argument. Where f is the function that constitutes the equations of motion (or equivalently the RHS of my system of diff. The model I use is implemented as follows: k1 = h * f(vec_old,omega,phi) that describe the equations of motions, from which I want to calculate the trajectory of the ball, and compare results for different orders of the RK method. ![]() I have a 6-dimensional vector of position and velocity of a ball, and a corresponding system of diff. I don't know how much of the code I should show, so I'll describe the problem in detail, and please guide me as to what I should add/remove to/from the post to make it answerable. ![]() I'm trying to implement the Runge Kutta method in Fortran and am facing a convergence problem. ![]()
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |