I was curious to know if one can extract a specific solution from a LPSolve routine. 

As an example, consider the following output to a constrained linear problem. The objective value is 8 and the decision variable values (binary) are given.

Sol := [8, [w[1, 1] = 1., x[0, 0, 1] = 0, x[0, 1, 1] = 1, x[0, 2, 1] = 0, x[1, 0, 1] = 0, x[1, 1, 1] = 0, x[1, 2, 1] = 1, x[2, 0, 1] = 0, x[2, 1, 1] = 0, x[2, 2, 1] = 0, y[0, 0] = 0., y[0, 1] = 0., y[1, 1] = 2.]]

I am interested to know if we can isolate any variable value from this solution. I know that Sol[1] will return 8, and Sol[2] will return the remaining terms. But what if I wanted, say, x[1,2,1] alone?

Thanks for reading!

A bit of an annoyance. 

typing 'numerator' / 'denominator' generates an error.  The first quote never gets automatically grouped as it should.

Two workarounds.  The first is to move to, and delete the first quote and re-enter it again in front of numerator.
The second is to use brackets, although one shouldn't have to.

Let us consider 

restart; 
MultiSeries:-limit(sin(n)/n, n = infinity, complex);
0

The answer is wrong: in view of the Casorati-Weierstrass theorem the limit does not exist. Let us try another limit command of Maple

limit(sin(n)/n, n = infinity, complex);


(lim) (sin(n))/(n)

which fails. Therefore, Maple user does not obtain the correct answer. 

Hi

I have this problem that my computer crashed during Maple 2016.1 to 2016.2 update.

Now it won't start at all.

I have tried reinstalling Maple several times, tried not updating to 2016.1 and running the program as vesion 2016.0.

I get this error in my eventlog, when I try to run the program.

Name of program errors: javaw.exe, version: 8.0.1020.14, timestamp: 0x576af3f4
Name of the module with error: nvSCPAPI64.dll, version: 7.17.13.7849, timestamp: 0x588218a5
Exception code: 0xc0000409
Displacement by mistake 0x0000000000034b2f
Process ID 0x166c
The program's start time 0x01d279add4515202
Program path: E: \ School Program Files \ Maple 2016 \ jre \ bin \ javaw.exe
Module path: C: \ Program Files (x86) \ NVIDIA Corporation \ 3D Vision \ nvSCPAPI64.dll
Report ID: d52f952d-8fd2-4130-890a-196e6a7bdc92
Full name of the program with the error:
Relatively application ID for program errors:

So as my titel suggests I need help to remove all files associated with Maple so I can install the program anew. It has worked fine for over a year so far. Only stopped working after that crash.

I have even tried installing the older version, Maple 2015 and it doesn't work either, so I am guessing the files which are corrupted aren't removed during a normal uninstall process, and they are shared among the versions of the program.

Thanks in advance :)

Hi everybody, 

The sequence 
    writeto(MyFile)
    showstat(Myproc)
keeps printing the result on the screen

How is it possible to redirected the output of showstat(Myproc) into a file ?

Thanks in advance

Hello,

I've just started to use Maple 16, and I can't seem to get a neat result, when I use the solve function to solve 4 equations with 4 unknown constants.

I've posted my calculations.

How can I get a simpler result? I found a simplification command, but that didn't help

Where I made a mistake I got this error in for loop and how to fix it?

Download quest.mw

I have a need to calculate base 2 math as would be done in an integrated circuit. Math will be done using 15 bit 2's complement mantissa and 8 bit exponent. (the exponent is always assumed to be negative) We need to perform multiplication and addition, where each internal operation is represented by such a number.

The goal is to best represent the errors associated with such calculations and export the resulting code to a hardware description language for implementation on an integrated circuit.

Is there a package in Maple that can do that? Any advice on how te proceed?

How to do for to built animation Cylinder inscribed in cone?( in optimization problem)

I'm trying to plot the velocity of a ball thrown upwards with air resistance proportional to v^2 and also some simpler forms of this.

But the solution to v^2 returns root of and the plot stops for some specific time value. How can I proceed this plot to let's say 10 sec?

Staffan

Download tal_3.9_sid_66_b.mw

Dear all,

Is there any way to fill the upper region in plot3d?

