I'd like to find out what number Maple can still calculate the isprime for, and what number will make it go "No thanks, this is madness - ERROR".

I tried isprime (10^100,000,000 + 2), and it can still calculate the result. No prime of course, that's the whole point of the +2.

I tried isprime(10^10ˆ9 + 2) and it crashed.

But what is the tipping number?

I decided I should write a loop, and let maple write away the last number till it crashes. Is this possible? I think I'll need something like

> for n from 5 to ??? do isprime (10ˆn + 2) and ??? end do

This isn't suitable at all, but I cant figure out how to program "to error" or "write away this number"...

I hope the question is clear and look forward to ideas on how to make this work.

Dear Community,

I have difficulties evaluating a rather complicated infinite integral (see attachment). The result should be something like 12.319. Any idea how to go get this?

Tx for the kind help in advance,

best regards

Andras

Infinite_Integral_2.mw

Plz help me！ How to accelerate the calculation speed the following the program？

restart;
with(PDEtools); declare((u, W)(x, t)); KN := 10;
AFP := proc (C1, C2, C3, C4, H, KN, N) local ADM1, n, lambda, F, i, A0, A, U, W, u, L, R, NL, w, PDE1, IC1, d, Eq1, Eq2, Eq3, LI, trL, tr1, trN, Apr, AprK, trSol, Sol, AD; declare((u, W)(x, t)); ADM1 := proc (n) options operator, arrow; convert(subs(lambda = 0, value((Diff(F(Sum(lambda^i*U[i], i = 0 .. n)), `$`(lambda, n)))/factorial(n))), diff) end proc; A0[0] := F(U[0]); for n to KN do A0[n] := ADM1(n) end do; for n from 0 to KN do A[n] := convert(C1*(diff(subs({seq(U[i] = W[i](x, t), i = 0 .. KN)}, expand(subs(F(U[0]) = U[0]*exp(U[0]), A0[n]))), x)), diff)+convert(C2*subs({seq(U[i] = Diff(W[i](x, t), x), i = 0 .. KN)}, expand(subs(F(U[0]) = U[0]^2, A0[n]))), diff) end do; L := proc (w) options operator, arrow; diff(w(x, t), t) end proc; R := proc (w) options operator, arrow; C3*(diff(w(x, t), x)) end proc; NL := proc (w) options operator, arrow; C1*(diff(w(x, t)*exp(w(x, t)), t))-C2*(diff(w(x, t), x))^2 end proc; PDE1 := proc (w) options operator, arrow; L(w)-R(w) = -NL(w) end proc; IC1 := u(x, 0) = sum(2*(int(sin((d+1/2)*Pi*x/H), x = 0 .. H))*exp(-C4*(d+1/2)^2*Pi^2*t/H^2)*sin((d+1/2)*Pi*x/H)/H, d = 0 .. N); LI := proc (w) options operator, arrow; Int(w(x, t), t = 0 .. t) end proc; tr1 := u-rhs(IC1); Eq1 := LI(lhs(PDE1(u))) = LI(rhs(PDE1(u))); Eq2 := simplify(subs(lhs(Eq1) = tr1, Eq1)); trL := u = add(u[j](x, t), j = 0 .. KN); trN := LI(NL(u)) = Int(Sum(A[i], i = 0 .. KN), t = 0 .. t); Eq3 := subs(trL, lhs(Eq2)) = subs(trN, rhs(Eq2)); Apr[0] := u[0](x, t) = rhs(IC1); AprK := u[k+1](x, t) = -(Int(AD[k], t = 0 .. t)); for i from 0 to KN do Apr[i+1] := value(subs({seq(Apr[m], m = 0 .. i)}, subs({seq(W[m] = u[m], m = 0 .. i)}, subs(k = i, AD[i] = A[i], AprK)))) end do; trSol := {seq(Apr[i], i = 0 .. KN)}; value(subs(trSol, trL)) end proc;

Download AFP.mw

Am here again, pls help me check out this adm code, cant get a result.

below is the attached file

 adm_2_method.mw

Ok i have used a diverse range of programming IDE now, but i am of course alot more comfortable the maple interface simply because of the amount of time i have spent on it, also it was my first "real" introduction to programming, but i was keen to trial and review of all the major CAS out on the market, but for a common purposed package within them they provide, for example the number theory package.

So the question is quite broad, im looking for suggestions for the most unbias way to go about this.

it would need to be a coauthored review with a representative writer for the seasoned user for each one right?

