FORTRAN 90 - DVERK? - initial condition solver

Click For Summary
FORTRAN users seeking a routine similar to MATLAB's ODE45 for solving initial condition problems in ODEs are advised to consider coding their own solutions. Simple methods like Euler's method and 4th order Runge-Kutta are recommended for their ease of implementation and accuracy. However, both DVERK and ODE45 utilize adaptive step-size techniques, making them more complex to replicate. A reference to RKSUITE is provided as a potential resource for further exploration. Custom coding may be necessary for those requiring specific functionalities in FORTRAN 90.
FrogPad
Messages
801
Reaction score
0
Hi,

I am new to FORTRAN and am trying to work with a system of ODE's. Does someone know of a routine like DVERK for solving initial condition problems for FORTRAN 90?

In MATLAB I would simply use ODE45. I just need something similar.

Thank you!
 
Mathematics news on Phys.org
Your best bet would be to code your own. Start with a Eulers method, this is very simple though not so precise. A 4th order Runga Kutta routine is also easy to code and is quite accurate. Do a bit of web searching and have a go at it.
 
Yeah both the Euler and the fixed step size Runge-Kutta are simple implementations. The DVERK and ode45 both use adaptive step-size routines, and it gets much more complicated there.
 

Thread 'Area of Overlapping Squares'

Here is a little puzzle from the book 100 Geometric Games by Pierre Berloquin. The side of a small square is one meter long and the side of a larger square one and a half meters long. One vertex of the large square is at the center of the small square. The side of the large square cuts two sides of the small square into one- third parts and two-thirds parts. What is the area where the squares overlap?
View full post »

Similar threads

Fortran FORTRAN 90 - Insert Matrices into Other Matrices
  • · Replies 2 ·
Replies
2
Views
1K
Fortran programming to solve ODEs
  • · Replies 1 ·
Replies
1
Views
4K
Large-scale eigen solver package for Fortran?
  • · Replies 1 ·
Replies
1
Views
2K
Physics Scattering / inverse scattering industry jobs?
  • · Replies 3 ·
Replies
3
Views
3K
Fortran Solving Magnetic Field Using Floating Random Walk
  • · Replies 6 ·
Replies
6
Views
2K
How to read a column of data into Fortran without arrays? (Fortran 77)
  • · Replies 7 ·
Replies
7
Views
3K
Graduate How to work around this equation giving infinities in the numerical calculation?
  • · Replies 4 ·
Replies
4
Views
4K
Graduate Forcing conditions in a solution of an under-determined system
  • · Replies 2 ·
Replies
2
Views
2K
Undergrad Can initial conditions for an ODE be given by functions instead of constants?
  • · Replies 10 ·
Replies
10
Views
2K
Undergrad Make Intelligent Initial Guesses for Newton-Raphson Programmatically
  • · Replies 16 ·
Replies
16
Views
4K