Surface Complex function

using GLMakie
let
    x = -2:0.005:2
    y = -2:0.005:2
    f(z) = (z^2 + 1)/(z^2 - 1)
    fvals = [f(u+1im*v) for u in x, v in y]
    fvalues = abs.(fvals)
    fargs = angle.(fvals)
    indxCut = fvalues .> 3
    fvalues[indxCut] .= 3.01
    fig, ax, pltobj = surface(x, y, fvalues, color = fargs,
        colormap = :diverging_rainbow_bgymr_45_85_c67_n256,
        colorrange = (-π,π), backlight = 1f0, highclip = :black)
    fig
end
Dependencies
      Status `~/Desktop/BeautifulMakie/_assets/Project.toml`
  [e9467ef8] GLMakie v0.4.4

Last Updated: July 25, 2021