Fuller's Problem (C)

The differential equations for Fuller's problem in C code read as follows

double myu = -1 + 2*u[0];
 
rhs[0] = xd[1];
rhs[1] = myu;
rhs[2] = xd[0]*xd[0];