Question: Animation optimization (phase portrait)

hi! lately I'm thinking about animating complex phase portraits so I tried something like this:

restart;
with(plots);
f := (a, z) -> (z*a - 1)/(z^2 + z + 1);
display(seq(densityplot(proc(x, y) local res; res := evalhf(argument(f(1/100*a, y + x*I))); return res; end proc, -4 .. 4, -4 .. 4, axes = boxed, view = [-2 .. 2, -2 .. 2], restricttoranges, colorstyle = HUE, style = surface, grid = [201, 201]), a = 1 .. 200), insequence = true);
 

but this takes a long time to load. Do you how can we optimizate this code?

Please Wait...