You may try setting optimalitytolerance:

res:=LSSolve(Res,{x1<=1,x2<=1,x3<=1,x4<=1,x5<=1,x6<=1,x7<=1,x8<=1,x9<=1,x10<=1,T>=Ti},assume=nonnegative,iterationlimit=5000,optimalitytolerance=0.0001);
[0.00517923021408238112, [T = 2782.09143970567812, n = 10.1585637591277820,

  x1 = 0.00872658980230692216, x10 = 0.0748602181855385214,

  x2 = 0.0766666505398920090, x3 = 0.579521666596851892,

  x4 = 0.0000115210794556045040, x5 = 0.0901704459988488805,

  x6 = 0.0852556575745576917, x7 = 0.00577381343641316520,

  x8 = 0.00427475756541988076, x9 = 0.0642723480725257840]]

See the help page: ?Optimization, Options

Preben Alsholm

You may try setting optimalitytolerance:

res:=LSSolve(Res,{x1<=1,x2<=1,x3<=1,x4<=1,x5<=1,x6<=1,x7<=1,x8<=1,x9<=1,x10<=1,T>=Ti},assume=nonnegative,iterationlimit=5000,optimalitytolerance=0.0001);
[0.00517923021408238112, [T = 2782.09143970567812, n = 10.1585637591277820,

  x1 = 0.00872658980230692216, x10 = 0.0748602181855385214,

  x2 = 0.0766666505398920090, x3 = 0.579521666596851892,

  x4 = 0.0000115210794556045040, x5 = 0.0901704459988488805,

  x6 = 0.0852556575745576917, x7 = 0.00577381343641316520,

  x8 = 0.00427475756541988076, x9 = 0.0642723480725257840]]

See the help page: ?Optimization, Options

Preben Alsholm

You may also want to scale the equations.

Here I have scaled the last one:

Res:=map(lhs-rhs,[eq1, eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,eq10,eq11,eq12*10^(-2)]):
res:=LSSolve(Res,{x1<=1,x2<=1,x3<=1,x4<=1,x5<=1,x6<=1,x7<=1,x8<=1,x9<=1,x10<=1,T>=Ti},assume=nonnegative,iterationlimit=5000);
Warning, limiting number of major iterations has been reached
[0.00517609062184478092, [T = 2782.44660329790214, n = 10.1582899935774478,

  x1 = 0.00873618614514319950, x10 = 0.0748388854712784686,

  x2 = 0.0767030979399348362, x3 = 0.579548942356043040,

  x4 = 0.0000115713864132603826, x5 = 0.0901633563755876583,

  x6 = 0.0852311188569541273, x7 = 0.00577787937886131958,

  x8 = 0.00427369013816442536, x9 = 0.0642550496232529162]]

Preben Alsholm
 

You may also want to scale the equations.

Here I have scaled the last one:

Res:=map(lhs-rhs,[eq1, eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,eq10,eq11,eq12*10^(-2)]):
res:=LSSolve(Res,{x1<=1,x2<=1,x3<=1,x4<=1,x5<=1,x6<=1,x7<=1,x8<=1,x9<=1,x10<=1,T>=Ti},assume=nonnegative,iterationlimit=5000);
Warning, limiting number of major iterations has been reached
[0.00517609062184478092, [T = 2782.44660329790214, n = 10.1582899935774478,

  x1 = 0.00873618614514319950, x10 = 0.0748388854712784686,

  x2 = 0.0767030979399348362, x3 = 0.579548942356043040,

  x4 = 0.0000115713864132603826, x5 = 0.0901633563755876583,

  x6 = 0.0852311188569541273, x7 = 0.00577787937886131958,

  x8 = 0.00427369013816442536, x9 = 0.0642550496232529162]]

Preben Alsholm
 

This is interesting.

I think op(RS1S) should be op([1,1..-1],RS1S).

Preben Alsholm

This is interesting.

I think op(RS1S) should be op([1,1..-1],RS1S).

Preben Alsholm

I'm using Firefox myself, and mine just asks me if I want to open mulleralg.mw in Maple 13, to which I just hit OK.

Then the worksheet is opened in Maple 13.

I do have Internet Explorer on my computer. If I just execute the link to the worksheet, the screen is filled with xml-code in which the mw-worksheet is written.

So another suggestion is:

Go to the shorter link

http://www.cbu.edu/~wschrein/media/ANA/

which is a directory of a lot of files.

Go to mulleralg.mw.

Right-click on it. Select "Save destination as" (or something very similar; mine speaks Danish).
Then edit the file name so it becomes mulleralg.mw (not the suggested mulleralg.xml) and change file type to "All files".

Now save the file in a place, where you can find it.

Finally open it in Maple and continue as before.

Preben Alsholm

I'm using Firefox myself, and mine just asks me if I want to open mulleralg.mw in Maple 13, to which I just hit OK.

Then the worksheet is opened in Maple 13.

I do have Internet Explorer on my computer. If I just execute the link to the worksheet, the screen is filled with xml-code in which the mw-worksheet is written.

So another suggestion is:

Go to the shorter link

http://www.cbu.edu/~wschrein/media/ANA/

which is a directory of a lot of files.

Go to mulleralg.mw.

Right-click on it. Select "Save destination as" (or something very similar; mine speaks Danish).
Then edit the file name so it becomes mulleralg.mw (not the suggested mulleralg.xml) and change file type to "All files".

Now save the file in a place, where you can find it.

Finally open it in Maple and continue as before.

Preben Alsholm

The link should be to mulleralg.mw, i.e. this one

http://www.cbu.edu/~wschrein/media/ANA/mulleralg.mw

Now after executing the procedure (or in this case the whole worksheet) define p do muller(p, -13, -12, -11, 0.1e-5, 100, r1);

Preben Alsholm

The link should be to mulleralg.mw, i.e. this one

http://www.cbu.edu/~wschrein/media/ANA/mulleralg.mw

Now after executing the procedure (or in this case the whole worksheet) define p do muller(p, -13, -12, -11, 0.1e-5, 100, r1);

Preben Alsholm

You could let infinity be 5 (or whatever you like):

f := piecewise(t < 0, 0,t < 1, 1, t < 2, (-2+t)^2, t=4,5,0);
plot(f,t=-10..10,sample=[$-10..10]);
 

Preben Alsholm

You could let infinity be 5 (or whatever you like):

f := piecewise(t < 0, 0,t < 1, 1, t < 2, (-2+t)^2, t=4,5,0);
plot(f,t=-10..10,sample=[$-10..10]);
 

Preben Alsholm

You didn't execute the procedure muller.

In the worksheet mulleralg.mw execute the muller procedure, i.e. put the cursor on it and press Enter.
Then define p by etc. ....

If you executed the muller procedure and then a restart, then the muller procedure is not known to Maple.

Preben Alsholm

You didn't execute the procedure muller.

In the worksheet mulleralg.mw execute the muller procedure, i.e. put the cursor on it and press Enter.
Then define p by etc. ....

If you executed the muller procedure and then a restart, then the muller procedure is not known to Maple.

Preben Alsholm

The following should work, and works for me:

In the worksheet mulleralg.mw execute the muller procedure, i.e. put the cursor on it and press Enter.
Define p by

p:=x^5+11*x^4-21*x^3-10*x^2-21*x-5;
then do

muller(p,-13,-12,-11,.000001,100,r1);

The result is then the output and is printed to the screen as usual. The procedure is made to assign the result to the last argument (r1), which therefore should be a name. It is best to replace r1 by 'r1' as in

muller(p,-13,-12,-11,.000001,100, 'r1');

Preben Alsholm