Applications of differential equations to characterize the base of warped product submanifolds of cosymplectic space forms