10進basicで分割数を求める

mathematician.blog.jp行列式で分割数を求めています詳しくは上のサイトから

INPUT   n
DIM  a(200,200)
MAT  a=ZER
DIM   d(n,n)

FOR    i=1  TO   n
   FOR  j=1  TO   n
      FOR  k=1  TO  30
         IF  i=j    THEN   LET a(i,j)=1  
         IF    j=i+1  THEN    LET a(i,j)=-1
         IF   i=j+(k*(3*k-1)/2-1)  THEN   LET a(i,j)=(-1)^(k+1)
         IF   i=j+(k*(3*k+1)/2-1)  THEN   LET a(i,j)=(-1)^(k+1)
      NEXT K
       
   NEXT J
NEXT I
REM  -------------------------------------------
FOR    b=1   TO   n
   FOR   c=1  TO n
      LET d(b,c)=  a(b,c)
      PRINT  d(b,c);
   NEXT C
   PRINT
NEXT B
REM-----------------------------------
PRINT   DET(d)

END