Download 219_fsolveproblem.mpl
View file details

Are you defining Res as I did?

Did you execute all lines previous to the ones mentioned last?

Preben Alsholm

Download 219_fsolveproblem.mpl
View file details

Are you defining Res as I did?

Did you execute all lines previous to the ones mentioned last?

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 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