matantihadamard(n)={ matrix(n,n,i,j, my(t=j-i+1); if(t<1,t%2,t<3) ); } fib(n)=matdet(matantihadamard(n))