Hola
Yo te explico, disculpame que me meta.
x + r*cos(phi)) esto es el valor de X más la proyección de r sobre ese mismo eje de coordenadas,
y + r*sin(phi)) es el valor de y más el valor de r proyectado sobre el otro eje.
En resumen, te da las coordenadas de un punto (i,j), son valores enteros por eso usa trunc para tomar solamente la parte entera de las ecuaciones anteriores.
Está usando coordenadas polares me parece si mal no recuerdo
Mirá
aquí
Saludos
Disculpame coso