Hi Tomleslie here is the link to my code: NNMMM.mw

@dharr 

Thank you 

@tomleslie 

Thank you for the explanation

@tomleslie 

I'm kind of confused. Please in mydata, which column corresponds to xi =0 and xi =0.1. I want to differentiate them for plotting purposes.

@tomleslie

 Thank you very much. I really appreciate your help

@tomleslie 

Maybe I attached the wrong file. In the file: Case1_Revised_080922.mw the line has: 

ExcelTools:-Export(op([1, 1], plots:-display( convert( ans1[12..14,3], list))));

You may be wondering why am I doing this. While I wanted to export my data from maple to Matlab and compare it with my result in Matlab. I tried to export it as a Matlab script but couldn't work. So I had to use ExcelTools. So, I exported only my plots which are okay for only one value of xi. But for several values of xi, it couldn't work. If there are other alternatives, I need your suggestion.

@Carl Love

Thank you 

@mmcdara 

Unfortunately it didn't work on my maple version. I recieved this error message:

Error, invalid input: rand expects its 1st argument, r, to be of type {posint, integer .. integer}, but received 0. .. 1.
Error, invalid input: no implementation of Color matches the arguments in call, 'Color:-ModuleApply([r(), r(), r()])'

@dharr 

For instance how do I plot a function with singularities such that the asymptotes don't appear

f(x,y) = y/(x+1) + 1/(y^2-y) 

@dharr 

Unfortunely, plot3d don't accept discont=true as an argument.

@dharr 

Thanks. I actually use this as an instance. 

@Carl Love 

Thank you 

@acer 

Oh I see

@acer Okay thank you. I quite understand my project. I tried to do the same work in Matlab. I was having a different result which made me check my equations and I discovered there is a change of sign in the g equation. I used maple to check the analytical solutions and Matlab for numerical solutions. That was when I discovered the error. The only problem I have is the fact that when I change the sign from (-) to (+) in g equation maple only returns the real roots instead of the real and imaginary parts. When I solve it with Matlab, I discovered the function is purely imaginary and maple is not finding those roots. 

@acer 

While I'm passing judgment because the function r(q) is a complex function that is purely imaginary and Maple can't find all the roots. If there is any way I can buyout the process, I will appreciate it. 

Secondly, why is maple taking a long time to compute this: (gridji,rngj,rngi):=[81,81],0..3,-2.0..-0.0;

It took me more than 24 hours to do this (81 by 81 grids instead of 5 by 5) on my machine. Is there any way to shorten the time?