MaplePrimes Questions

Search Questions:

It is possible to perform the simplest QFT calculations with second quantization, in Maple? Bosons in a box. See attached example. bosons_in_a_box.mw

Sure any general purpose programming language is capable of performing this task with enough effort. What I am interested in is if the physics tools has a standard way of dealing with these calculations. The general impedement when attempting the calculation is that integrations are perfomed by replacements with delta functions or kronecker delta functions, and its not clear how to force the Maple Physics package to recognize this or if that's possible. Part of the problem is that integrations in maple are defined in one dimension at a time where as in QFT the integration element is almost always atleast three dimensional, d^3x or dxdydzy, the later of which can get extremely cumbersome with even a small number of fields under consideration. I don't find much of what I am refering to mentioned in the help pages and I doubt these types of QFT calculations are possible to perform in Maple without addressing these issues.

bosons_in_a_box.mw

Hello everyone!
I have had an issue for the past weeks, where it seems like Maple has a problem identifying the correct units and also sometimes having an issue with defining a variable.

In the picture below you can see I have defined rho, m and tried to Solve V. I get m^3, which is fine but I cannot change the unit in the right bar (see picture).

I even tried to just take square root of my V to see if I could then change units. It seems like it thinks I'm playing with weight.
I even had 2 teachers trying to help me find a solution, without luck.

I have tried executing the whole document and also only bits of the document without luck.

Does anyone know this problem and has a solution?

Thanks in advance! 

- Jacob

Hi 
How can I plot f3=0 for different values of Tch and qc?

1.mw

For the derivation of the results that I need for my research, I use advanced calculus, special functions etc.

Consider a derivation of the form LHS = RHS1 = ... RHSn.

Suppose that I have a LaTeX file describing this derivation, with fine details of the derivation in notes written using pencil on paper. 

To check this derivation for errors, such as a wrong sign, I could retype LHS and RHSn into Maple and check for equality using symbolic algebra.

To provide a more detailed check of this derivation, I could type LHS, RHS1, ..., RHSn into Maple and check the sequence equalities LHS = RHS1, ..., RHS(n-1) = RHSn using symbolic algebra.

However, this is a laborious and, itself, an error-prone process.

Is there some way of submitting my LaTeX file (or its compiled version) and/or an image of my handwritten notes for checking by Maple, line-by-line, using symbolic algebra?

It's disappointing that after all these years, Maple still hasn't implemented a real sinc() function, in which an expression with a trig expression(s) can be entered as an argument, with a return value that is expressed as a sinc() function(s), not a return of sin() functions, which is what I get when I try to define my own sinc() function. I find that the Sinc function in the Dynamic Systems package to not be very useful for general Maple calculations.

The following code converts an expression containing a trig function to a sinc-function expression. How can this code (or some other code) be implemented in Maple by defining a real sinc function that accepts an expression in its argument and returns an evaluated sinc function, not a return value that has sin functions?

eval(expr, {sin = (x -> x*sinc(x)), 
              cos = (x -> (x+Pi/2)*sinc(x+Pi/2)),
              tan = (x -> x*sinc(x)/(x+Pi/2)/sinc(x+Pi/2))});

Attribute: This code was originally posted by  Robert Israel 6577

Thanks for any assistance.

i don't  know  why my graph make a problem and what is issue i did plot  but this time make issue for me which i don't know where is problem there is anyone which can help and even modified the plot?

explore-chaotic.mw

i try to use explor option but is not work like i want and i try to use other way but they so classical i want to ask there is any other way for do plotting of this kind of function?
in second part even plot not working?

 

Bgraph1.mw

Hi everyone

I have the following list:

L:=[[F[1], F[2], F[3], F[4]], [F[1], F[2], F[3], M[1]], [F[1], F[2], F[3], M[2]], [F[1], F[2], F[3], M[3]], [F[1], F[2], F[3], M[4]], [F[1], F[2], F[4], M[1]], [F[1], F[2], F[4], M[2]], [F[1], F[2], F[4], M[3]], [F[1], F[2], F[4], M[4]], [F[1], F[2], M[1], M[2]]

How do I select only the subsets in which the four elements have different indices? For example, I don't want [F[1], F[2], F[3], M[1]] Regards,

Oliveira

How do I fully expand an expression a polynomial with a constant c and variable m ?

i did substitution but my result is so different from the author i think he just take the linear term of theta but i didn't do that so how take just linear term of that function and find unknwon , and how afeter replacing eq(12) inside eq(11) we can remove thus exponential and find w? also i think author did a mistake which the equation 12 is theta(x,t) not Q(x,t)

restart

with(PDEtools)

undeclare(prime, quiet)

declare(u(x, t), quiet); declare(U(xi), quiet); declare(V(xi), quiet); declare(theta(x, t), quiet)

pde := diff(u(x, t), `$`(t, 2))-s^2*(diff(u(x, t), `$`(x, 2)))+(2*I)*(diff(u(x, t)*U^2, t))-(2*I)*alpha*s*(diff(u(x, t)*U^2, t))+I*(diff(u(x, t), `$`(x, 2), t))-I*beta*s*(diff(u(x, t), `$`(x, 3)))

diff(diff(u(x, t), t), t)-s^2*(diff(diff(u(x, t), x), x))+(2*I)*(diff(u(x, t), t))*U^2-(2*I)*alpha*s*(diff(u(x, t), t))*U^2+I*(diff(diff(diff(u(x, t), t), x), x))-I*beta*s*(diff(diff(diff(u(x, t), x), x), x))

(1)

T := u(x, t) = (sqrt(Q)+theta(x, t))*exp(I*(Q^2*epsilon*gamma+Q*q)*t); T1 := U = sqrt(Q)+theta(x, t)

u(x, t) = (Q^(1/2)+theta(x, t))*exp(I*(Q^2*epsilon*gamma+Q*q)*t)

 

U = Q^(1/2)+theta(x, t)

(2)

P := collect(eval(subs({T, T1}, pde)), exp)/exp(I*(Q^2*gamma*`ε`+Q*q)*t)

diff(diff(theta(x, t), t), t)+(2*I)*(diff(theta(x, t), t))*(Q^2*epsilon*gamma+Q*q)-(Q^(1/2)+theta(x, t))*(Q^2*epsilon*gamma+Q*q)^2-s^2*(diff(diff(theta(x, t), x), x))+(2*I)*(diff(theta(x, t), t)+I*(Q^(1/2)+theta(x, t))*(Q^2*epsilon*gamma+Q*q))*(Q^(1/2)+theta(x, t))^2-(2*I)*alpha*s*(diff(theta(x, t), t)+I*(Q^(1/2)+theta(x, t))*(Q^2*epsilon*gamma+Q*q))*(Q^(1/2)+theta(x, t))^2+I*(diff(diff(diff(theta(x, t), t), x), x)+I*(diff(diff(theta(x, t), x), x))*(Q^2*epsilon*gamma+Q*q))-I*beta*s*(diff(diff(diff(theta(x, t), x), x), x))

(3)

 

TT := Q = alpha[1]*exp(I*(k*x-t*w))+alpha[2]*exp(-I*(k*x-t*w))

Q = alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w))

