Let's say you make a tabulate of plots (using tabulate, because you don't want equal space for your plots as in array-plot etc). Then how should one save or export the final result as one eps or pdf file? When I right click and choose export, it only save one cell of the tabulate alone.

A083758 is the lexicographically earliest infinite sequence of distinct primes such that the concatenation of the first n terms is prime for all n>=1. It starts: 2,3,11,7,41,31..

The first ~ 1000 terms are already known from other codes (Mathematica, PARI), and it so far primes 13,47,61... have not yet been seen in the data. Progress beyond this point is difficult because of the rapid growth of the concatenations and uncertainties about the primality of such large numbers.

I would like to know if it could be possible to have a Maple code (operable in the 2017 version) which could produce terms significantly the first 1000 terms, and including a mechanism for flagging (at the end of any given run) the smallest prime not yet seen in the sequence?

Best regards

David. 

Hi, 

I have an ODE system with a free parameter P, that I can only solve numerically.
To simplify let's suppose this parameter has only 2 values P1 and P2.
Let SOL1 the solution of this system when P=P1 and SOL2 its solution for P=P2:

Again for a sake of simplification let's suppose the independent variable is t and the dependent one is x(t).

How can I plot the difference x1(t)-x2(t) where x1(t) is the solution corresponding to SOL1 and x2(t) the solution corresponding to SOL2
(sorry I screwed something up and changed the config of my keyboard making the interrogation point unaccessible)

Here is a detailed example

Thanks in advance
 

Download odeplot_differences.mw

Hi!

Does it exist a command in the LieAlgebras package to find the list of normal subalgbras of a given finite dimensional Lie algebra?

Thank.

Hi!

I used the command "Decompose" for the following trivial 2 dimensional abelian Lie algebra:

This command returns an error:

Alg:=[D_x,D_y]:
Alg_L:=LieAlgebraData(Alg):
DGsetup(Alg_L):
Decompose();

Error, (in print/_DG) final value in for loop must be numeric or character

Why this issue?

Thanks

Hi

How to associate the Pi curve with its Qi name, ( for example P1 is the name of the curve Q1) starting from the RandomTools random command

Thanks for ideas

Question.mw

Hello, I have a sipmle wave PDE with dirichlet BC,I am trying to plot the solution , so I us the following code try to solve it at first,

pde := diff(u(x, t), t $ 2) - diff(u(x, t), x $ 2) = 0:

h := piecewise(2 < x and x <= 3, 6 - 2*x, 0):

ic := u(x, 0) = h, D[2](u)(x, 0) = 0:

bc := u(0, t) = 0:

solution := pdsolve([pde, ic, bc])

The problem is the following code give me a zero solution, but it should't be according to the d'Alembert formula.

I am not sure where is problem in my code, could someone help me on this? thanks.

I have the first and second fundamental forms of a surface. The Bonnet theorem tells me that this is enough to determine a surface in Euclidean space. Has this theorem been implemented in Maple?

Dear friends, please, I would like to ask for your help with the following problem. 

When I try to compile a code like the following 

P:= proc(n)
m:= floor(log[10](n));
A:= Array(1..m);
end proc:

I get the error "Error, (in CodeGeneration:-IssueError) cannot analyze non-integer range boundary m". 

Could you please assist me in overcoming this problem? Many thanks in advance. 

This is my first time using Maple and I want to use Maple to check the exact solution to a pde.

My pde is  -e(u_xx + u_yy) + a*u_x = x(1-p^2 *x)*sin(ly) and the exact solution is u(x,y) = (Ax^2 + Bx +C)*sin(ly)

where A = -p, B = (1-2aA)/p, C = (2eA - aB) and p = el^2

I tried the following in Maple,

pde := -e*(diff(u(x, y), x, x) + diff(u(x, y), y, y) )+ a*diff(u(x, y), x) - x*(-p^2*x + 1)*sin(l*y) = 0

and i get the result below which I don't understand.

How do I do the right thing in order to check my exact solution.

Thanks

1. Write a procedure that enables to compute the volume of a cone. Remember that a cone volume can be calculated using the following equation:

𝑉𝑜𝑙𝑐𝑜𝑛𝑒 =1/3𝜋𝑟2ℎ
with r and h generated randomly, both in the range from -100 to 100.

If Maple returns negative figures for r and h, ask it to replace these values with their absolute values.

If Maple returns 0 for r and h, Maple asks directly the user to input desired values.

2) With reference to point 2, you have the following matrix:

[
−65 1 32
1 −12 3

12 −2 2
]
Write a procedure that chooses randomly two numbers from the above matrix and assigns them to r and h. Then, it compares the two cone volumes and returns the largest value.

I was plotting some 3d stuff and wanted to animate it. I used orientation in the plot3d to get the orientation I wanted. When trying to animate them the orientation is ignored.

Specifically though I'm suing Pk = plot3d(...);
then animate(display, [P1, ..., Pn]);
 

(more or less, my real code is more complex but that is the idea).

Where can I tell animate to use a specific orientation? I tried adding it to the options of display but it didn't work either. Seems animate and display completely ignore orientation.

How can I remove the extra whitespace around a 3D plot and also scale it so it takes up quite a bit of the window so I can see more detail?

If I use a density plot then I can set the size using

_SP := proc(P,sz::[posint,posint] := [1200,500])
  op(0,P)(remove(type,[op(P)],'specfunc(ROOT)')[],
          ROOT(BOUNDS_X(0),BOUNDS_Y(0),
               BOUNDS_WIDTH(sz[1]),BOUNDS_HEIGHT(sz[2]))):
end proc:

to get exactly what I want. With 3D plot there is a lot of whitespace. I'd rather use plot3d since it is more confirgurable but I need to be able to remove all that additional whitespace(which I guess has to do with rotation in 3D, but here I'm not rotating it, although I might want to in some cases).

I am trying to solve a differential equation numerically. But I am facing some difficulties. It is saying that ODE system has a removable singularity at r = 0. How can I remove the singularity?  Please check the attached figure and try to help.

hello i am trying to remove "[ ]" from an expression so it can be used inside a solve function.

the problem is (eqs) has an "[ ]" outside the expression so it lookslike this "[a = c, a = d, 4*a = e, b = 2*e, b = c + 2*d + 2*f]"

how can i remove [] from this expression?

code

eq := a*(K + Mn + 4*O) + b*(H + Cl) = c*(K + Cl) + d*(Mn + 2*Cl) + e*(2*H + O) + 2*f*Cl;
elems := [K, Mn, O, H, Cl];
eqr := collect(expand(rhs(eq)), elems);
eql := collect(expand(lhs(eq)), elems);
eqs := zip(`=`, map2(coeff, eql, elems), map2(coeff, eqr, elems));
solve([eq, eqs, a = 1]);

thanks