My Maple is apparently getting rusty!  Question: how do I recognize a finite constant?

type(Pi, finite) is false (by design, apparently), while type(infinity, constant) is true.  So now I don't know quite what to 'ask' maple!

i use a not good example's polynomials to illustrate the idea

u1:=3*a*c+b^2+c;
u2:=7*a^5*b*c^3;
u3:=a+3*b^2;
T := lexdeg([a,b,c],[e1,e2, e3]);
GB := Basis([e1-u1, e2-u2, e3-u3],T);
result := NormalForm(a*b+b*c, GB, T);

now result is to express a*b+b*c in terms of e1, e2, e3 which represent u1, u2, u3.

is it possible to use preimage to find possible u1,u2,u3 if unknown u1,u2,u3 and given known eqx?

How to use preimage to find possible eqx if given known u1,u2,u3 to find eqx in NormalForm(eqx, GB, T)?

what i confused in code below is that if i know it in terms of -1*e1+2*e2+*e3

it already can be used to find eqx, it seems reasonable to put source1list as unknown to find eqx  or find unknown eq1, eq2, eq3 if given known eqx.

with(RegularChains):
with(ConstructibleSetTools):
source1 := PolynomialRing([e1, e2, e3]);
target1 := PolynomialRing([a, b, c]);
source1list := [...];
target1list := [eq1, eq2, eq3];
cs := PolynomialMapPreimage(target1list, source1list, source1, target1);
Info(cs, source1);
[[e1-e2-e3], [1]]
-1*eq1+2*eq2+*eq3;

Maple 15 does not allow me to do even small changes after executing even the smallest bit of command. Is there a way to overwrite this default setting and force Maple to undo? Do the newer versions of Maple provide this functionality and what is the purpose of undo getting deactivated after a command, at least it should be enabled for the input?

Hi, I am using Maple 18 and struggling with plotting Newton's Method.

I am wanting use the function f(x)=x^3 +cx + 1 where c is a parameter and uses 100 parameter values between -2 and 0, with 100 iterations of each parameter.

Any help would be brilliant.

Thanks in advance,

Neil

deqn4 := {4*eta^4*(diff(f(eta), eta, eta, eta))*(diff(f(eta), eta))-12*eta^3*(diff(f(eta), eta, eta, eta))*f(eta)+12*eta^4*(diff(f(eta), eta, eta))*(diff(f(eta), eta))-36*eta^3*(diff(f(eta), eta, eta))*f(eta)-24*eta^3*(diff(f(eta), eta))^2+12*eta^2*(diff(f(eta), eta))*f(eta)-36*eta*f(eta)^2-4*eta^3*(diff(f(eta), eta, eta))*(diff(f(eta), eta))-12*eta^2*(diff(f(eta), eta, eta))*f(eta)+12*eta^3*f(eta)*(diff(f(eta), eta))+36*eta^2*f(eta)^2-4*eta^3*(diff(f(eta), eta, eta, eta, eta))-12*eta^3*(diff(f(eta), eta, eta, eta))+12*eta^2*(diff(f(eta), eta, eta))-12*eta*(diff(f(eta), eta))+24*f(eta) = 0};
/
| 4 / d / d / d \\\ / d \
< 4 eta |----- |----- |----- f(eta)||| |----- f(eta)|
| \ deta \ deta \ deta /// \ deta /
\

3 / d / d / d \\\
- 12 eta |----- |----- |----- f(eta)||| f(eta)
\ deta \ deta \ deta ///

4 / d / d \\ / d \
+ 12 eta |----- |----- f(eta)|| |----- f(eta)|
\ deta \ deta // \ deta /

3 / d / d \\
- 36 eta |----- |----- f(eta)|| f(eta)
\ deta \ deta //

2
3 / d \ 2 / d \
- 24 eta |----- f(eta)| + 12 eta |----- f(eta)| f(eta)
\ deta / \ deta /

2
- 36 eta f(eta)

3 / d / d \\ / d \
- 4 eta |----- |----- f(eta)|| |----- f(eta)|
\ deta \ deta // \ deta /

2 / d / d \\
- 12 eta |----- |----- f(eta)|| f(eta)
\ deta \ deta //

3 / d \ 2 2
+ 12 eta f(eta) |----- f(eta)| + 36 eta f(eta)
\ deta /

3 / d / d / d / d \\\\
- 4 eta |----- |----- |----- |----- f(eta)||||
\ deta \ deta \ deta \ deta ////

3 / d / d / d \\\
- 12 eta |----- |----- |----- f(eta)|||
\ deta \ deta \ deta ///

2 / d / d \\ / d \
+ 12 eta |----- |----- f(eta)|| - 12 eta |----- f(eta)|
\ deta \ deta // \ deta /

\
|
+ 24 f(eta) = 0 >
|
/
init4 := { {f(0) = 1, (D(f))(0) = 0, ((D@@2)(f))(0) = 0, ((D@@3)(f))(0.1e-2) = 0} };
{{f(0) = 1, D(f)(0) = 0, @@(D, 2)(f)(0) = 0,

@@(D, 3)(f)(0.001) = 0}}
dsol4 := dsolve(deqn4 union init4, numeric, method=rkf45);

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

Hi,

I used to use windows 32 bit, but I have 64 bit windows now. I installed the MapleSim 6.4 and Maple 18 and try to run the simulations that I created with the 32 bit windows. I have this error 'Unable to compile (rc=1), please try again, and if that fails verify your Windows compiller installation'. Could you please let me know what should I do to eliminate this problem?

Best

Onder

Hi,

I have a eight different size lists.

