FORTRAN 90 - DVERK? - initial condition solver

Click For Summary
SUMMARY

This discussion focuses on finding a FORTRAN 90 routine similar to MATLAB's ODE45 for solving initial condition problems in ordinary differential equations (ODEs). Users suggest that while DVERK is not directly available, implementing basic methods like Euler's method or a 4th order Runge-Kutta routine is feasible for beginners. For more advanced adaptive step-size solutions akin to DVERK and ODE45, the RKSUITE library is recommended as a valuable resource.

PREREQUISITES
  • Basic understanding of ordinary differential equations (ODEs)
  • Familiarity with FORTRAN 90 programming
  • Knowledge of numerical methods, specifically Euler's method and Runge-Kutta methods
  • Experience with adaptive step-size algorithms
NEXT STEPS
  • Research the implementation of Euler's method in FORTRAN 90
  • Learn how to code a 4th order Runge-Kutta routine in FORTRAN 90
  • Explore the RKSUITE library for adaptive step-size ODE solvers
  • Investigate the differences between fixed and adaptive step-size methods in numerical analysis
USEFUL FOR

This discussion is beneficial for FORTRAN developers, numerical analysts, and anyone interested in solving initial condition problems for ODEs using 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.
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 6 ·
Replies
6
Views
2K
Replies
7
Views
3K
Replies
4
Views
4K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
10
Views
2K
  • · Replies 16 ·
Replies
16
Views
4K