00001 void myAcadoDifferentialEquation( double *x, double *f, void *user_data ){ 00002 00003 // x[0] -> time t 00004 // x[1] -> v 00005 // x[2] -> s 00006 // x[3] -> m 00007 // x[4] -> u 00008 00009 double h = 0.01; 00010 00011 double t = x[0]; 00012 double v = x[1]; 00013 double s = x[2]; 00014 double m = x[3]; 00015 double u = x[4]; 00016 00017 f[0] = s + h*v; 00018 f[0] = v + h*(u-0.02*v*v)/m; 00019 f[0] = m - h*0.01*u*u; 00020 00021 }