[ create a new paste ] login | about

Link: http://codepad.org/PW2a9CmW    [ raw code | output | fork ]

C, pasted on May 14:
#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 */


Output:
x = 0.000000, y = 0.064493
x = 1.000000, y = 0.124827
x = 2.000000, y = 0.181269
x = 3.000000, y = 0.234072
x = 4.000000, y = 0.283469
x = 5.000000, y = 0.329680
x = 6.000000, y = 0.372911
x = 7.000000, y = 0.413354
x = 8.000000, y = 0.451188
x = 9.000000, y = 0.486583
x = 10.000000, y = 0.519695
x = 11.000000, y = 0.550671
x = 12.000000, y = 0.579650
x = 13.000000, y = 0.606759
x = 14.000000, y = 0.632120
x = 15.000000, y = 0.655846
x = 16.000000, y = 0.678042
x = 17.000000, y = 0.698806
x = 18.000000, y = 0.718231
x = 19.000000, y = 0.736403
x = 20.000000, y = 0.753403
x = 21.000000, y = 0.769307
x = 22.000000, y = 0.784185
x = 23.000000, y = 0.798103
x = 24.000000, y = 0.811124
x = 25.000000, y = 0.823306
x = 26.000000, y = 0.834701
x = 27.000000, y = 0.845362
x = 28.000000, y = 0.855335
x = 29.000000, y = 0.864665
x = 30.000000, y = 0.873393
x = 31.000000, y = 0.881558
x = 32.000000, y = 0.889197
x = 33.000000, y = 0.896343
x = 34.000000, y = 0.903028
x = 35.000000, y = 0.909282
x = 36.000000, y = 0.915133
x = 37.000000, y = 0.920606
x = 38.000000, y = 0.925726
x = 39.000000, y = 0.930517
x = 40.000000, y = 0.934998
x = 41.000000, y = 0.939190
x = 42.000000, y = 0.943112
x = 43.000000, y = 0.946781
x = 44.000000, y = 0.950213
x = 45.000000, y = 0.953424
x = 46.000000, y = 0.956428
x = 47.000000, y = 0.959238
x = 48.000000, y = 0.961867
x = 49.000000, y = 0.964326
x = 50.000000, y = 0.966627
x = 51.000000, y = 0.968779
x = 52.000000, y = 0.970793
x = 53.000000, y = 0.972676
x = 54.000000, y = 0.974438
x = 55.000000, y = 0.976087
x = 56.000000, y = 0.977629
x = 57.000000, y = 0.979072
x = 58.000000, y = 0.980422
x = 59.000000, y = 0.981684
x = 60.000000, y = 0.982866
x = 61.000000, y = 0.983971
x = 62.000000, y = 0.985004
x = 63.000000, y = 0.985972
x = 64.000000, y = 0.986876
x = 65.000000, y = 0.987723
x = 66.000000, y = 0.988514
x = 67.000000, y = 0.989255
x = 68.000000, y = 0.989948
x = 69.000000, y = 0.990596
x = 70.000000, y = 0.991203
x = 71.000000, y = 0.991770
x = 72.000000, y = 0.992301
x = 73.000000, y = 0.992798
x = 74.000000, y = 0.993262
x = 75.000000, y = 0.993697
x = 76.000000, y = 0.994103
x = 77.000000, y = 0.994483
x = 78.000000, y = 0.994839
x = 79.000000, y = 0.995172
x = 80.000000, y = 0.995483
x = 81.000000, y = 0.995775
x = 82.000000, y = 0.996047
x = 83.000000, y = 0.996302
x = 84.000000, y = 0.996541
x = 85.000000, y = 0.996764
x = 86.000000, y = 0.996972
x = 87.000000, y = 0.997168
x = 88.000000, y = 0.997350
x = 89.000000, y = 0.997521
x = 90.000000, y = 0.997681
x = 91.000000, y = 0.997831
x = 92.000000, y = 0.997971
x = 93.000000, y = 0.998101
x = 94.000000, y = 0.998224
x = 95.000000, y = 0.998338
x = 96.000000, y = 0.998446
x = 97.000000, y = 0.998546
x = 98.000000, y = 0.998640
x = 99.000000, y = 0.998727


Create a new paste based on this one


Comments: