@Hummingbird Please attach your worksheet (rather than cut-and-paste) using the green up arrow, which is the last icon on the second row of the toolbar in the MaplePrimes editor.

@mehdi jafari Mehdi, the range is already set by the boundary conditions to -50..50. The range option is only needed when the boundary conditions do not give both endpoints of the space range.

@Hummingbird

In order for

solve({seq(eq[i], i= 1..6)}, convert(phi4N, set));

to work, the eq[i] must be those that you created in the for loop, not those that you entered manually.

The code

solve({seq(eq[i], i= 1..6)}, convert(phi46N, set));

is wrong because the second argument only uses one of the six variables. If you want to solve for fewer variables than there are equations, use eliminate instead of solve.

@Kitonum Clearly the OP did not actually use that .. syntax. There would have been an error message if it had been input like that.

Please show explicitly the code that does not work.

@sarra Remove the range option! That was not in your code when I made my recommendation. Your range is already set at -50..50 by your boundary conditions. Please don't change your original code and then say that my fix doesn't work!

@mahmood180 Please post this as a seperate Question. It seems unrelated to the original.

@c4meleon Change the dsolve command to

res:=dsolve({diff(y(x),x)=f(x,y(x)),y(0)=1}, {y(x)});

You must be using a significantly older version of Maple. Which is it?

@c4meleon Please show your dsolve command that is not working.

@Preben Alsholm There have been several versions of this problem posted by the OP in recent weeks, with various amounts of specificity as to what the symbols mean. The suggestion that A is a linear operator with spectral radius less than 1 was put forth by me. However, looking at the other versions of this problem, it is clear that A is not a linear operator. In one version, A is given by a procedure. The f's are NxN matrices, and A operates on some elements of interior of f, replacing them with the average of the four horizontally and vertically adjacent elements.

@Preben Alsholm My square brackets after D were unintentional. I'm also surprised that they worked. Thanks for spotting that. I'll fix the original.

@Axel Vogt Continuity of the integrand is obvious over a nonnegative interval. But Maple refuses to do the definite integral even when specific nonnegative limits of integration are supplied. Here's a simplified example:

Maple does the indefinite integral:

int(t^(5/2)*exp(t^(5/2)), t);

Maple does not do the definite integral:

int(t^(5/2)*exp(t^(5/2)), t= 0..1);

@dohashi I agree with itsme that Grid:-Map should just go ahead and copy the full memory state of the main kernel.

Regarding that fourth paragraph of ?Grid,Map : Could you show an example where Map is not "called outside the grid, on a non-compute node, not as part of a parallel computation"? I.e., an example where it is called as it "is designed to be called on all the nodes of a parallel computation." Would the last example at ?Grid,Map , where Launch is used, be an example of what the fourth paragraph refers to as "robust" use?

@itsme: Does your wrapped and Mapped procedure use global variables? I wonder if that is the issue. If it is, then maybe you could put the procedure in a module.

@mehdi jafari Numeric solution of PDEs without a time component is not currently implemented.

@Edinburgh Please explain what you are trying to plot. Is it simply the bezier curve with control points P[0]...P[5]? Why are you using complex numbers?