sbs->sbs4 = sbsNew4(ta);
sbs->sbs4 = sbsNew4(treeaClone(ta));
struct SpringBallSystem *sbsNew(struct TreeAdaptor *ta) {  struct SpringBallSystem *sbs;  int howManyNodes = treeaNodeCount(ta);  sbs = clCalloc(sizeof(*sbs), 1);  sbs->sbs4 = sbsNew4(ta);  sbs->sys.function = func;  sbs->sys.jacobian = NULL;  sbs->sys.dimension = howManyNodes * 2 * sbs->sbs4->d;  sbs->sys.params = sbs->sbs4;  sbs->modelSpeed = 1.0;  sbs->t = sbs->modelSpeed * cldatetimeStaticTimer();  sbs->T = (gsl_odeiv_step_type *) gsl_odeiv_step_rk4;  //sbs->T = (gsl_odeiv_step_type *) gsl_odeiv_step_bsimp;  sbs->s = gsl_odeiv_step_alloc(sbs->T, sbs->sys.dimension);  sbs->c = gsl_odeiv_control_y_new(1e-1, 0.0);  sbs->e = gsl_odeiv_evolve_alloc(sbs->sys.dimension);  sbs->h = 0.1;  return sbs;}