TOM: a Matlab code for the Numerical Solution of Boundary Value Problems