<<sm.mat;							(*In[1]*)
V=SeifertMatrix[2,{1,1,1}]					(*In[2]*)
V=SeifertMatrix[3,{1,2,1,2,1,2,1,2}];				(*In[3]*)
Factor[Det[ V - t*Transpose[V] ]]				(*In[4]*)
De={1,2,3,4,5,6,7,1,2,3,4,5,6,1,2,3,4,5,1,2,3,4,1,2,3,1,2,1};	(*In[5]*)
b1={-4,-3,-5,-4}; b2={-6,-7,4,5}; bb2={-5,-4,6};		(*In[6]*)
b3[3]={-3,4,-3,-4,3}; b3[5]={-5,4,-5,-4,5};			(*In[7]*)
b=Join[b1,{-5,-5,-5,-5,-5,3,3,3,3,3,3,3,-4,-5},b2,
Table[-6,{14}],bb2,{5,-4,-4},b3[5],De];	V=SeifertMatrix[8,b];	(*In[8]*)
Factor[Det[ V - t*Transpose[V] ]]				(*In[9]*)
ssmW[2,{{n,1}}]							(*In[10]*)
W=ssmW[2,{1,1,{n-2,1}}]						(*In[11]*)
Det[W]								(*In[12]*)
alpha1=13; alpha2=6;						(*In[13]*)
b=Join[b1,{{1-e1,5},{e1,3}},b2,{{1-alpha1,6}},bb2,
{{-e2,5},{1+e2,3},{-beta1,4}},b3[3],{{1-beta2,4}},De];		(*In[14]*)
W=ssmW[8,b];							(*In[15]*)
d=Det[W//.{beta1->3,beta2->3,e1->12-alpha2-e2}]/16		(*In[16]*)
d=Det[W//.{beta1->1,beta2->5,e1->12-alpha2-e2}]/16		(*In[17]*)
Union[Table[ Mod[d//.e2->n,25], {n,25} ]]			(*In[18]*)
Union[Table[Mod[n*n,25],{n,25}],Table[Mod[-n*n,25],{n,25}]]	(*In[19]*)