e.g. if I put the filled option in 

plot3d(y^3+x^2, x = 0 .. 1, y = -1 .. 1, filled)

then the region between the x-y plane and my plot will be filled. What if I want to fill the upper region?

The reason is that I want to show the upper region is acceptable for me, but I couldn't find any other way. 

If you have any solutions for me, I appreciate it.

Hi,

I am a maths teacher and would like to create an ebook for my students with explanations, examples and quizes in the same time. They should be able to access and work through it, even on mobile devices, using the free Maple player. I am new in Maple and don't know if all this is possible. I would like to see some samples of code to learn from and of course, some advice. I would highly apprecaite your support. Thank you.

One of the things I like about Maple, is that I am able to look at most of the source code, and learn from it. But since I am newbie in Maple, sometimes I find it hard to know which next call I need to showstat() in order to see what happens next and follow the code.

Sometimes I do not know how to find the proc() being called in the listing shown by showstat().

For an example, I wanted to see how Maple implements AiryAi(x).

So I did  showstat(`AiryAi`) and I see the listing of AiryAi. The top level of the function. But then I see it calling, which must be some other level of `AiryAi` function? For example, at line 16, it says

     res := ('AiryAi')(x)

Since I am looking at AiryAi itself, this must be either recursive call or calling other internal AiryAi?

But if it recursive, then I do not see where the actual implementation is? I wanted to see if it uses the integral form of AiryAi or series form (not likely, uses GAMMA). I see the same recursive call in many places, such as line 35 and line 37 and othere places.

But looking at all the 42 lines, I do not see where the actual calculation of AiryAi(x) is done. Did I miss it?

There is a lot of error checking and looking at special cases and such. But how does one know from looking at this listing, which other function they need to showstat() in order to see the actual implementation? There must be more to it than those 42 lines?

This is on Maple 2016.2

Hi all,

I am using Maple 2016.

I have defined 5 polynomials: f1, f2, f3, f4 and f5 with 5 unknowns q1,q2 ,q3, q4 and lamda.

After this, I generated the Gröbner basis. But when I try to find the normal set I got an error.

with(Groebner);

