Summary: The numerical solution of problems with singularities presents special difficulties for most methods. Adjustments to standard methods are typically made for only a special type of singularity, usually known a priori. The family of sine numerical methods is naturally suited for general singular problems. Here, the methods are extended and applied to two-dimensional, elliptic first-order systems of mixed boundary value problems with singularities of the form \(x^\alpha\).