MaplePrimes Questions

Search Questions:

Latest Questions Latest Questions Feed

Hi. I have a commutation table, that I've got as output of Greg Reid algorithm.

Can I use to define Lie algebra for components DifferentialGeometry & Lie Algebras?

Worksheet is attached...


In the following problem though b and c are same (except the way denominator 2 is hanfled), command ' a-b ' readily answers zero, but a-c not so. Why? Only on condition of assumption real it gives zero!

a := (1/2)*(kappa*omega^2+omega^3)*(Y+(1/2)*(-sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(kappa*omega^2+omega^3))^2/omega:

b := (1/2)*(kappa*omega^2+omega^3)*(Y+(1/2)*(-sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(kappa*omega^2+omega^3))^2/omega:




c := (1/2)*(kappa*omega^2+omega^3)*(Y+(-sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(2*(kappa*omega^2+omega^3)))^2/omega:








Why the answer is not given as zero?





What difference therms b and c make for Maple? Are they not same?

Ramakrishnan V

On the Physics Research and development page there are only 3 versions of the Physics package available.  What happened to the earlier final updates for Maple 16, 15 etc...

What is the latest one that can work on Maple12?  Was there one for M12?  What is the earliest versions of Maple the Research & Development Physics packages can work on?

are these quations in wave form for harmonic ocillator?

i can not plot them

plot(exp(I*2*t) - exp(-I*2*t), t=0..10);

plot(cos(t) - I*sin(t), t=0..10);


why same solution but their plot graph are different?

plot(cos(t) + sin(t), t=0..10);
complexplot(exp(I*t), t=0..10);


Hello all, first time posting here.  Anyways, I am trying to make a velocity vector plot from known u and v velocity components at known x and y locations i.e. at (0.1, 0.1), u=-0.5 m/s and v=0.43 m/s.  My calculations have been performed with a FORTRAN code, and the results output in matrix format.  What I am left with are 4 vectors: x location, y location, u velocity component, and v velocity component.  I know that MATLAB has a specific function, called "quiver," for creating a velocity vector plot from known data points (through the use of vectors containing those points).  I have trouble believing that there is not an equivalent function, however convoluted it may be, in Maple to create these plots.  I have tried to manipulate the fieldplot, vectorplot, and arrow functions to no avail.

Any ideas on how to create a velocity vector plot from known data points?  I am using Maple 18.

Here is a link to an example of what I am looking for:

Thanks, Brandon.

I have used data and the curve fitting procedure to find a curve and I now have to put error bars on the graph to complete my task? I know that there is an errorplot function but I cant seem to find a way to do it with the curvefitting thing?


I have attached the whole worksheet, I hope someone can help! Thanks


Even if there was some method of exporting the graph to another medium that would also be very helpful


I use Maple18 on unix 64. I have the MAPLE and the LD_LIBRARY_PATH properly set up to use OpenMaple. When I link my code to OpenMaple library I get segmentation fault when I try to debug using gdb. In order to isolate the problem, I created a simple file that contains the following:


#include <string>

#include <stdio.h>


int main()


int x = 10;

std::string abc;

abc = "testString";

printf("this is it %s\n", abc.c_str());



When I compile the code using g++ -g -o temp I can run ./temp and step into the code using gdb with no problems.


If I link to Maple Library g++ -g -o temp -L/opt/maple18/bin.X86_64_LINUX –lmaplec


I can run the code ./temp but I can not step into it using gdb. I get the error as shown below when I try to step over std::string abc.




GNU gdb (GDB) 7.8

Copyright (C) 2014 Free Software Foundation, Inc.

License GPLv3+: GNU GPL version 3 or later <>

This is free software: you are free to change and redistribute it.

There is NO WARRANTY, to the extent permitted by law.  Type "show copying"

and "show warranty" for details.

This GDB was configured as "x86_64-unknown-linux-gnu".

Type "show configuration" for configuration details.

For bug reporting instructions, please see:


Find the GDB manual and other documentation resources online at:


For help, type "help".

Type "apropos word" to search for commands related to "word"...

Reading symbols from temp...done.

(gdb) b main

Breakpoint 1 at 0x40091f: file, line 6.

(gdb) r

Starting program: xxxxxxxxxxxxxxxxxxxxxxxxxxx

[Thread debugging using libthread_db enabled]

Using host libthread_db library "/lib/".


Breakpoint 1, main () at

6         int x = 10;

(gdb) n

7       std::string abc;

(gdb) n

BFD: reopening /tmp/xf-dll/xf-1274259bc6684076beea8652f625b03bf2133.tmp: No such file or directory


BFD: reopening /tmp/xf-dll/xf-1274259bc6684076beea8652f625b03bf2133.tmp: No such file or directory


BFD: reopening /tmp/xf-dll/xf-1274259bc6684076beea8652f625b03bf2133.tmp: No such file or directory


Can't read data for section '.eh_frame' in file '/tmp/xf-dll/xf-1274259bc6684076beea8652f625b03bf2133.tmp'

0x00000000004007a8 in std::basic_string<char, std::char_traits<char>, std::allocator<char> >::basic_string()@plt


(gdb) n

Single stepping until exit from function _ZNSsC1Ev@plt,

which has no line number information.

Segmentation fault

 with small step rather than two points?

I had writen my questions into this code And I think there's maybe something wrong with A=[0..130]

I think you can understand what I want by reading this code but if you are confused please let me know.

Dear all

I have query related to generating commutator table from given Lie algebra spanning symmetries of PDE, I have ten vectors that I need to commute with each other to constitute commutator table. I have calculated all the commutations manually using "symmetrycommutator" command.

My queries is, can I convert this operation into short loop type program which automatically calculate cummutation of vectors from 1 to 10.

The Maple file is attaches with this query for reference.

Hello All,


I am doing research for my master's thesis and I am trying to type a PDE into maple and have been struggling for about 2 weeks now. I am new to Maple but have had a lot of Matlab experience. I've seen a couple of ways to do variable changes and things but I'm still.


These are the substitutions I'm trying to do.

into this equation, but and some of these terms are 0 (x*d/dt are 0, P, and M_e, and q_ye).

 So it is basically a Euler-Bernoulli beam with a free vibration.

I tried to use dchange to do the substitutions but I am having trouble defining ...

 or doing both



 I see why there is an error but I don't know how to fix it.


I've also tried

Which is close but still doesn't look quite right. because d/dt(omega*t/2) -> omega/2


Any thoughts?


Thanks in advance

Dear all,

I am using Maple to perform numerical integrations. When the final index in the loop is set to 5, the computation is fast and the results are quickly delivered. When I set a number higher than that, even 6, the program gets really slow and often crashes.

I herewith attach the script I use to generate the results. I guess that there should be a problem of memory management and I tried to use gc() as suggested in some forms but without success. I would appreciate it if someone here could explain the reason behind the problem.

Thank you,

restart; Ts := 1.; sigma := 1.; C := 1.; B := 2./(1+C); with(inttrans); beta := B*Ts*omega; assume(Tb > 0); assume(u >= 0); FzzS := -(3/2)*u^3*((2*u+I*beta)*(exp(2*sigma*u)+sigma^2*exp(2*u))-4*sigma*u*((1+sigma)*u-1))/((2*u+I*beta)^2*exp((2*(1+sigma))*u)-4*u^2*((1+sigma)*u-1)^2); InvFzzS := simplify(invfourier(FzzS, omega, t)); logTimeMin := -2; logTimeMax := -1; NumSteps := 6; logTimeStep := evalf(1.0*(logTimeMax-logTimeMin)/NumSteps); curdirectory(); A := matrix(NumSteps, 2); T1 := Array(1 .. NumSteps); AF := Array(1 .. NumSteps); for i to NumSteps do logTime := evalf(i*logTimeStep+logTimeMin); curTime := evalf(10^logTime); A[i, 1] := curTime; A[i, 2] := evalf(Int(eval(InvFzzS, t = curTime), u = 0 .. infinity, epsilon = 10^(-5))); T1[i] := A[i, 1]; AF[i] := A[i, 2] end do







here a comlicated formula,how i simplify

thanks  a lot.


f := (kappa*omega^2+omega^3)*(Y+(-sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(2*(kappa*omega^2+omega^3)))^2/(2*omega)+(-kappa*omega^2+omega^3)*(X+(sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(2*(-kappa*omega^2+omega^3)))^2/(2*omega)+(Omega*N*cos(theta[2])*omega+Omega*N*cos(theta[1])*omega-P__X^2*kappa+P__X^2*omega+P__Y^2*kappa+P__Y^2*omega)/(2*omega)-(sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)^2/(8*omega*(-kappa*omega^2+omega^3))-(-sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)^2/(8*omega*(kappa*omega^2+omega^3))







    f is a complicated function,i want to make it more simplify,but i want to keep square style,

 let coefficients of X and Y keep one unit,and simplify terms  containd special symbol of omega



it what i wanted.

How using Maple one can obtain commutator and adjoint table for given set of infinitesimal generators spanning Lie symmetries for PDE?



The plot command outputs a graph of symbols where the line is made up of "H"s, vertical lines on the axes are "+"s, horizontal lines on the axes are "-"s, and intercepts are "*"s.


How do I fix this problem so that the graph displays a graph normally?

3 4 5 6 7 8 9 Last Page 5 of 1234