#include <chplot.h> #include <numeric.h> #define NVAR 2 #define POINTS 256 void func(double t, double y[], double dydt[], void *param) { double mu; mu = *(double*)param; dydt[0] = y[1]; dydt[1]=mu*(1-y[0]*y[0])*y[1] - y[0]; } int main() { double t0=1, tf=30, y0[NVAR] = {1, 0}; double t[POINTS], y[NVAR][POINTS]; double mu = 2; oderk(func, t0, tf, y0, &mu, t, y); plotxy(t, y, "The solution for the van der Pol equation", "t (seconds)", "y1 and y2"); }

---