This is a python implementation for non-linear wave simulations using the FDTD method. It is solved using the method of successive approximations.