Motzkin numbers: Difference between revisions
add parigp
(add parigp) |
|||
Line 812:
40 66_368_199_913_921_497 false
41 192_137_918_101_841_817 false</pre>
=={{header|PARI/GP}}==
<lang parigp>
M=vector(41)
M[1]=1
M[2]=2
for(n=3, 41, M[n]=(M[n-1]*(2*n+1) + M[n-2]*(3*n-3))/(n+2))
M=concat([1], M)
for(n=1, 42, print(M[n]," ",isprime(M[n])))
</lang>
=={{header|Perl}}==
| |||