Integrating a differential-functional equation

I need to solve a differential-functional equation:

$\partial_t x_m(t,s) = \sum_n A_{mn} x_n(t,s) + \int_0^t \sum_n \sum_{n'} B_{mnn'}(t,s') x_n(t,s) x_{n'}(t,s') ds'$

with $t > s$ and initial condition $x_m(t,t) = C_m$.

What is the best numerical procedure?

-

Right, sorry. So your equation just doesn't make any sense, because there is no constraint on the s-evolution of x. Think first of the simplest case : in one variable, A = 1, B = 0, $\partial_t x(t,s) = x(t,s)$. Impose an equation for $\partial_s x$ or a constraint. –  Antoine Levitt Dec 4 '11 at 12:00
I have solved it in a trivial way: using Euler discretization in the $t$ direction and trapezoid quadrature in $s$ direction, with equal steps. Works quite well for the data I use.