@vv Thanks, I was mistaken, but now I understand due to your example. 

@ecterrab The member @nm has a long-standing habit of deleting their own Questions (that haven't been Answered), Answers (once another Answer appears), and Comments. I've tried to discourage them from doing that, but I don't think they had any response. In my many years of being a Moderator here, I don't recall ever deleting anything by @nm nor seeing anything that I thought should be deleted. Indeed, I consider @nm to be the all-time most-prolific poster of good-quality Questions, and a substantial proportion of their 6000+ reputation points come from me voting up their Questions.

I suspect that the vast majority of cases where a Reader notices that a worthwhile Question has been deleted fall into one of these categories:

1. The author deleted their own Question.
2. The Question is a duplicate of (or substantially the same as) another recent Question that the Reader didn't notice.
3. The Question is a followup to an Answer to a recent Question by the same author, and the Reader hasn't noticed that (yet).

Note that there is no mechanism whatsoever on MaplePrimes for Moderators to do any of the following things; their options are only Edit, Reply, or Delete:

1. To move a Question into another thread (such as to the thread for which it is a followup) where it rightfully belongs;
2. To move material to a holding area for possible further review.

I'd be happy to use these things (or any other StackExchange-like features) if MaplePrimes implemented them.

@Alex0099 I suspect that the remaining errors are deeply dependent on the underlying mathematics, which I know nothing about unfortunately. The semicolon error that I corrected was just a simple syntax error.

@vv You wrote:

• For an abelian group, ... if the order is even, the product (of the elements) may be <>1. 

I think that your "may be" can be strengthened to "is". Isn't the product necessarily equal to the unique nonidentity element that is its own inverse (such as -1)?

@Alex0099 To correct the "Error, unterminated loop", you need a semicolon at the end of the line tmp:= RHS(x)*pn(i, 1, t).

@ijuptilk What I gave you is only the code for the procedure definition of simple_deflate. To get results, you need to apply it to arguments, as in simple_deflate(A[6], u). To do that, you need the vector u. The top of the first page of your PDF excerpt shows that u is constructed from w. How to get w is presumably on an earlier page of the PDF. If you show the part of the PDF that explains how to get w, I can probably condense the entire algorithm to 2 or 3 lines of code.

@ijuptilk There is some invisible invalid character(s) in the whitespace at beginning of the 2nd line of simple_deflate. If I backspace over the whitespace (including the line break) and then reinsert it, then it works.

@ijuptilk Please post a worksheet showing the error as an attached file. Use the green up arrow on the toolbar of the MaplePrimes editor.

@ijuptilk Your problem is likely caused by the lines "Approximate eigenvalues" and "Precise eigenvalues". These lines are comments and you're treating them as code. To keep these lines, use # to tell the parser that they are comments:

#Approximate eigenvalues
#Precise eigenvalues

@ijuptilk The procedures pindex and simple_deflat from the PDF can both be replaced by this simple procedure:

simple_deflate:= (A::Matrix, u::Vector)->
    (p-> A-u.A[p]/u[p])(max[index](abs~(u)))
:

This code also addresses the "ill-conditioning" issue mentioned in the last line of the PDF by selecting the entry of maximal absolute value rather than any nonzero entry.

If your Maple version is too old to have max with the index option, let me know, and I can modify the above. You'd be able to tell immediately from the help page ?max.

@Earl Yes, Kitonum's Answer assumes a uniform density of 1. Regardless of whether the density is uniform, it can be accounted for by including it as a factor of the integrand in all of the integrals (2 for the center of mass and 1 for the moment of inertia). If that factor is constant (uniform density), it can be pulled out of the integrals. The dimensions of moment of inertia are mass*length^2.

@mnovaes The command subs has much more general uses than substituting values for variables. It can substitute any expression for any other expression provided that the expression being searched for exists as a distinct entity within the expression being modified. I fairly often substitute one number for another. About half of the time that I do this, the number that I want to change is represented by a variable. If subs treated x as a plain name in the following code, then I wouldn't be able to do that number changing:

x:= 3:
subs(x= 4, 3);

Maple can generate a "series" solution for your ODE, meaning the exact values of any finite number of leading terms of something akin to a power series for b(x). Specifically, it returns _C1*P1(x) + _C2*P2(x) + _C3*ln(x)*P3(x) + P4(x) where the _Ck are constants of integration and the Pk are power series with rational coefficients truncated to some finite order that you can specify.

Order:= 20: #desired truncation order
dsolve(
    diff(b(x),x$3)/2 + diff(b(x),x$2)/x + 
        2*diff(b(x),x)/(x^2+1)^2 - 4*x*b(x)/(x^2+1)^3, 
    b(x), series,
    x= 0 #expansion point (any complex number or `infinity`)
);
Order:= 6: #reset to default value

There are other cases where it's possible to get a "formal series" solution, meaning explicitly expressed as infinite sums, but I can't do it with this ODE, despite the fact that the reformulated ODE is homogenous and linear with integer polynomial coefficients.

@Angie7 In the 2nd of my 3 answers, change x(1) to D(x)(1).

I deleted your earlier version of this Question. I see many errors in your code, but I'm on my phone (no Maple) right now. I'll Answer soon if someone else doesn't get to it first. The command fsolve is very good with numeric polynomials.