codepad
[
create a new paste
]
login
|
about
Language:
C
C++
D
Haskell
Lua
OCaml
PHP
Perl
Plain Text
Python
Ruby
Scheme
Tcl
#include <stdio.h> double f(double x, double y) { return (1.0 - y) / 15; } #define h (double)1.0 #define A (double)0.0 #define B (double)100.0 #define C (double)0.0 int main() { double x, y, k1, k2, k3, k4, dy; y = C; for (x = A; x < B; x += h) { k1 = f(x, y) * h; k2 = f(x + h / 2.0, y + k1 / 2.0) * h; k3 = f(x + h / 2.0, y + k2 / 2.0) * h; k4 = f(x + h, y + k3) * h; dy = (k1 + 2.0 * k2 + 2.0 * k3 + k4) / 6.0; y += dy; printf("x = %lf, y = %lf\n", x, y); } return 0; } /* end */
Private
[
?
]
Run code
Submit