(4)

S := eval(subs(TT, P))

diff(diff(theta(x, t), t), t)+(2*I)*(diff(theta(x, t), t))*(gamma*epsilon*(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^2+(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))*q)-((alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^(1/2)+theta(x, t))*(gamma*epsilon*(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^2+(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))*q)^2-s^2*(diff(diff(theta(x, t), x), x))+(2*I)*(diff(theta(x, t), t)+I*((alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^(1/2)+theta(x, t))*(gamma*epsilon*(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^2+(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))*q))*((alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^(1/2)+theta(x, t))^2-(2*I)*alpha*s*(diff(theta(x, t), t)+I*((alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^(1/2)+theta(x, t))*(gamma*epsilon*(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^2+(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))*q))*((alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^(1/2)+theta(x, t))^2+I*(diff(diff(diff(theta(x, t), t), x), x)+I*(diff(diff(theta(x, t), x), x))*(gamma*epsilon*(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))^2+(alpha[1]*exp((k*x-t*w)*I)+alpha[2]*exp(-I*(k*x-t*w)))*q))-I*beta*s*(diff(diff(diff(theta(x, t), x), x), x))

(5)

Download steps.mw

or for this equation 

steps-2.mw

Prove:
If a is an irrational number, then the function
y(x) = cos(a*x) + cosx is not periodic.

Further question:
Since every irrational number can be represented as the limit of a sequence of rational numbers and all software "only" works with finite representations of numbers, it would be interesting to know at what number length for a this becomes in examples apparent in the plot for y(x).

I'm not sure exactly when this started but maple on one of my linux machines won't print more than one row of a matrix. I get dots like this:

Note that interface(rtablesize) is [10,10]. This happens with both maple2024 and maple2025. 

On any of my other machines, maple works just fine. I have also tried deleting ~/.maple and ~/.maplesoft

Any idea why Maple dsolve can't find solution to this ode? From textbook

The strange thing, it solves if it asked for implicit solution. But the default, will give no solution.

Is this a defect? Should it not have returned the book solution automatically?   How is a user supposed to know the ode has a solution or not, if default call returns no solution?

restart;

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

SupportTools:-Version();

`The Customer Support Updates version in the MapleCloud is 29 and is the same as the version installed in this computer, created June 23, 2025, 10:25 hours Eastern Time.`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1877 and is the same as the version installed in this computer, created 2025, July 11, 19:24 hours Pacific Time.`

restart;

ode:=v(x)*diff(v(x),x) = g;
ic:=v(x__0) = v__0;
sol:=dsolve([ode,ic]);

v(x)*(diff(v(x), x)) = g

v(x__0) = v__0

restart;

ode:=v(x)*diff(v(x),x) = g;
ic:=v(x__0) = v__0;
sol:=dsolve([ode,ic],'implicit');

v(x)*(diff(v(x), x)) = g

v(x__0) = v__0

-2*g*x+v(x)^2+2*g*x__0-v__0^2 = 0

#why did not default call return this?
PDEtools:-Solve(sol,v(x))

v(x) = (2*g*x-2*g*x__0+v__0^2)^(1/2), v(x) = -(2*g*x-2*g*x__0+v__0^2)^(1/2)

Download dsolve_gives_no_solution_sept_2_2025.mw

I have bought Maple Flow 2025 and trying to use it in an effective way. While doing this I noticed some little issues that could be improved.

1. When typing text it would be nice that the cursor goes automatically to the next line when at the border of chosen paper size (e.g. A4).

2. It would be nice that you could separate the containers when they overlap by using 1 button. Now you can align containers vertically but it would be nice if you can do that also horizontally.

I also have a question, is it possible to place a section in another section as you can in Maple or is this not possible? 

I test a lot of them but some of them make a problem i  don't know i am do it in wrong way or the author did wrong i need verifying thus solution of odes specially in case 4 when we have not equal sign how use that?

and case 5 is Weierstrass elliptic function which i don't know how set up and use i think is a on kinf of odes but why they use that sign for this function?

ode-17.mw

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