ring s = 0,(x,y),ds;
poly f = y2-2x28y-4x21y17+4x14y33-8x7y49+x56+20y65+4x49y16;
LIB "sing.lib";
LIB "classify.lib";
classify(f);

==>

The singularity is R-equivalent to A[2260]: y2+x2261 Milnor number(f) = 2260 modality(f) = 0