def test_solve_fixed_trajectory():
 (t0, k0) = (0, np.array([5.0]))
 results = _compute_fixed_length_solns(model, t0, k0)
 for (integrator, numeric_solution) in results.items():
  ti = numeric_solution[:, 0]
  analytic_solution = solow_analytic_solution(ti, k0, *valid_params)
  np.testing.assert_allclose(numeric_solution, analytic_solution)

def test_solve_fixed_trajectory ( ) : ( t0 , k0 ) = ( 0 , np . array ( [ 5.0 ] ) ) results = _compute_fixed_length_solns ( model , t0 , k0 ) for ( integrator , numeric_solution ) in results . items ( ) : ti = numeric_solution [ : , 0 ] analytic_solution = solow_analytic_solution ( ti , k0 , * valid_params ) np . testing . assert_allclose ( numeric_solution , analytic_solution )
testing computation of fixed length solution trajectory .