Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi Maple community. How to find the any real root of the complex equation which has been mention in my worksheet.

 

Real_root.mw

Hi. is there a site for using maple online?

 

 

How to solve this integral equation?

rho:=1/(2);

mu:=1/(2);

T:=1/(4)*(t-x)^(1/(2));

E[rho,mu](T):=Sum((T^(k))/(GAMMA(k*rho+ mu)),k=0..5);
                       
eq1 := h(t) = (8/3)*(int((t-x)^(1/2)*E[rho, mu](T)*(1/(sqrt(Pi)*x^(1/2))-(3/32)*x^2+(3/16)*exp(-x)-1), t = -1 .. 1));       

 

Download 1.mw

   

Hi.

I have 7 polynomials r[9], r[6, 6], r[7, 9], r[8, 6], r[9, 6], r[6, 9, 7], r[6, 9, 6] each of 7 variables a[3],...,a[9] and I generate an ideal QQ usinig these polynomials. What Maple tells me quite fast is the follows:

IdealMembership(a[3], QQ) 

true

Also using Groebner as follows
F := [r[9], r[6, 6], r[7, 9], r[8, 6], r[9, 6], r[6, 9, 7], r[6, 9, 6]];
G := Basis(F, tdeg(a[3], a[4], a[5], a[6], a[7], a[8], a[9]));
provides me with the basis quite fast and a[3] is part of it. 

---------------

However, Mathematica is stuck and I have noticed that if I change a[3],...,a[9] to x,y,z,,v,w,t then Maple with the same polynomials is evaluating very long time. 

My question is: what am I doing wrong with subscripts a[3],...,a[9] in Maple? Does Maple 2019 have a problem with it?

I have attached my worksheet. 
Maple_Problem.mw
Would appreciate any help.

BRINKMAN_TRIAL.mw

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/BRINKMAN_TRIAL.mw .
 

Download BRINKMAN_TRIAL.mw

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/BRINKMAN_TRIAL.mw .
 

Download BRINKMAN_TRIAL.mw

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/BRINKMAN_TRIAL.mw .
 

Download BRINKMAN_TRIAL.mw

 

 

hi, can anyone help me?

i am trying to solve a set of boundary layer equation.  when i executed the file, the following message appear.

Suppose I specify a metric in say the differential geometry package or Physics package that has an arbitary function f(x,y) of some of the coordinates as a component. Can I specify Maple to treat the function as real valued when Maple is asked to compute curvature etc for the metric? In some examples I have performed Maple sometimes returns expressions in f(x,y) for components that are in fact zero once one treats f(x,y) as real valued. It would be preferable to have Maple actually return '0' in such cases.

 

I tried the 'assume' command but Maple complained it was being applied to a protected name and I tried using

use RealDomain in simplify (5) end use

after Maple's computations, where (5) lablled Maple's output, but Maple didn't perform any simplification.

 

Hi, see above image, how do I get maple to do the chain rule automatically, for example, if I simply put in the alias V(h), instead of V(h(x)), this last command generates zero. 

 

Also, what is D(V)(h) ? is this equivalent to dV/dh, and if so how would I get maple to present it in this form?

 

Many thanks for your help ! 

 

I am using code edit region for parameters but getting this error, can someone help me to solve this issue.

 

Can anyone correct me, what's wrong with it.

 

Help_solution.mw

 

 

Hello.

 

I have the ode 

x' = 4 * lambda * x^4 - x^6

and I want to plot the phase portrait with x be the horizontal axis and lambda be the vertical axis.

Also, I want to plot the arrows showing the stability of the fixed points as lambda is changing values.

Thanks !

I have two functions, f1(x,y,z) and f2(x,y). I want to use Grid package to parallelize the runs of two functions. I have tried Grid:-Seq() but this way the maple does not seem to finish the computation.

Another thing I have tried so far was using Grid:-Run(), like so:

Grid:-Setup("local", numnodes=2):

Grid:-Run(0, f1, [x,y,z]):

Grid:-Run(1, f2, [x,y]):

It seems that this way i create background processes but I never see them finish. On their own the functions f1 and f2 work fine if used sequentially.

Is there any way to run time in parallel?

Hello,

I'm experimenting the RSA encryption algorithm and in one example, I need to calculate (m^e) mod n where :

n=4142074788597595058238933087635402958555715420375128103710780511

m=65108105099101044032

e=3467995563410655813589059428902496618444213189544232426283402949

and I receive : Error, numeric exception: overflow

Even though these are large numbers, modular arithmetic algorithms should be able to deal with that.

Why is it not the case ?

NB : One of the mathApps in the category Computer Science is "RSA Encryption" and proves to be fairly comfortable with such large numbers, no overflow!

 

hi, why my dsolve command is not working

dsolve.mw

Here is a simple example which generates a loss of the connection to the kernel.

A few words to help you understand what this simplified code is aimed to do.
The inputs are:

  • a matrix (10 by 2 in the example given above), 
  • a "slave" column p,
  • a "master" column q,
  • a list of rangesr.

Let nr the number of elements of r.
The goal is to separate the elements of M[.. ,p] in nr sorted lists L1, ..., Lnr according to the rule

  • M[n ,p] is in Lk  if M[n ,q] is in r[k]  

A notional example is given in the snippet code below.
The variable num contains nr lists, the kth one contains the values of such that M[n ,q] is in r[k] .

As soon as one element of num is empty list the line S := ... generates the Connection to Kernel lost message.

I have made sure that this no longer happens by coding more carefully.

However, since some recent posts have mentioned this type of connection loss (usually difficult to trace), I thought it might be useful to publish this example (a census of connection loss situations might help make Maple more robust)

restart
interface(version)
Standard Worksheet Interface, Maple 2015.2, Mac OS X, December 21 2015 Build ID 1097895

LTB := ListTools:-BinaryPlace:
 
M  := LinearAlgebra:-RandomMatrix(10, 2):
p  := 1:
q  := 2:
r  := [-100..-1, -1..1, 1..100]:
nr := numelems(r):

Y_x, Y_ix := sort(M[.., q] , output=[sorted, permutation]):  
bounds    := map(u -> 1+LTB(Y_x, op(1, u))..LTB(Y_x, op(2, u)), r);
Y_nb      := map(u -> op(2, u)-op(1, u)+1, bounds);
            [2, 0, 8]

num       := [ seq(Y_ix[bounds[k]], k=1...nr) ];
            [[9, 1], [], [3, 2, 6, 7, 10, 4, 8, 5]]

S := [ seq(sort(convert(M[num[k], p], Vector)), k=1..nr) ];

KernelConnectionLost.mw

How can I expand exp(I x) to sin(x)+I*cos(x)

1 2 3 4 5 6 7 Last Page 1 of 1777