@mmcdara If you check the subscribe box when you post/answer/reply, then you will get email notifications and the thread will remain alive. But perhaps you won't notice this and find out.  :)

In terms of the possibilty of spacecurve coloring, when I found spacecurve gave me an error I assumed it couldn't be done. I was too lazy to check the plots,structure help page, but I used the word "think" not to be definitive.

So now I "think" only implicitplot3d's ISOSURFACEs can't take a color function. That one has frustrated me for a long time since I want to add phase colouring to those surfaces.

@mehdibgh Your Eig method doesn't work because the pseudoinverse is singular. 

I attach my version of your LV method, which avoids taking an inverse and makes some use of the fact that M is diagonal, though that could be avoided. I only did operations that kept the matrices symmetric, and avoided some of the slow conversions you had. Accuracy probably cannot be improved, because some data in M11 are small and inaccurate, so the inverse square root now has large inaccurate values. As you tried, some scaling can help but the matrices are all ill-conditioned. I'm guessing the input data does not have enough sig. figs.

Problem2.mw

@ecterrab Actually, I stole my ranges from DLMF. I debated putting a link to the nice rotatable graphics there, but I am trying to promote Maple!

@janhardo Your plot of the integral is just a horizontal line because the evaluated integral is just the number Pi*cosh(1)-Pi*sinh(1), and the default plot range is -10..10. So it is like plot(1.8,x=-10..10).

cos(z) is not exp(I*z). Try convert(cos(z),exp) to see what it is.

Suggest you use complexplot to visualize some of the functions you are using.

@acer @Carl Love Thanks for the info. I'm a little surprised that the default would become extended but then the special function rules would be off by default - it seems many probably share the OPs confusion.  

@Kevin Dragnet  I think you should contact technical support - something definitely wrong here. https://www.maplesoft.com/support/ [Edit: see @Carl Love's answer below]

MeijerG typesets in mathematical notation; I don't have FoxH in my version.

Download typesetting.mw

@mmcdara I'd probably combine these two, so 1/(sqrt(1-x)) = (1-x)^(-1/2) =1-(-1/2)x+... approx= 1+(1/2)x for x small, using the binomial expansion, (1+x)^n = 1+n*x+... or (1-x)^n = 1-x*n

@The function See the help page ?pochhammer, which says

pochhammer(z, n) = z*(z+1) .. z+n-1

so there are n factors going up in steps of 1, starting at z.

@Carl Love convert(sqrt(x), FormalPowerSeries, x = 1) will work directly here.

@Katatonia See edit to my answer below.

@Katatonia Perhaps my version of Maple is too old, but I only see the definition of the piecewise function, and not the dsolve command. But you seem still to have p*T=0, which as an initial condition should be p(0)=0.

@nm I was not suggesting that the equivalent interface options would work - I tried that myself and found they did not. My point is that since these plot options appear under "interface" that they may be part of the "communication with the user interface" (see the interface help page), which is turned off when the worksheet is run "headless". From the RunWorksheet page: "The invoked worksheet runs "headless", meaning that it will not appear with a user interface." (my italics).

@FDS To convert Pa to mmHg use convert(27., 'units', 'Pa', 'mmHg'). You can just use x=0..Bloodpath as the odeplot range, with useunits=true.

To find the value at Bloodpath use sol(Bloodpath). I guess in this case you need to remove units.

Please upload your Maple worksheet using the green up-arrow option in Mapleprimes.

I don't know much about this, but in Maple 2015 FunctionAdvisor("branch_cuts",HeunC); returns "The location of the "branch cuts" for HeunC is unknown to the FunctionAdvisor"; compare FunctionAdvisor("branch_cuts",BesselJ); which returns "[BesselJ(a, z), And(a::(Not(integer)), z < 0)]". Likewise, MathematicalFunctions:-Get("branch_cuts",HeunC); returns empty. 

There is further information in other parts of the MathematicalFunctions package such as Evalf that may be useful.