Solution of IVP using eigenfunction expansions