Describir: An efficient approximation method for nonlinear singular value problems arising in astrophysics: an operational matrix approach