void f7 ( double t, double u[], double dudt[] ) { double g = 386.09, l = 10.0; dudt[0] = u[1]; dudt[1] = - ( g / l ) * sin ( u[0] ); return; }