We present a novel iterative scheme for restoring uneven illumination in grayscale images. Our approach solves, at each global iteration, a nonlinear elliptic equation for an auxiliary field u and then updates the illumination via an explicit time-marching step driven by the logarithmic potential of u. We establish existence and uniqueness of the discrete subproblems, analyze convergence of the fixed point iteration, and demonstrate performance on standard test images, reporting PSNR and SSIM gains over classical methods. Future work will address color coupling and GPU acceleration.