A new high accuracy finite difference discretization for the solution of 2D nonlinear biharmonic equations using coupled approach