Fuller's Problem (Muscod)

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

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

rhs[0] = xd[1];
rhs[1] = myu;
rhs[2] = xd[0]*xd[0];