/* $Id: nc3.c 1208 2006-12-30 02:22:56Z SYSTEM $ */

/* Numeral Computing with gnuplot: 3 */

#include <stdio.h>

int main(void);

int main(void) {
  int n;            /* ステップ数 */
  double x;         /* ステップ n における x の値 */
  double y;         /* x における y の値 */
  double dx = 0.1;  /* ステップ毎のxの増減*/

  /* 最初の値 f(-10) を求める */
  x = -10.0;
  y = (x*x*x)/2 - x*x + 20*x;
  printf("%f %f\n", x, y);
  for (n=1; n<=200; n++) {
    y = y + (3*dx*x*x)/2 - 2*dx*x + 20*dx;
    x = x + dx;
    printf("%f %f\n", x, y);
  }

  return 0;
}