PROGRAM MAPEL 18

JULIASET with(plots); JuliaSet := proc (a, b) local z1, z2, z1s, z2s, m; z1, z2 := a, b; z1s := z1^2; z2s := z2^2; for m to 30 while z1s+z2s < 4 do z1, z2 := z1s-z2s, 2*z1*z2+.55; z1s := z1^2; z2s := z2^2 end do; m end proc; densityplot(JuliaSet, -1.5 .. 1.5, -1.4 .. 1.4, colorstyle = HUE, grid = [200, 200], style = patchnogrid, axes = none);