FORTRAN 90 - DVERK? - initial condition solver

AI Thread 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 'Fermat's Last Theorem'

Fermat's Last Theorem has long been one of the most famous mathematical problems, and is now one of the most famous theorems. It simply states that the equation $$ a^n+b^n=c^n $$ has no solutions with positive integers if ##n>2.## It was named after Pierre de Fermat (1607-1665). The problem itself stems from the book Arithmetica by Diophantus of Alexandria. It gained popularity because Fermat noted in his copy "Cubum autem in duos cubos, aut quadratoquadratum in duos quadratoquadratos, et...
View full post »

Thread 'Non-orthogonal bases'

Suppose ,instead of the usual x,y coordinate system with an I basis vector along the x -axis and a corresponding j basis vector along the y-axis we instead have a different pair of basis vectors ,call them e and f along their respective axes. I have seen that this is an important subject in maths My question is what physical applications does such a model apply to? I am asking here because I have devoted quite a lot of time in the past to understanding convectors and the dual...
View full post »
Thread 'Imaginary Pythagorus'

Thread 'Imaginary Pythagorus'

I posted this in the Lame Math thread, but it's got me thinking. Is there any validity to this? Or is it really just a mathematical trick? Naively, I see that i2 + plus 12 does equal zero2. But does this have a meaning? I know one can treat the imaginary number line as just another axis like the reals, but does that mean this does represent a triangle in the complex plane with a hypotenuse of length zero? Ibix offered a rendering of the diagram using what I assume is matrix* notation...
View full post »

Similar threads

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

Hot Threads

Recent Insights

Back
Top