Sorry for the late answer.

I still don't get along with those two files, but I doubt that I am allowed to post it in the internet.

I receive the same result in both files for a variable omega:

> Omega[1];
                                   3    2
                                   -- x1
                                   16    
i.e. 3/16*x1^2

but the command int(Omega[1],x1=0..1); in wholecode gives

> int(tempOmega[1], x1 = 0 .. 1);
                                   3    2
                                   -- x1
                                   16    
which is obviously wrong. While the same command in test.mw gives the correct result of 1/16.

Apparently Maple doesn't recognize the variable x1 and interprets it as a constant.

I still have no explanation for the restart thing neither.

map is fine :)

psi is defined as this:

psi := <v(x1, x2, x3), w(x1, x2, x3), u(x1, x2, x3)>

where v,w,u are vectors containing functions of x1,x2,x3 in its entries:

t[1]:=unapply(E1[13],x1,x2,x3):

E1[i] is something like this: x1+x2+x3 or x1/4+x2/2....

I'd still be interested to find out where the problem came from, as it wasn't the first time that it came up

i'm just wondering how I can produce a Matrix of vectors v(-x1, x2, x3), -w(-x1, x2, x3), u(x1, x2, x3)?

v(-x1,x2,x3) gives a correct result but if i call <v(-x1, x2, x3), -w(-x1, x2, x3), u(x1, x2, x3)> i receive an error message : Error, (in Matrix) this entry is too wide or too narrow: ....

As I don't have that problem with psi := <v(x1, x2, x3), w(x1, x2, x3), u(x1, x2, x3)>, i assume that it comes from "-"

I found out that the best way to solve the problem was to store the piecewise function in an array and modify the integralfunction. I wrote a script which now integrates the entries of the matrix piecewise. 

like this:

int(int(int(T1[11], x1 = -x2+1 .. x2), x2 = (1/2)*1 .. -x3+1), x3 = 0 .. (1/2)*1)+
 int(int(int(T1[12], x2 = x1 .. -x1+1), x1 = x3 .. (1/2)*1), x3 = 0 .. (1/2)*1)+
 int(int(int(T1[9], x1 = x2 .. -x2+1), x2 = x3 .. (1/2)*1), x3 = 0 .. (1/2)*1)+
 int(int(int(T1[10], x2 = -x1+1 .. x1), x1 = (1/2)*1 .. -x3+1), x3 = 0 .. (1/2)*1)+

Hopefully my comment helps users having the same problem.

Cheerio,

Nerdpole

Thanks! This does the job!