Polynomiography is a fusion of Mathematics and Art, which as a software results in a new form of abstract art. Rendered images are through algorithmic visualization of solving a polynomial equation via iteration schemes. Images are beautiful and diverse, yet unique. In short, polynomiography allows us to draw unique and complex-patterned images of polynomials which be re-colored in different ways through different iteration schemes. In the modern age, polynomiography covers a variety of applications in different fields of art and science. The aim of this paper is to present polynomiography using newly constructed root-finding algorithms for the solution of non-linear equations. The constructed algorithms are two-step predictor corrector methods. For reducing computational cost and making the algorithm more effective, we approximate the second derivative via interpolation technique. These methods have been derived by employing Househölder's method, interpolation technique and Taylor's series expansion. The convergence criterion of the newly developed algorithms has been discussed and proved their sixth-order convergence which is higher than many existing algorithms. To analyze the accuracy, validity and applicability of the proposed methods, several arbitrary and engineering problems have been tested and the obtained numerical results certify the better efficiency of the suggested methods against the other well-known iteration schemes given in the literature. Finally, we present polynomiography through the constructed iteration schemes and give a detailed comparison with the other iteration schemes which reflects the convergence properties and graphical aspects of the constructed algorithms.