? vector(30,j, matsolvemod(matrix(j,1,x,y,1), primes(j)~, vector(j,i,i-1)~)[1])

%55 = [0, -2, -8, 52, -788, -788, -210998, 4383592, -34415168, -703693778, 5765999452, 2211931390882, -138782093170508, -5006786309605868, -253579251611336438, -12551374903381164638, -142908008812141343558, 40235059344426324076912, -2891974474640747950504838, 169991099649125127278835142, -18242056294531940584645872728, -1443780877268247785856392194178, -85102544828125743391232590601558, -9165296078263504689482449202356418, 941584379775558526136539054851975982, -109725677889609325862012364017072315378, -1972624592757588213062518899061131219938, 861481022448550626004372260093860171043622, 67587260079918535402473915458935723029299842, 8739372161433421556929034933493037822954213972]

? vector(30,j, matsolvemod(matrix(j,1,x,y,1), primes(j)~, vector(j,i,i)~)[1])

%56 = [1, -1, -7, 53, -787, -787, -210997, 4383593, -34415167, -703693777, 5765999453, 2211931390883, -138782093170507, -5006786309605867, -253579251611336437, -12551374903381164637, -142908008812141343557, 40235059344426324076913, -2891974474640747950504837, 169991099649125127278835143, -18242056294531940584645872727, -1443780877268247785856392194177, -85102544828125743391232590601557, -9165296078263504689482449202356417, 941584379775558526136539054851975983, -109725677889609325862012364017072315377, -1972624592757588213062518899061131219937, 861481022448550626004372260093860171043623, 67587260079918535402473915458935723029299843, 8739372161433421556929034933493037822954213973]