def sphere_conserved ( t, y ): #*****************************************************************************80 # ## sphere_conserved() evaluates a conserved for the sphere ODE. # import numpy as np h = np.sum ( y**2, axis = 1 ) return h