sylvester[n]:= if n=0 then sylvester[n]:2 else sylvester[n-1]^2-sylvester[n-1]+1$
/* Test cases */
makelist(sylvester[i],i,0,9);
apply("+",1/%);