How does one import large polynomials with 10k terms or more (from Mathematica) into Maple? I don't think there is a compatiable file format between Mathematica and Maple.

I tried copying and pasting my two-variable Laurent polynomial from Mathematica into Maple. It does create a polynomial in Maple but just 1000 terms is enough to consume 5GB RAM and causes Maple to hang unresponsive. For some reason, just copying/pasting is enough to cause evaluation to occur in Maple without hitting the enter key. The "ready" dot at the bottom screen switches to "evaluating" dot automatically when I paste the polynomial in.

An old Maple book from 2011 states Maple is capable of handling 68,000 term polynomials -- I am sure this number has gone higher up now in 2020. 

with(RegularChains):
with(ChainTools):
with(MatrixTools):
with(ConstructibleSetTools):
with(ParametricSystemTools):
with(SemiAlgebraicSetTools):
with(FastArithmeticTools):
R := PolynomialRing([x,y,z,a,b]):
sys := [a*x^2+b*x+c]:
N := []:
P := [a]: 
H := [x,y,z]:
dec := RealTriangularize(sys,N,P,H,R):
Display(dec, R);

how to write : a <> 0 and (some x)(a*x^2+b*x+c) for quantifier elimination in maple?

is RealTriangularize parameters are after quantifier elimination or it accept quantifer expression?

if not, how to input above statement into real or lazy triangularize and expect return b^2 + 4*a*c;

Hi,

I'm coming from Mathematica since Mathematica uses single-core to compute multiplications of puiseux series, so I'm looking for any software that might do this faster. I don't know Maple -- I hope it can use multiple cores and runs faster.

My calculation is basically multplying out thousands of terms of the following type:

How should I do this? Maple doesn't seem to have multivariate puiseux series data structure. But observe my series is simple -- I have a power series in q and a Laurent series in y.

I could

1. Multiply using 'expand' . My worry is, if I expand product of two series in q^10000 then this is probably slower than Mathematica since 'expand' would compute far more terms than necessary instead of throwing away the extra terms beyond q^10000 in some O(q^10000) notation like Mathematica does.

2. I could use 'series' in q. 

This has the undesirable feature of not expanding out in y, returning 

.

I could then call expand on this, but using both series and then expand makes me worry about runtime.

Any advice?

Here is the minimal working example: 

with(Physics):
Setup(noncommutativeprefix={W,V});
simplify(W^(-1)*(-W^(-1)-W*V^(-2)));

It gives the result  below, which is wrong. Or am I misinterpreting? (Maple 2018)

(-V^2-W^2)*(1/V^2)*(1/W^2)

how i can solve a system of integral equations? thanks for the help.
 

Download fsolve.mw

Hello,

I have a procedure, SoS that returns the value for sos. After compiling the procedure, I get float(undefined) for the return value. How can this be fixed?

Thanks so much for your time.

Download Compile_proc.mw

Hey, to all!

First post here and after searching the web for answers I thought I'll ask it here.

haven't used Maple for a long time so I'm quite new in it.

I've solved a PDE in a numeric and in a analytic way and I want to the find the difference between the results.

I've come up with this:

1.2.mw

Had few tries with a few commands but for no success.

I seem to be able to find the specific value of the analytical solution but not the numeric one although it is solved.

Also I can't subtract to find the difference.

I will appreciate help in any kind,

Thank you!

I'm not sure if it's possible to find a formula with predictive qualities, but I thought I'd ask anyway ...

IOD.mw

         I'd like to draw a graph with its  vertices  of 3d style.  Just like following image：

  I did not find the desired choices for this graph.

Hello All,

I have an autonomous system of ordinary nonlinear equations. In order to investigate it near the fixed points, I would like to reduce it to the normal or hyper-normal form. I see in the literature that some authors developed unique algorithms for the very specific differential systems (e.g., the systems with two/three eigenvalues), which is not my case.  

Maybe you know, if Maple has any specific commands that may conduct such type of reduction?

Thanks in advance,

Dmitry

Hi,

I’m trying to animate the volume generated by rotating the function curve around the Ox axis. the Solidofrevolution command with animation output generate a disk volume

Thanks

QuestionAnimationVolumeOfRevolution.mw

I am solving coupled PDEs under various BCs and ICs on the (x,t) domain, to obtain momentum p((x(t),t) solution plots of the within a MAPLE worksheet. 

However, I am getting JAVA memory overload (> 2GB), which brings MAPLE to a halt when plots are displayed in the worksheet.  Is there a way of writing the graphical plots directly to a file on my computer, so that the plots can later be view separately on my computer screen, thus avoiding JAVA memory overload within MAPLE.

If this is possible, please could someone send me MAPLE code to create external plots...perhaps using simple examples in MAPLE?

Melvin .

with(DEtools); with(plots); with(plots, implicitplot);

LC := [diff(x(t), t) = x(t)*(1-x(t)^2-y(t)^2)+y(t)*((1-x(t))^2+y(t)^2), diff(y(t), t) = y(t)*(1-x(t)^2-y(t)^2)-y(t)*((1-x(t))^2+y(t)^2)]

p1 := DEplot(LC, {x(t), y(t)}, t = 0 .. 100, [x(0) = 1, y(0) = 1])

Error, (in DEtools/DEplot/CheckInitial) too few initial conditions: [x(0) = 1]

Download mm.mwmm.mw

Hello,

I am trying to generate a compiled c function

dc_test := proc( a :: Array(datatype=float[8]), b :: Array(datatype=float[8]), c :: Array(datatype=float[8]) ) 

c[..,..] := a . b;

return 0

end:

cdc_test:=Compiler:-Compile(dc_test, optimize):

a := Array(1..2, 1..2, fill=1.0, datatype=float[8]);
b := Array(1..2, 1..2, fill=1.0, datatype=float[8]);
c := Array(1..2, 1..2, datatype=float[8]);
                           
dc_test(a, b, c);

cdc_test(a, b, c);

This doesn't work because it doesn't recognize the "." symbol. Also if I were to take the transpose, for example Transpose(a), then it doesn't recognize that either. Are there commands for this that I'm unware of or do you have to program them from scratch?

Hi,

How to animate a plot with Shadebetween command?

Thanks

AnimationShadebetween.mw