The Tau method for numerical (approximate) solutions of differential equations