This paper studies the nonlinear two-dimensional difference system with multiple delaysΔ(xn+p1nxn-τ1)+f1(n,xa1n,…,xahn,yb1n,…,ybkn)=q1n,Δ(yn+p2nyn-τ2)+f2(n,xc1n,…,xchn,yd1n,…,ydkn)=q2n,n≥n0. Using the Banach fixed point theorem and a few new analysis techniques, we show the existence of uncountably many bounded positive solutions for the system, suggest Mann iterative algorithms with errors, and discuss the error estimates between the positive solutions and iterative sequences generated by the Mann iterative algorithms. Examples to illustrates the results are included.