Maple 18 Questions and Posts

These are Posts and Questions associated with the product, Maple 18

Using the learning sequence as an alternative to learn problems related to "balance of a body" is shown in this video; thanks to the kindness that Maple offers us in its fundamental programming syntax.

Lenin Araujo

Ambassador Of Maple



I have a 2 dimensional array with variable length. I want to extract the number of elements of each row. For example What is the number of element in x(1,i)?

Hi everyone, I am trying to plot field lines but I face problem in plotting field plot. kindly help me if anyone can do it. Find the code in the attached file. Thank you.

tomleslie 2300Mohsen 20 @Carl Love @Markiyan Hirnyk @maple fan 175 @asa12 405  @Adri van der Meer 1253 

Hi, i want to investigate  chaos for the problem , cantilever beam under random narro band excitation, but the code has errors .the code is this:

      restart:with(plots):      h:=1: Omega:=(0..376):alpha1:=617.2:alpha2:=1.02*10^(8): c:=.002:k:=18.4:  step:=0.1:imax:=376:  for i from 0 to imax do;  Omega[i]:=i*step:   f:=evalf(solve({((-a*Omega[i]^(2)+alpha1*a+3/(4)*alpha2*a^(3)+1/(4)*k*Omega[i]^(2)*a^(3)-(3)/(4)*k*Omega[i]^(2)*a^(3))^(2)+(c*Omega[i]*a^())^(2))=h^(2),a>0}));  ff[i]:=((rhs(f[1]))^(2))/(2):  end do:   l1:=[[Omega[n],ff[n]] $n=0..imax]:  p1:=plot(l1, x=0..3,y=0..1,  style=point,symbol=solidcircle,symbolsize=4,color=red):    jmax:=914: f1:=array(377..914):f2:=array(377..914):f3:=array(377..914):Omega1:=array(377..914):  for j from 377to jmax do;  Omega1[j]:=j*step:   fff:=evalf(solve({((-a*Omega1[j]^(2)+alpha1*a+3/(4)*alpha2*a^(3)+1/(4)*k*Omega1[j]^(2)*a^(3)-(3)/(4)*k*Omega1[j]^(2)*a^(3))^(2)+(c*Omega1[j]*a^())^(2))=h^(2),a>0}));  f1[j]:=((rhs(fff[1,1]))^(2))/(2):f2[j]:=((rhs(fff[2,1]))^(2))/(2):f3[j]:=((rhs(fff[3,1]))^(2))/(2):  end do:   ll1:=[[Omega1[n],f1[n]] $n=377..jmax]:  pp1:=plot(ll1, x=0..10,y=0..1,  style=point,symbol=solidcircle,symbolsize=4,color=red):    ll2:=[[Omega1[n],f2[n]] $n=377..jmax]:  pp2:=plot(ll2, x=0..10,y=0..1,  style=point,symbol=solidcircle,symbolsize=4,color=red):    ll3:=[[Omega1[n],f3[n]] $n=377..jmax]:  pp3:=plot(ll3, x=0..15,y=0..1,  style=point,symbol=solidcircle,symbolsize=4,color=red):       plot({  seq(seq(p1), seq(seq(pp1),seq(seq(pp2),seq(seq(pp3))  },style=point,title=`Pitchfork Diagram`);  Thanks for your help

Dear all,
When I try to simplify the expresion z+conjugate(z) I get, as expected, 2*Re(z).

The problem arises when I try to do the same with a more complicated expression as z1*conjugate(z2)+conjugate(z1)*z2. Trying to simplify this Maple only give me back the expression with nochange. The expresion is equivalent to the first example so I expected to get 2Re(z1*conjugate(z2)).

Does anybody know how to solve this problem?
Thank you in advance!

Dear All

Is it possible to change the position of labels so that they do not overlap with curves of plots? Please see the attached file in which I want to shift label for vertical axes to slightly left and label for horizontal axes to the slightly right.


Hi Dears,

I installed Maple 15 and 18 on my PC. But, when I want to use the PolyhedralSets package the following error is appeared. Is this package only on Maple 2015?


With the bests.

When I use listcontplot, the tickmarks on the axis show the ordinal of the point plotted, so if I have 20x50 points, it shows (1 to 21)x(1 to 51).

Is there any way to reescalate the axis to show the actual units of the points?, so for example (0.3 to 0.7)x(2.2 to 3.7).


hi every one

 im new in maple, i have this problem , the output result (nod) dont disply after the for loop how can i disply it ?( sorry for this question, im new in this prog.)

the 2nd question is the (Typesetting:-delayDotProduct) how can i resolve this???

plz help 


for i from 2 to 3 do for j from 2 to 3 do Typesetting[delayDotProduct](rod, [i, j], true) := s*imp*l*i*f*y*subs*({cx^2 = 1-cy^2-cz^2}, Typesetting[delayDotProduct](rod, [i, j], true)) end do end do;
Typesetting[delayDotProduct](rod, [1, 1], true) := collect(rod[1, 1], cos(theta));


Implementation of Maple apps for the creation of mathematical exercises in

In this research work has allowed to show the implementation of applications developed in the Maple software for the creation of mathematical exercises given the different levels of education whether basic or higher.
For the majority of teachers in this area, it seems very difficult to implement apps in Maple; that is why we show the creation of exercises easily and permanently. The purpose is to get teachers from our institutions to use applications ready to be evaluated in the classroom. The results of these apps (applications with components made in Maple) are supported on mobile devices such as tablets and / or laptops and taken to the cloud to be executed online from any computer. The generation of patterns is a very important alternative leaving aside random numbers, which would allow us to lose results
onscreen. With this; Our teachers in schools or universities would evaluate their students in parallel on the blackboard without losing the results of any student and thus achieve the competencies proposed in the learning sessions.
In these apps would be the algorithms for future research updates and integrated with systems in content management. Therefore what we show here is extremely important for the evaluation on the blackboard in bulk to students without losing any scientific criteria.


Lenin Araujo Castillo

Ambasador of Maple


There seems to be a bug: in some of the cases, the intersection command works, in others it doesn't.

I'd appreciate any hints as to why this is the case.

Hello everyone, 

please need some brain power here to help me with this problem. I'm stacked here and I cannot find a solution. 

I would like to find the intersection between a plane and a quadric surfaces. Please notice that the two equations are written like f(x,y) and g(x,z) respectively. My first problem is that I cannot visualize the two surfaces correctly as implicitplot3d of the quadric surfaces displays nothing in return. 

I want then to find the centre of the ellipse and plot the point and the two surfaces all together in the same 3d plot. Somehow I canont make it. What is wrong with me? 

I am pasting here my MWE. 

Any kind of help would be really appreciated. Thank you so much. 





with( Optimization ):

my_plane := 2024.30587449691-.341275505799078*x-3.89936179341114*y;
my_quadric := -10595.4104931095+6.73749956241827*x+42.1022818380012*z-0.654818649508000e-1*x^2-0.421174257681115e-1*z^2;





my_quadric_plot_3d := implicitplot3d(my_quadric, x=-10000..10000,
                                       y = -100..100,
my_plane_plot_3d := implicitplot3d(my_plane, x=-10000..10000,
                                       y = -100..100,

my_ellipse := CompleteSquare( subs(y = solve(z = my_plane, y), my_quadric));



ellipse(e1, my_ellipse, [x,z]):
XY_centre_e1     := coordinates(center(o,e1)):
my_center_point := [XY_centre_e1[1],XY_centre_e1[2],
                       subs(x=XY_centre_e1[1], y=XY_centre_e1[2], my_plane)]:
my_center_point_plot := pointplot3d(my_center_point, color = "red"):

display(my_center_point_plot, my_plane_plot_3d, my_quadric_plot_3d);






Does any body have example of maple code for solving optimal control problem in deterministic model using Pontryagin's maximum (or minimum) principle?

First 8 9 10 11 12 13 14 Last Page 10 of 64