@Aakanksha 

It's simple:

plots[logplot]([result1, result2, result3]);

The reason that Markiyan's suggestion did not work is that your Maple is a very old version. What version is it?

@Aakanksha 

That should be plots[pointplot](result1), not plots[pointplot(result1)].

@quann169 

Please see the response by Alejandro Jakubi above and follow the link that he provided. That link leads to a post that has a video that explains how to deal with this situation. Unfortunately, you cannot use $include in Maple IDE, but apparently there is a viable alternative.

@k20057 5 The method that you describe for converting a Maxwell-Boltzmann sample to a Bose-Einstein sample does not make sense to me when you consider a large sample. For 9 quanta and 6 particles, the Maxwell-Boltzmann model has a population of 2002 different states and the Bose-Einstein model has 26 different states. If you have a sample size of 100,000, then every Maxwell-Boltzmann state is occupied about 50 times. If you convert that to a Bose-Einstein sample in the method that you describe, then all 26 states will be used with probability very close to 1. You might as well just look at the entire population of 26 states rather than a sample.

In other words, after making the conversion as you describe, there is no longer anything random about the sample.

What would make sense to me is taking the Maxwell-Boltzmann sample and reducing it by removing any members that are not in the Bose-Einstein population, and then further reducing it by removing any members that are not in the Fermi-Dirac population.

@H-R I was using Maple 18. The great increase in speed that I experienced was due to Maple 18's hardware float evaluation of Bessel functions.

Where is the N that you want to set to 100? I found an N in your code, and I set it to 100 with no problem. I also set k to 100 with no problem.

@quann169 I don't know anything about IDE, but did you remember to put the $ in the first column?

The code in your file is a statement defining a procedure. When that statement is executed by the read, it defines the procedure at the top level rather than as part of the module.

@derekcowley Oh, I think that I understand now what you want: You want to express a polynomial as a sum of terms each of which is factored a particular way. That's much more difficult than just factoring. I can't immediately think of a way to do that. Certainly the factor command cannot do it. It may be possible using simplify with side relations. The side relations would be p1 - p2 = p12, p1 - p3 = p13, etc., (every possible pair of p's is equated to a unique variable). I may be able to work on it tonight.

@derekcowley 

You seem to expect that your expression T2 ((24) or (25)) will be divisible by p1-p2. I am saying that it is obviously not divisible by that.

@Markiyan Hirnyk 

What are you talking about??? The only thing that I did to make the OP's code work was to change the fsolve command as indicated. The answer is complete and true.

I get sqrt(14) for both norms. What do you get?

@Markiyan Hirnyk 

Here is my executed worksheet. I get similar results with Maple 18.00, with Digits = 15, and/or with alfa set to 1/5 rather than 0.2.

fsolve_problem.mw

It is natural to think that specfunc(name, `+`) would be the correct typespec. However, x+y is not considered a function by Maple (so specfunc doesn't apply), even though it is the result of applying the function `+` to arguments of type name. It's kind of like sin(0) is not of type specfunc(anything, sin) because sin(0) is not a function after it's evaluated (it's just 0).

`+` and `*` are their own fundamental types.

@Preben Alsholm 

This can also be written

map(`/`, c, d);

which avoids two calls to user arrow procedures.

@Carl Love I still don't see your attached file.

Update: Okay, I got the attached file now.