I want to put it into the Array with simplest and clear way as possible respect to first column and containg symbols.
I enclose workseet.

Could you help somone?.
Thank you in advanced

wzel

fill_the_array.mw

i want to plot a arc that starting angle and finishing angle are variable

but it didn't show the arc there is my code

@Markiyan Hirnyk 

how to calculate for module case

M := [[x*y,y,x],[x^2+x,y+x^2,y],[-y,x,y],[x^2,x,y]];

i find maple help using regular chain library, but it has ideal source and ideal target,

when compare with singular, seems different, i am not sure singular code whether correct or not.

so, would like to ask how to do this in maple

if you know singular, i would like to know too.

ring r = 32003, (x,y), lp;
setring r;
ideal Z;
ideal i = x-x2-2*x*y+2*x2*y-2*x*y2+y+y2, 1-x2-2*x*y+2*x2*y-2*x*y2+y2;
map phi = r, i;
ideal i1 = preimage(r,phi,Z);
i1;
ideal i2 = preimage(r,i,i);
i2;
ideal i3 = preimage(r,i,Z);
i3;
ideal i4 = preimage(r,Z,i);
i4;

would like to apply to find Cohen Maculay

While(NewKer <> 0)

Mapping = Basis(M – N) ^ K[M]

NewKer = ker(Mapping, M, N)

N = M

M = NewKer

If IsCM(NewKer) = true then

    NewKer

End if

Do

Hello,

      I would like to solve a system of 9 nonlinear equations, with the constraints on all 9 variables to be that they are nonnegative. How can I do this?

My code is below - I am trying NLPSolve and have tried solve, but am getting stuck.

with(Optimization);

restart; eq1 := 531062-S/(70*365)-(.187*(1/365))*(H+C+C1+C2)*S/N = 0;eq2 := (4/365*(T+C))*S/N-(.187*(1/365))*(H+C+C1+C2)*T/N-(1/(70*365)+1/(5*365))*T = 0; eq3 := (.187*(1/365))*(H+C+C1+C2)*S/N-(4/365)(T+C)*H/N-(1/(70*365)+1/(4*365))*H = 0; eq4 := (.187*(1/365))*(H+C+C1+C2)*T/N+(4/365*(T+C))*H/N-(1/(70*365)+3/(8*365)+.2*(1/365)+.1)*C = 0; eq5 := .1*C-(1/(70*365)+1/(4*365)+1/60+.5)*C1 = 0; eq6 := (1/60)*C1-(1/(70*365)+1/(4*365)+1/210+.5)*C2 = 0; eq7 := .5*C1-(1/(70*365)+1/60+0.1e-2)*CT1 = 0; eq8 := .5*C2-(1/(70*365)+1/210+(1/9)*(0.1e-2*7))*CT2+(1/60)*CT1 = 0; eq9 := N-S-T-H-C-C1-C2-CT1-CT2 = 0; soln := NLPSolve({eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9}, {C, C1, C2, CT1, CT2, H, N, S, T}, assume = nonnegative);

Dear All,

I need your help to plot the phase portrait using DEtools[DEplot]  these are the lines of the code. But when I make RUN, there is an error. I need your help to fix the error. Many thinks.

r1:=1; r2:=1; q2:=2; q1=0.5;  a1:=1;

Sys1 := {diff(N(t),t) = r1*N*(1-N/q1)-P*N/(1+N), diff(P(t),t) = r2*P*(exp(-a1*P)-q2)+P*N/(1+N)};

DEtools[DEplot](Sys1,[N(t),P(t)],t=-10..10,N=0..2,P=0..2);

Many thinks

The problem is The square root of 16-x^2 over the interval [0,-4]  0 being the upper bound, -4 being the lower bound.  I have solved 3/4s of this problem but I don't understand what they mean by "Solve the definite integral exactly by geometry". 

Update:

thanks for the answers below, I can now get listing of content of an .mla like this:

restart;
with(LibraryTools):
LibLocation:=cat(kernelopts(mapledir),"/lib/maple.mla");
ShowContents(LibLocation);

But could not yet find a command that actually lists the standard *.mla libraries that come with Maple. I have to go to C:\Program Files\Maple 18\lib and copy the name of each file with .mla extention in there to get its content. Is there a way to obtain a list of all the .mla files or the standard libraries so that to automate all of this? Otherwise, I'll just copy the names of all the mla files in the above folder and make a list of them manually.

original question:

in Maple, I can do ?packages to get a help page showing list of packages available. But I need a way to do this in code, and obtain result as list or set. Next I'd like to go over each package and also obtain list of all symbols and functions/proc in this package.

For example, once I get list of packages, I'd like to get list of functions in say algcurves package, so I should get ["AbelMap", "algfun_series_sol",.....].  The list can be just strings. That is ok, I am just trying to enumerate all packages and build-in commands and symbols in Maple. May be make a small report of it for documentation purposes and to get better idea what is in Maple.

What are the commands I need to get such lists?

I am using Maple 18.02 on windows.

Hello, Maple wizards.

I have a problem with Tasks for parallel computations. In each task my programm should execute this command:

>B:=SignalProcessing:-DFT(A)

"A" and "B" are Arrays (1..m). "m" is integer. I have 8 tasks. Each task is a loop ("for i from 1 to n do..." n is integer), where my command is executing. "A" and "B" are local for each task. In this case each task should execute my command n times.

So, I obtain this:

"Error, (in SignalProcessing:-Engine:-DFT) attempting to assign to `FwdDFT` which is protected"

What is it? In Maple Help I read that ">SignalProcessing:-DFT" is Thread-Safe Function.