# “Fermat's last Python script”

## fermat.py

```
def fermat(n):
"""Returns triplets of the form x^n + y^n = z^n.
Warning! Untested with n > 2."""
from itertools import count
for x in count(1):
for y in range(1, x+1):
for z in range(1, x**n+y**n + 1):
if x**n + y**n == z**n:
yield x, y, z
```

changed September 16, 2007