calc(N,Sum,Sum):- 2*(N-1)+1>1001, !. calc(N,Sum,Result):- Sum1 is Sum+(4*(2*N-1)^2-6*((N-1)*2)), N1 is N+1, calc(N1,Sum1,Result). main:-calc(2,1,Ans),write(Ans).