f1 := lamda*q1-(3380075947548081*q1*(1/140737488355328)-259050600068343*q2*(1/140737488355328)-1826834460600733*q3*(1/1125899906842624)+4414049272733425*q4*(1/9007199254740992))*(q2*(8289619202186977*q1*(1/9007199254740992)+3380075947548081*q2*(1/281474976710656)-4414049272733425*q3*(1/18014398509481984)-1826834460600733*q4*(1/2251799813685248))+q3*(1826834460600733*q1*(1/2251799813685248)-4414049272733425*q2*(1/18014398509481984)+843667886835955*q3*(1/70368744177664)-215663898201129*q4*(1/9007199254740992))-q4*(4414049272733425*q1*(1/18014398509481984)+1826834460600733*q2*(1/2251799813685248)+431327796402257*q3*(1/18014398509481984)+843667886835955*q4*(1/70368744177664))-q1*(3380075947548081*q1*(1/281474976710656)-259050600068343*q2*(1/281474976710656)-1826834460600733*q3*(1/2251799813685248)+4414049272733425*q4*(1/18014398509481984)));
f2 := lamda*q2+(259050600068343*q1*(1/140737488355328)+3380075947548081*q2*(1/140737488355328)-4414049272733425*q3*(1/9007199254740992)-1826834460600733*q4*(1/1125899906842624))*(q2*(8289619202186977*q1*(1/9007199254740992)+3380075947548081*q2*(1/281474976710656)-4414049272733425*q3*(1/18014398509481984)-1826834460600733*q4*(1/2251799813685248))+q3*(1826834460600733*q1*(1/2251799813685248)-4414049272733425*q2*(1/18014398509481984)+843667886835955*q3*(1/70368744177664)-215663898201129*q4*(1/9007199254740992))-q4*(4414049272733425*q1*(1/18014398509481984)+1826834460600733*q2*(1/2251799813685248)+431327796402257*q3*(1/18014398509481984)+843667886835955*q4*(1/70368744177664))-q1*(3380075947548081*q1*(1/281474976710656)-259050600068343*q2*(1/281474976710656)-1826834460600733*q3*(1/2251799813685248)+4414049272733425*q4*(1/18014398509481984)));
f3 := (1826834460600733*q1*(1/1125899906842624)-4414049272733425*q2*(1/9007199254740992)+843667886835955*q3*(1/35184372088832)-862655592804515*q4*(1/18014398509481984))*(q2*(8289619202186977*q1*(1/9007199254740992)+3380075947548081*q2*(1/281474976710656)-4414049272733425*q3*(1/18014398509481984)-1826834460600733*q4*(1/2251799813685248))+q3*(1826834460600733*q1*(1/2251799813685248)-4414049272733425*q2*(1/18014398509481984)+843667886835955*q3*(1/70368744177664)-215663898201129*q4*(1/9007199254740992))-q4*(4414049272733425*q1*(1/18014398509481984)+1826834460600733*q2*(1/2251799813685248)+431327796402257*q3*(1/18014398509481984)+843667886835955*q4*(1/70368744177664))-q1*(3380075947548081*q1*(1/281474976710656)-259050600068343*q2*(1/281474976710656)-1826834460600733*q3*(1/2251799813685248)+4414049272733425*q4*(1/18014398509481984)))+lamda*q3;
f4 := lamda*q4-(4414049272733425*q1*(1/9007199254740992)+1826834460600733*q2*(1/1125899906842624)+862655592804515*q3*(1/18014398509481984)+843667886835955*q4*(1/35184372088832))*(q2*(8289619202186977*q1*(1/9007199254740992)+3380075947548081*q2*(1/281474976710656)-4414049272733425*q3*(1/18014398509481984)-1826834460600733*q4*(1/2251799813685248))+q3*(1826834460600733*q1*(1/2251799813685248)-4414049272733425*q2*(1/18014398509481984)+843667886835955*q3*(1/70368744177664)-215663898201129*q4*(1/9007199254740992))-q4*(4414049272733425*q1*(1/18014398509481984)+1826834460600733*q2*(1/2251799813685248)+431327796402257*q3*(1/18014398509481984)+843667886835955*q4*(1/70368744177664))-q1*(3380075947548081*q1*(1/281474976710656)-259050600068343*q2*(1/281474976710656)-1826834460600733*q3*(1/2251799813685248)+4414049272733425*q4*(1/18014398509481984)));
f5 := q1^2+q2^2+q3^2+q4^2-1;
ord := tdeg(q1, q2, q3, q4, lamda);
                  tdeg(q1, q2, q3, q4, lamda)
G := Basis([f1, f2, f3, f4, f5], ord);

IsZeroDimensional(G);
                             false
ns, rv := NormalSet(G, ord);
Error, (in Groebner:-NormalSet) The case of non-zero-dimensional varieties is not handled.

Any help please ?

Thank you.

Hi!

I have seen th following procedure to compute the image of the points of [0,1] under the so called Peano space-filling curve (sorry, I have to pasted the code in "text plane" mode):

P[0] := (x, y) -> ((1/3)*y, (1/3)*x);

P[1] := (x, y) -> (-(1/3)*x+1/3, (1/3)*y+1/3); 

P[2] := (x, y) -> ((1/3)*x, (1/3)*y+2/3); 

P[3] := (x, y) -> ((1/3)*x+1/3, -(1/3)*y+1);

P[4] := (x, y) -> (2/3-(1/3)*y, 2/3-(1/3)*x); 

P[5] := (x, y) -> ( (1/3)*x+1/3, 1/3-(1/3)*y));

 P[6] := (x, y) -> (1/3)*x+2/3, (1/3)*y);

P[7] := (x, y) -> (-(1/3)*x+1, (1/3)*y+1/3);

P[8] := (x, y) -> ((1/3)*x+2/3, (1/3)*y+2/3);

peano := proc (t::numeric, depth::integer)

local q, r; global P;

if depth = 0 then return 0, 0 end if;

q := floor(9*t); r := 9*t-q;

return P[q](peano(r, depth-1))

end proc;

Now, I need to use the procedure "peanofun" as a function. For instance, if we define f:=(x,y)->x+y, I need to use (plot, compute, etc) for instance, the function f(peanofun(t,5))

Can you help me with this issue, please?

Many thanks for your time!