#include #include /* euler.c : 1999 3 28 by T.Oguni */ #define NF 200 double f(double t, double y) {return y-20.; } euler(tvals,yvals,t0,tf,y0,h) double tvals[NF+1], yvals[NF+1],t0,tf,y0,h; { double y, t, y1, t1, s; int i, j, n; y=y0; t=t0; yvals[0]=y0; tvals[0]=t0; s=(tf-t0)/h; n=(int)s+1; for (j=1; j