Return to site

Rk4 Fortran Code

broken image


  1. Runge Kutta Fehlberg Fortran Code
  2. Rk4 Fortran Code Examples
In the following table, each line/entry contains the program file name, the page number where it can be found in the textbook, and a brief description. Click on the program name to display the source code, which can be downloaded.
Chapter 1: Introduction
first.f906-7First programming experiment
pi.f90 8Simple code to illustrate double precision
Chapter 2: Number Representation and Errors
oct.f90 49Print in octal format
hex.f90 50Print in hexadecimal format
numbers.f90 60-61Print internal machine representation of various numbers
xsinx.f90 77-79Example of carefully programming f(x) = x - sinx
Chapter 3: Locating Roots of Equations
bisection.f9094-95Bisection method
rec_bisection.f9095-96Recursive version of bisection method
newton.f90106-107Sample Newton method
secant.f90127-128Secant method
Chapter 4: Interpolation and Numerical Differentiation
coef.f90152-155Newton interpolation polynomial at equidistant pts
deriv.f90185-186Derivative by center differences/Richardson extrapolation
Chapter 5: Numerical Integration
sums.f90200Upper/lower sums experiment for an integral
trapezoid.f90207Trapezoid rule experiment for an integral
romberg.f90223-224 Romberg arrays for three separate functions
Chapter 6: More on Numerical Integration
rec_simpson.f90241Adaptive scheme for Simpson's rule
Chapter 7: Systems of Linear Equations
ngauss.f90270-271Naive Gaussian elimination to solve linear systems
gauss.f90285-287Gaussian elimination with scaled partial pivoting
tri.f90301-302Solves tridiagonal systems
penta.f90304Solves pentadiagonal linear systems
Chapter 8: More on Systems of Linear Equations
Chapter 9: Approximation by Spline Functions
spline1.f90385Interpolates table using a first-degree spline function
spline3.f90404-406Natural cubic spline function at equidistant points
bspline2.f90427-428Interpolates table using a quadratic B-spline function
schoenberg.f90430-431Interpolates table using Schoenberg's process
Chapter 10: Ordinary Differential Equations
euler.f90448-449Euler's method for solving an ODE
taylor.f90451Taylor series method (order 4) for solving an ODE
rk4.f90462-463Runge-Kutta method (order 4) for solving an IVP
rk45.f90472-473Runge-Kutta-Fehlberg method for solving an IVP
rk45ad.f90474Adaptive Runge-Kutta-Fehlberg method
Chapter 11: Systems of Ordinary Differential Equations
taylorsys1.f90489-490Taylor series method (order 4) for systems of ODEs
taylorsys2.f90491Taylor series method (order 4) for systems of ODEs
rk4sys.f90491-493,496Runge-Kutta method (order 4) for systems of ODEs
amrk.f90510-513 Adams-Moulton method for systems of ODEs
amrkad.f90513Adaptive Adams-Moulton method for systems of ODEs
Chapter 12: Smoothing of Data and the Method of Least Squares
Chapter 13: Monte Carlo Methods and Simulation
test_random.f90562-563Example to compute, store, and print random numbers
coarse_check.f90564Coarse check on the random-number generator
double_integral.f90574-575Volume of a complicated 3D region by Monte Carlo
volume_region.f90575-576Numerical value of integral over a 2D disk by Monte Carlo
cone.f90576-577 Ice cream cone example
loaded_die.f90581Loaded die problem simulation
birthday.f90583Birthday problem simulation
needle.f90584Buffon's needle problem simulation
two_die.f90585Two dice problem simulation
shielding.f90586-587Neutron shielding problem simulation
Chapter 14: Boundary Value Problems for Ordinary Differential Equations
bvp1.f90602-603Boundary value problem solved by discretization technique
bvp2.f90605-606Boundary value problem solved by shooting method
Chapter 15: Partial Differential Equations
parabolic1.f90618-619 Parabolic partial differential equation problem
parabolic2.f90620-621Parabolic PDE problem solved by Crank-Nicolson method
hyperbolic.f90633-634Hyperbolic PDE problem solved by discretization
seidel.f90642-645 Elliptic PDE solved by discretization/ Gauss-Seidel method
Chapter 16: Minimization of Functions
Chapter 17: Linear Programming
Runge kutta fehlberg fortran code

Addditional programs can be found at the textbook's anonymous ftp site:

Runge Kutta Fehlberg Fortran Code

[Home] [Features] [TOC] [Purchase] [Sample Code] [Web] [Manuals] [Errata] [Links]
Last updated: 5/20/2003

Rk4 Fortran Code Examples

You should compile your code for debugging. If you use gfortran, use -g -fbacktrace -fcheck=all -Wall. For ifort use -g -traceback-check -warn. Other compilers have similar options. – Vladimir F Dec 1 '14 at 19:30. Fourth Order Runge-Kutta. Intro; First Order; Second; Fourth; Printable; Contents Introduction. In the last section it was shown that using two estimates of the slope (i.e., Second Order Runge Kutta; using slopes at the beginning and midpoint of the time step, or using the slopes at the beginninng and end of the time step) gave an approximation with greater accuracy than using just a single.





broken image