@Preben Alsholm 

Your bcsy has boundary conditions specified at x=0, x=1, and x=inf.

You asked the same question five days ago, here, and I answered it. Did you have some problem using that answer?

@nm There is no way in Maple to change an operator's precedence (see ?precedence ). &= has a rather high precedence, and := has the absolute lowest. I can change it so that the grouping symbol is unevaluation quotes:

x &= 'r1 + r2';

or any other grouping symbol, but there must be some grouping symbol.

@Ratch And where did you put the sort command?

@Markiyan Hirnyk 

I agree that int and Int should give infinity and Float(infinity) as answers.

While the overall series is not clear, it is clear that the coefficient of the (x-Pi)^(-1) term is not zero and has modulus 1.

@Markiyan Hirnyk Please clarify your statement about the general solution. And what do you mean by "exception"? It seems that the general solution cannot be produced by dsolve.

@BobWild Do you mean that you want to solve part b, i.e., get the analytical solution? That is just the same answer that I gave before with different specific values.

Or do you mean that you want to solve part c, i.e., get the matrix/vector solution using a stepsize of h=Pi/8? That's more complicated. To help you with that, I would need to see the page that precedes equation 11.18 in your book. Specifically, I need to see how I the mesh points are defined. I don't know if they are using a midpoint or endpoint method. You will need to set up the matrix and vector for that, because the mesh options to dsolve(..., numeric) will not let you make a mesh that small.

@BobWild Correct: You do not need to define the mesh points. That is done automatically. Often it helps to adjust the number of mesh points. This is done with the options initmesh and maxmesh. See ?dsolve,bvp for precise details, or post a specific BVP problem here. But it is best to not add any options unless you initially get an error message without them. An error message usually gives a suggestion for which option to adjust and which direction to adjust it.

There is no need to create a matrix and no option for it either. That is all handled internally. You can pass in a vector for an approximate solution. And you can receive the output in matrix/vector form, but the procedure form of output is usually more convenient, with interpolations handled automatically.

Please read the help page ?dsolve,bvp .

@BobWild That is correct; only the boundary conditions are needed.

@Ratch I see what you mean. The first two terms should be 1+x instead of x+1, and when the slider is set to 49 or 50, the last two terms are inverted also.

I can't find the underlying startup code. How do I get to it?

Could you please provide more details of your situation? some example values perhaps? or a worksheet?

@Carl Love I just significantly updated this Answer to use plots:-loglogplot rather than a simple point plot with logarithmic axes. So, it is more in line with your question.

@highsciguy After applying the first evalindets, hit it with

evalindets(%, identical(1.), 1);

Please, let me know how that goes!

@Joe Riel For me, your post that compares all the Cartesian product procedures is the top hit when I enter "Cartesian product" (without quotes) into the MaplePrimes search bar. It is dated April 27, 2010, which may be why you're having trouble finding it. There are hundreds (maybe thousands) of posts and questions with that date, so I guess that a lot of old posts had their date changed to that in some "update".

@Markiyan Hirnyk Is that entrance exam for all students entering the University, or just for those majoring in math?