Can Maple 13 solve nonlinear algebraic equations system numerically without using fsolve? for example Gauss siedel iterative method .

I wrote a procedure that solves a volterra integral equation numerically. This procedure  takes N , number of radial basis functions, as parameter  .the problem is:Error,  bad index into Matrix

How can I fix this problem?

here is my procedure:

I wrote a procedure that solves an optimal control problem governed by volterra integral equation. This procedure  takes N , number of lagrange interpolating polynomials, as parameter and call NLPSolve function .the problem is: for N>10 maple displays this Error: Error, (in Optimization:-NLPSolve) object too large. How can I fix this problem?

here is my procedure:

I am tring to solve a PDE system with four equations that its unknowns have three independent variable.

Is it possible to solve this system numerically?