Il programma ricalca quello dell'albero di Pitagora, tranne che alterniamo le rotazioni e rispettive contrazioni a seconda che del numero dei passi.

Per la procedura useremo le stesse variabili, ma aggiungeremo una variabile tipo integer che segni i passi del programma, a partire da 1

Sottoprogramma PitNat(n: intero, x1, y1, x2, y2, x3, y3, x4, y4, ra, rb: Double, p: intero)
Se n = 0 allora

...   'usa lo stesso codice

Altrimenti

Se p è dispari allora

'per il quadrato a destra si usa la rotazione di 150° in senso antiorario intorno a (x4,y4) e la contrazione di 0.5, cioè

   T = 2.618
   ra = 0.5

'per il quadrato a sinistra si usa la rotazione di 240° in senso antiorario intorno a (x3,y3) e la contrazione di 0.866, cioè