MaplePrimes Questions

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

 

I Want to plot only the intersection between this Surfaces 

 

z = sqrt(x^2+y^2)

y >= abs(x)

x^2+y^2+z^2 >= 4

x^2+y^2+z^2 <= 9

 

Anyone knows how ?

Thanks!

Greetings,

Maple 15 allows the following syntax

omega := sqrt(w0^2*(1+((z-zf)/z0)^2))

But it does not allow this one 

omega^2 := sqrt(w0^2*(1+((z-zf)/z0)^2))

Why is this so for functions and variables? Is there any way around this, I am really bugged by this issue. The Maple Math told me that the syntax was invalid so I had type the commands directly.

 

Hi,

 

I've encountered a problem with Maple.

 

Is it possible to create a procedure which both displays some plots and also returns a value?

 

When I write a simple example:

 

It plots a circle as expected.

But when I add return:

 

It returns value 2, but no plot is displayed.

Is it possible to do both, display and return a value?

I used to write two copies of procedure, one for plotting and one for returning value.

But thats not "pretty" solution.

 

Thanks for help.

I have a formula, which is called the cube law. It requires input between 0 <= x <= 1. When creating  a procedure I run into all kinds of problems and get nothing but error codes. Trying to use assume works, it puts a tilde mark on the variable x, but doesn't restrict the input. Here is the formula the way it is writtin in a book on calculus:

f(x) := x^3/x^3+(1-x)^3 (0 <= x <= 1).

Trying to restrict the input to be restricted as required by the relational statement has been proven fruitless.

Any help would be greatly appreciated, thanks,

Rolf.

How do I setup/plot this probability

As a thought example, time waiting in a line at cash register.  Let's say 40% of the time you are waiting for ~1 minute while 35% of the time you are served immediately and the rest of the time trailing off waiting up to 15 minutes.

Is this is like a skewed normal distribution?  and how could I set it up?  And then how can I randomly generate some values from this distribution?

I have a characteristic equation. some times It has polar roots . sometimes It has real roots and sometimes both of them.

I want to extract real roots and extract polar roots if they are.

for instance:

q:=m3*r^3+m2*r^2+m1*r+m0:

rot:=solve(q=0,r);

I want to know how can I use if in this part ?

Hallo, I have a list with derrivatives. I'd like to choose elements with highest derrivative.

Is there any procedure in Maple to determine order of derrivative?

Hi there,

I used the implicitplot command to plot an implicit equation as follows:

with(plots):

implicitplot(x^3+y^3-3xy=0,x=-3..3,y=-3..3, scaling=constrained)

but maple13 takes error and specifies {xy} as an extra unknown.

Your help will be appreciated.

Regards

First 1350 1351 1352 1353 1354 1355 1356 Last Page 1352 of 2434