In this paper we provide algorithms for computing the bidiagonal decomposition of Wronskian matrices monomial basis polynomials and exponential polynomials. It is also shown that these can be used to perform accurately some algebraic computations with matrices, such as calculation their inverses, eigenvalues or singular values solutions linear systems. Numerical experiments illustrate results.