An algorithm that exploits features of the bicomplex field to compute 4-dimensional Newton fractals for procedural texturing applications. The generated fractals are computed in a pixel shader only on the target surface to achieve real-time performance.