hi...how i can pdsolve this equation numerically or analyticlly?

this equation is time-fractional  equation with generalized Cattaneo model

where

 is the fractional derivative operator considered in the
Caputo sense.

FRACTION.mw

Download FRACTION.mw

Dear Community,

I would like to calculate a rather sophisticated function (see it in the attachmenet) Once I try to evaluate it at certain values I get the error message:

"Error (in pD) illegal use of an object as a name" 

What object Maple means? I guess its a minor mistake, yet I'm lost here, and the online help says "There is no help page available for this error"

Tx for the kind help in advance,

Andras

Bessel_plot.mw 

Is anyone running the 2016.1a patch on ubuntu 16.04?

A fresh installation starts up and works, but after I install the patch the installation seems broken:

cd /usr/local/maple2016/bin
./maple
/usr/local/maple2016/bin.X86_64_LINUX/mserver: error while loading shared libraries: libicuucmpl.so.56: cannot open shared object file: No such file or directory
maple: could not start mserver

it looks like my installation does not have libicuucmpl.so.56. I can see these however:
./bin.X86_64_LINUX/libicuucmpl.so.49.1.1
./bin.X86_64_LINUX/libicuucmpl.so
./bin.X86_64_LINUX/libicuucmpl.so.49

running xmaple has maple open, then stay unresponsive.

Here is the upgrade file, i've been trying to use:
md5sum ./Maple2016.1LinuxX64Upgrade.run
7502caaa65cc623d5d2574823eee9343  ./Maple2016.1LinuxX64Upgrade.run


Has anyone else encountered this problem?
thanks

Dear Pros, I'm a biginer so I have a question about my program.

I have a lot of arrays which are result from 2 while loop. Now, I want creat a matrix from them but i can't. So, could you help me to do it.

For detail: 

V[1]:=[ 1 2 3]

V[2]:=[2 3 4]

V[3]:=[3 4 5]

V[4]:=[2 6 7]

V[5]:=[7 8 9]

...

V[n]

with type of V[i] is a array.

I searched and found a solution by manual to create a matrix as follow:

V_matrix:=<V[1],V[2],V[3]>

but in this case i can't but manual with n=100

please help me to have a Matrix.

Thank a lots.

I wondering if it is able to create a new symbol in maple, like it is possible to create a custom snippet in the pallets. Speceficly I am asking how to create a new symbol/charachter for 2D math. So if anyone knows please tell me :)

Dear all,

I would like to solve the following non linear ODE with Maple, but I am no able. I do not know if it is possible, beccause it is nolinear.

I really appreciate any advice or help. This is the equation:

y'(x) - (Q - x*p0*(exp(alpha-beta*y(x)))/(1+exp(alpha-beta*y(x))))^2=0

thanks a lot

Hey everyone,

I'm using printf to print out the following:

printf("Absolute error = %.10f \nRelative error = %.10f", abs(fractionPart-binarySum), 100*abs(fractionPart-binarySum)/fractionPart)

however, it's tedious always having to edit the two occurrences of "10" decimal places. Is there a way to define a variable like decimals := 10 and then refer to it inside the quotation marks in printf?

Hello everyone,

Is it normal that commands #1 and #2 below do not return the same thing ?

L := [`Norman.Mailer`, `Richard.Brautigan`]

#1

map(u -> convert(u, string), L);

["Norman.Mailer", "Richard.Brautigan"]


#2

`convert/string`~(L)

["`Norman.Mailer`", "`Richard.Brautigan`"]


Side question : I am not really familiar with the tilde operator and I often use map instead.
Does it exist a better practice in these matters ?

Thanks in advance

Hello everybody!

Please help me to solve the attached partial differential equation. I am getting an error. I do have its analytical solution and that works fine.

The error is as follows
Error, (in pdsolve/numeric/plot) unable to compute solution for t>HFloat(0.0):
solution becomes undefined, problem may be ill posed or method may be ill suited to solution

The worksheet is attached hereshortsngle.mw

In feature "Procedures for two animations" posted by Kitonium appears this line, within a procedure.

T1:=(O2*~y+O1*~u)/~(y+u): T2:=(O3*~u+O2*~v)/~(u+v): T3:=(O4*~v+O3*~z)/~(v+z):

What is the meaning of "~y", for instance, and how does it differ from "y"?

I understand the conventional usage of ~ with a name in output to imply that that name is subject to an assumption.

I have been able to find nothing pertinent in the Maple Help information.  Please enlighten me.