In this article, block BS methods are considered for the numerical solution of Volterra integro-differential equations (VIDEs). Convergence and stability properties are analyzed. A new Matlab code for the solution of VIDEs, called VIDEBS, is presented. Numerical results using a variable stepsize implementation show the effectiveness of the proposed code.