read NEIL: Sipurd(3,6,15,n,N,20,1000,10); quit: