Solving the FULL! System on a square using fix-point iterations, plotting quivers for the velocity vector field and plotting the pressure