This thesis deals with the application of wavelet bases for the numerical solution of operator equations, as boundary value problems and boundary integral equations. The use of suitable wavelet bases has the advantage that the arising stiffness matrices are well-conditioned uniformly in their sizes, allowing for an efficient iterative solution, and, for integral equations, that they allow for a compression of the in that case densely populated stiffness matrices to truly sparse ones without reducing the convergence rates.