I am looking at the d'Almbert operator for the Minkowski metric in cartesian coordinates. When using the dAlmbertian command and converting to partial derivatives it gets the coefficients wrong, the negatives are not in the correct place based on the metric. But when I manually write the operator with the metric and derivative operator it returns the correct answer. Any thoughts on this? 

See attached, thank you. 

dAlmbert.mw

Suppose I have a metric g, and I want to perform a conformal transformation g = exp(2Phi(X))*g, is there a straightforward way to do this for curvature quantitieies (Christoffel, Ricci Scalar etc)? I was able to do it rather easily for the Christoffel symbols, as seen below, but it required me making a substitution for each index pair. While this isn't horrible, it would be nice if there was a way to do it without this procedure.

** Edited to make it Phi(X) 

Any thoughts appreciated, thank you!

MyConformal.mw

I have a procedure that I am trying to run that would be an improvement/more sophisticated way of solving a problem that I have previously solved. When I try and run my procedure I am getting an error, and from what I gather with the error is that there are some values when inserted into my procedure that cannot be evaluated. Just for context it is a procedure that contains numerical solutions to a system of DEs and and contains inequalities. 

I would like to know is there an easy method to figure out what values are giving me this error? 

Or a follow up, is there something wrong with my procedure that is giving me this error? I have included some commentary in my workshet as well to hopefully make everything clear. 

Thanks. 

Proc_Error.mw

I have a differential equation which I am looking for a series solution in inverse powers of r. I am doing this by matching inverse powers between the RHS and LHS of the differential equation and then finding the coefficient in the series solution that would solve this. I have written a little procedure which returns the coefficients, I am just curious if there is a command that will do this for me or if i can improve on this technique? I tried the Solve command in PDETools but it just kept evaluating and did not return an answer. 

Thanks in advance. 

EDIT: More detail provided in the first reply. 

coefficient_question.mw

I have a thirder order ODE with non polynomial coefficients and I naively thought to try dsolve for fun to see what happens and Maple returned DESol with a second order differential equation and an arbitrary coefficient. I know Maple outputs DESol when it cannot find a solution similar to RootOf but the arbitrary constant is what is throwing me off. 

I am unsure how to interpret this, if a particular solution is found I could reduce the order and see how I could get with the second order ODE but maple doesn't produce a particular solution when I run that command. 

DESol_Question.mw