We consider a class of equations in divergence form with a singular/degenerate weight −div(|y| a A(x, y)∇u) = |y| a f (x, y) or div(|y| a F (x, y)) . Under suitable regularity assumptions for the matrix A and f (resp. F ) we prove Hölder continuity of solutions which are even in y ∈ R, and possibly of their derivatives up to order two or more (Schauder estimates). In addition, we show stability of the C 0,α and C 1,α a priori bounds for approximating problems in the formas ε → 0. Finally, we derive C 0,α and C 1,α bounds for inhomogenous Neumann boundary problems as well. Our method is based upon blow-up and appropriate Liouville type theorems.