FUNCTION :   cuspmake1 -  set of inequivalent cusps for Gamma[1](N)

CALLING SEQUENCE :  cuspmake1(N)

PARAMETERS :      N - positive integer

GLOBAL VARIABLES : 

SYNOPSIS :   Returns a set of inequivalent cusps
             for Gamma[1](N).
             Each cusp a/c in the list is represented by [a,c]
             so that oo is represented by [1,0]

EXAMPLES :   

> cuspmake1(40);
{[0, 1], [1, 3], [1, 7], [1, 9], [1, 11], [1, 13], [1, 17], [1, 19], [1, 2], [1, 6], [1, 14],

    [1, 18], [1, 4], [3, 4], [1, 12], [7, 12], [1, 5], [3, 5], [2, 5], [4, 5], [1, 15], [13, 15],

    [7, 15], [4, 15], [1, 8], [3, 8], [7, 8], [5, 8], [1, 16], [3, 16], [7, 16], [13, 16], [1, 10],

    [3, 10], [7, 10], [9, 10], [1, 20], [3, 20], [7, 20], [9, 20], [1, 0], [3, 40], [7, 40],

    [9, 40], [11, 40], [13, 40], [17, 40], [19, 40]}

DISCUSSION :

SEE ALSO :