On the approximate solutions of systems of odes by legendre operational matrix of differentiation