u is dimensioned as array[1..mp,1..np]
however in the code,

Have
For i := 1 to Ndata do
...
For j:= 1 to ma do
u[i,j] := etc.

This sends the first index out of bounds since ndata>mp.
Looks like indexing is wrong way round.