@bfathi I cannot find a freely available copy of that file on the web. It is part of that book. Is there a CD attached to the book? Does the text of the file appear in the book? You're in luck that I happen to have a copy of that book, the third edition (Walter Gander and Jiri Hrebicek, Solving Problems in Scientific Computing Using Maple and Matlab, 3rd ed., Springer 1997). If you tell me the section of the book that the file appears in, I may be able to help you more.

Your English is okay. I had no trouble understanding any of what you were saying.

@jahan The shoot package is not part of standard Maple. You'll need to download it somewhere. I have seen it before. It is very old, and may be difficult to use in today's Maple. Note that Maple does have BVP solvers built in. 

@bfathi 

Use the fully qualified file name, for example,

read `C:/Users/Carl/Desktop/dsnumsort.m`;

@14Robin I don't think that that's possible from the keyboard alone. And for me the following is definitely too much effort to be worth it. But if you really want to get rid of the [, go to the Format menu and select Create Document Block. Then press uparrow, control-Del, control-T.

@Markiyan Hirnyk Thank you. The misspelling "Newoton" in the original title probably pushed this further down the list of search results than it otherwise would've been.

@Markiyan Hirnyk 

The OP's only issue is with the prettyprinted display of %limit, not with its programmatic functionality.

This does not go back to Maple Vr4. For a discussion of forming inert functions by prepending %, see ?value.

@Markiyan Hirnyk The OEIS claims that 34 is the correct answer for order 5 graphs. It is amazing that Categorize and IsIsomorphic get this right when they are so wrong for order 4 graphs.

@Markiyan Hirnyk 

The MaplePrimes search is horrible. And Google does not properly index MaplePrimes: All searches point to the User's Answers page, which is worthless in my case because I have over 1800 Answers. So, I can't find my own Answer. So could you please post a link to it?

Note that I can easily use Google to find my Answers on StackOverflow, Yahoo groups, comp.soft-sys.math.maple, etc. This situation that MaplePrimes is not properly indexed by Google is deplorable. What a waste of valuble information!

@Andriy Yes, I believe it was added in Maple 18.

@Markiyan Hirnyk 

Yes, of course: Your x-values are not evenly spaced. Further hint: Represent that spacing by h.

@KimonoMyHouse After you Replied, I added information about your second error to my Answer above.

@Preben Alsholm

Here is the version with limit:

power:= (E::algebraic, t::name)->
     limit(simplify(expand(t*diff(expand(ln(E)), t))), t= infinity)
:

When used with symbolic exponents, this reveals some interesting limitations/bugs of limit:

power(t^a+t^2, t) assuming a>2;

                 a

That's fine. But...

power(t^a+t^2, t) assuming a<2;

                 limit((t^a*a+2*t^2)/(t^a+t^2), t = infinity)

And...

power(t^a+t^2, t) assuming a<2, a>0;

                 a

That's wrong. It should be 2.

@Preben Alsholm 

It works when the power of t is the same in all terms in the sum (unless you have a counterexample). The OP asked for the power of t, not the highest power of t. If the sum is of terms with different powers of t, then power returns an expression that depends on t that may have some significance that I haven't figured out yet. Perhaps I should take the limit as t -> infinity of this expression. That would give the highest power of t.

@Kitonum 

Your two procedures produce the same results, but the second is far more efficient. It can be made even more efficient and easier to read by using the structured type identical(t)^integer instead of `^` (see ?type,structure) and by using map2 to avoid calls to your anonymous arrow procedure s-> op(s)[2]. Like this:

power:= (E::algebraic, t::name)-> max(map2(op, 2, indets(E, identical(t)^integer))):

@Kitonum 

You need to at least specify the variable with respect to which the degree is being measured.