equation icon  un+1 = un - Dt-(^f(un ,un ,...,un ) - ^f(un ,un ,...,un )), (97) j j Dx j-r j-r+1 j+q j-r-1 j- r j+q-1