It is possible the plots created do not match the numeric values they state. My results match yours and do not match theirs even for the gamma=0.1 case, I only get the gamma=0.1 behaviour they show if my value is around 0.2. 

I got no errors, answer returned

Dsolve_Reply.mw

@ecterrab 

Thank you as well. 

@C_R 

I assumed that it was something like this, good to know my suspicions were correct. 

You are using multiplications between your partials and your potentials, you need to use d_[mu](A) if you want to have the partial actually evaluate on A. Secondly,  in other spots you are not multiplying two terms together just writing a single term (example gr, should be g*r). 

I think I fixed your problems in the attached code. Also, shouldnt you have some of your indices in a one up one down configuation in F? 

SU(2)_Reply.mw

@segfault 

Good luck with whatever you are trying to accomplish.

@segfault 

I was going to post an answer on here for you before I saw this exchange. 

Considering @ecterrab is the head of Maple Physics I would suggest you be much more respectful to his replies in the future. I have had issue with Maple physics in the past where I thought there were bugs but he appropriately responded and showed that some of the issues were my fault and I was not computing things properly. 

Regarding "pre-publication work", I have not once had my work stolen and published before I was done simply by asking for help on this Forum. This is not the kind of place people go scoping for ideas for manuscripts, it is for questions about Maple. 

If you look at my attached sheet you can do what you want.. if I am reading your issues properly. 

Define_MP.mw

@Carl Love 

Will use that in the future. 

@vv 

Not sure how I forgot to mention numeric solutions as I knew that it worked. Thanks for adding that. 

It is all dependent on the system of differential equations you have. Some systems cannot be solved by "analytic" means and must be solved numericlly. 

1) If by connection you are referring to Christoffel symbols then you must remember that they are not tensors, so the error is referring to that portion I believe. 

2) Could you upload your worksheet with the green arrow or tell me your H tensor is?

3) I suggest you use the Physics package for this problem, the commands are much simpler to use in my opinion. See attached where I create an expression with a single index contraction of Christoffel and Ricci tensor for a arbitrary static spherically symmetric metric. 

Contraction.mw

@acer 

Good points, I suppose since I knew how the final answer should look I used symbolic since it would take care of the a terms without needing the positive assumption. 

It is quite dodgy sometimes - well most of the time for complex expressions - but the simplicity of this simplification is what made me decide to use it. 

Worksheet from when I updated to Maple 2024 and ran with the updated Physics Package. 

Can you upload your worksheet? 

Or at the very least post what your ODEs are and the initial conditions? 

@acer 

I would say that I would prefer accuracy over timing here, although a fast result is something desirable to see the rough result. Ofcourse the end goal would always be a quick and accurate code. 

I appreciate this other approach as it could provide some useful insights, although I might need to spend some more time looking at it to make sure i am understanding it correctly. I do enjoy the "smoothness" you mention regarding the fact that my N list is small and creates a course pointplot result. 

I will definetely play around with this and see what it does. 

Thank you.