Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I want to shade the area where the two polar curves overlap. The first curve is r=2 and the other curve is r=2(1-cos(theta)). How do I do this?

I want to find the surface area of this parametric curve revolving around the x-axis. I was able to plot the 2D rendering but I want to show the plot so that I can see the surface area. I suspect that would be the 3D rendering.

x=cos(t), y=2+sin(t), 0<=t<=2pi, x-axis

I see this error once in a while

Error, (in SolveTools:-CancelInverses) assertion failed, simplify should not leak _Z variables out from RootOfs 

When I run a long script. Unfortunately, I have not been able to make a small example to reproduce it. So it is random. I've seen it about 5-6 times in the last 2-3 weeks, over 100 runs during this time.

Just wanted to see if anyone saw this before. The strange thing it is random. I can re-run the same script, or restart it from the point where it failed,  with no changes anywhere, and this error do not show up again for days, and suddenly it shows up again.

Must be related how to Maple internal cache memory happened to be arranged at that time. I noticed that a number of things in Maple seem to happen at random times, such as a random hang of the server at different times and places. I could never understand why this happens in Maple.

If there is only one thing Maplesoft could do for the next version, it will be to improve the robustness of its software.

Update March 7, 2024.

FYI, examples giving this error using Maple 2024 using different functions


 

interface(version);

`Standard Worksheet Interface, Maple 2024.0, Windows 10, March 01 2024 Build ID 1794891`

restart;

kernelopts('assertlevel'=2):

ode:=(exp(x)*sin(y(x))-3*x^2)+(exp(x)*cos(y(x))+y(x)^(-2/3)/3)*diff(y(x),x)=0;
DEtools:-odeadvisor(ode,y(x))

exp(x)*sin(y(x))-3*x^2+(exp(x)*cos(y(x))+(1/3)/y(x)^(2/3))*(diff(y(x), x)) = 0

Error, (in SolveTools:-CancelInverses) assertion failed, simplify should not leak _Z variables out from RootOfs

ode:=diff(y(x),x)*y(x)/(1+1/2*sqrt(1+diff(y(x),x)^2))=-x;
dsolve([ode,y(0)=3],y(x))

(diff(y(x), x))*y(x)/(1+(1/2)*(1+(diff(y(x), x))^2)^(1/2)) = -x

Error, (in SolveTools:-CancelInverses) assertion failed, simplify should not leak _Z variables out from RootOfs

 


 

Download leak_Z_maple_2024.mw

 

I have reported this to Mapesoft.

 

Any chance to get a dataplot with units in one or both axis?

Or is this another thing where units are cause of trouble?

plotunits1.mw

Hi Everyone!

I would like your help again.

Considering a Free-Pinned-Pinned-Free beam (page 88 in the pdf file). In case of a matrix 12x12 how could I find the coefficients (a1, a2, a3, a4, b1, b2, b3, b4, c1, c2, c3, c4) using MAPLE in order to plot the mode shapes of the Figure 3.22 (a) (page 70 in the pdf file)? in case of a matrix 16x16 and  20x20, the procedure is the same?

I tried to plot the mode shapes but I failed because I believe they should be similar to Figure 3.22 (a) (page 70 in the pdf file)

 

sys_GE := {-a1 + a3 = 0, -a2 + a4 = 0, 0.9341161484*a1 + 0.3569692163*a2 + 1.519943120*a3 + 1.819402948*a4 = 0, 0.6669014188*b1 - 0.7451459573*b2 + 5.530777989*b3 + 5.620454178*b4 = 0, 0.9938777922*c1 - 0.1104849953*c2 + 62.17096851*c3 + 62.17901032*c4 = 0, -0.9341161484*a1 - 0.3569692163*a2 + 1.519943120*a3 + 1.819402948*a4 + 0.9341161484*b1 + 0.3569692163*b2 - 1.519943120*b3 - 1.819402948*b4 = 0, -0.3569692163*a1 + 0.9341161484*a2 + 1.819402948*a3 + 1.519943120*a4 + 0.3569692163*b1 - 0.9341161484*b2 - 1.819402948*b3 - 1.519943120*b4 = 0, 0.3569692163*a1 - 0.9341161484*a2 + 1.819402948*a3 + 1.519943120*a4 - 0.3569692163*b1 + 0.9341161484*b2 - 1.819402948*b3 - 1.519943120*b4 = 0, -0.7451459573*b1 - 0.6669014188*b2 + 5.620454178*b3 + 5.530777989*b4 + 0.7451459573*c1 + 0.6669014188*c2 - 5.620454178*c3 - 5.530777989*c4 = 0, -0.6669014188*b1 + 0.7451459573*b2 + 5.530777989*b3 + 5.620454178*b4 + 0.6669014188*c1 - 0.7451459573*c2 - 5.530777989*c3 - 5.620454178*c4 = 0, 0.7451459573*b1 + 0.6669014188*b2 + 5.620454178*b3 + 5.530777989*b4 - 0.7451459573*c1 - 0.6669014188*c2 - 5.620454178*c3 - 5.530777989*c4 = 0}

 

solve(sys_GE, {a1, a2, a3, a4, b1, b2, b3, b4, c1, c2, c3, c4});
        {a1 = 0.1014436637 c4, a2 = -0.1143870369 c4, 

          a3 = 0.1014436637 c4, a4 = -0.1143870369 c4, 

          b1 = 0.07095134140 c4, b2 = -0.1260395712 c4, 

          b3 = 0.1510591164 c4, b4 = -0.1737777468 c4, 

          c1 = 0.2102272829 c4, c2 = -0.2816561313 c4, 

          c3 = -1.003990622 c4, c4 = c4}

26_06_2020_Transcedental_equation_matrix_12x12_artigo_2.mw

Hello everyone, can anyone help me to solve this problem? noted I want to solve a couple of odes and need to find involved integration constant i.e c1 to c4. 
 

restart

diff(q(y), y, y)-A*q(y) = B*(P*y+c1)

diff(diff(q(y), y), y)-A*q(y) = B*(P*y+c1)

(1)

NULL

bc := q(-sigma) = 0, q(sigma) = 0

q(-sigma) = 0, q(sigma) = 0

(2)

(1+N)*(diff(u(y), y, y))+N*(diff(q(y), y)) = P

(1+N)*(diff(diff(u(y), y), y))+N*(diff(q(y), y)) = P

(3)

bc := u(-sigma) = 1, u(sigma) = -k

u(-sigma) = 1, u(sigma) = -k

(4)

``


 

Download help_couple_ode.mw

Hello

I want to know if there is any option to convert the Maple code into MATLAB version to work across platforms. I have tried the code generation option in Maple and I could not understand how to use the single line code generated. I am working on an already written Maple code and conversion of code to use on both platforms is kind of important. Please advice me on the same.

Regards!

I have a string which includes characters from the extended ASCII character set. I simply want to print the string. Nothing fancy. Hence my question is: how does one print characters from the extended ASCII character set so that they show up as the characters and not a bunch of empty boxes?  (A technique for UTF-8  character printing would be a bonus.)   

<Edit> For example, in integers, my string could be:

s := [56, 72, 157, 38];

Hence there printing of s as a string should be a 4 character / bytes entity: 

printf("%s", convert(s, bytes));

 

I work with polarization in optical fibers.  I want to (theoretically) scramble the polarization and cover the entire Poincare sphere with several thousand points using a sequence of two variable retarders.   So looking at the Mueller matrices in the attached image, as I vary the arguments of the sin & cosines, I’d like to see the polarization states traced out on the Poincare sphere.    Has this already been done in Maple, as in some kind of canned routine or add-on module?

Thanks!

Hi, 

While "numpoints=..." is said to be an admissible option of SurfaceOfRevolution (or VolumeOfRevolution, default value set to 50), the graphical results doesn't account for the value of numpoints.

Is this a known issue?

TIA

Surface0fRevolution.mw



 

Hello everyone, 

 

I am seeking for help!

You can freely access the mentioned file in any of the following links:

https://www.mdpi.com/2076-3417/9/15/2996

https://www.mdpi.com/2076-3417/9/15/2996/pdf

https://www.researchgate.net/publication/334711617_Transverse_Vibration_of_Clamped-Pinned-Free_Beam_with_Mass_at_Free_End

 

I have a/L a parameter. Give it a value and then beta can be calculated as in Table 1. The first five lowest physically possible values are listed only, but the number of beta-s tend to infinity because of the sine/cosine functions. (22) is clearly a nonlinear equation. As such, in general, it can't be solved analytically. (Of course, e.g., sin(x)=0 is also nonlinear but has analytical solution, but (22) is much more complicated.)
At first step, I give a/L  a value, 0.6 , then plot the left side of (22) and try to read the zeros. Here, beta on the horizontal axis is the independent variable. This method is really slow and inaccurate to get all the results of Figure 2. Instead, I should find the roots numerically by using a built-in solver in Maple.

Once I am done with the above, I should add a do/for loop for the parameter a/L. HOWEVER, I DO NOT KNOW HOW TO DO.

In conclusion, the main question is: How could I plot the graph using MAPLE in Figure 2 considering the transcendental Equation (22)?

Best Regards, 

Wallyson Thomas

restart;
with(GraphTheory);
with(SpecialGraphs);
with(StringTools);
GT := GraphTheory;
G := GT:-SpecialGraphs:-GridGraph(3, 3);
G := Graph(AdjacencyMatrix(G));
H := GT:-SpecialGraphs:-GridGraph(3, 2);
H := Graph(AdjacencyMatrix(H));
GH := CartesianProduct(G, H);
V := Vertices(GH);
P := [];
flag := 0
f := proc(K::string) local L; global P; L := Split(K, ":"); P := [op(P), cat("(", L[1], ",", L[2], ")")]; end proc
if flag = 0 then
    for i to NumberOfVertices(GH) do f(V[i]); end do;
end if
GH := RelabelVertices(GH, P)
DrawGraph(GH, stylesheet = "legacy")
 
I need the graph with those exact relabed vertices type i have done

Now the problem i want my graph a little better in output  i want my vertices a better visible the edges as they more not a issue atleast want the verices a little more visbile and distictive  kind help i apologize for the distubance kind help

the below is output i get i want it more better

Dear experts;

I would like to know the angle between lines P0-P1 and P0-P3. I believe I need to solve for all the angles to get this.

I have attempted to set this up in maple but... maple doesn't like what I've done. Can you tell me what I've done wrong?

geometry.mw

Hi,

It seems that plottools:-extrude doesn't support the "style" option: no error returned, just the extusion being always of surface style.
Am I correct ?

G := GridGraph(3, 3);
 DrawGraph(G);

Now in place  of (1,3) i want label 1 , in place  of (2,3) i want label 2, in place  of (3,3) i want label 3, in place  of (1,2) i want label 4, in place  of (2,2) i want label 5,in place  of (3,2) i want label 6,in place  of (1,1) i want label 7,in place  of (2,1) i want label 8,in place  of (3,1) i want label 9

And i want the graph to straight like the above and the edges should be correct and strainght kind help

 

Similary for below

 

G := GridGraph(2, 3);
 
DrawGraph(G);

 

 

Now in place  of (1,3) i want label 1 , in place  of (2,3) i want label 2, in place  of (1,2) i want label 3, in place  of (2,2) i want label 4, in place  of (1,1) i want label 5, in place  of (2,1) i want label 6

 And i want the graph to straight like the above and the edges should be correct and strainght kind help 

 

Kind help in such vertex labeled graph outputs and graphs should be strainght with edges and presented

If i can get both in legacy mode the above all in the above form in the legacy mode it will be more good as i want the color etc like that in my output exactly

 

kind help  i apologize for the disturbance kind help please

Vertices should not have have any of the (1,1),(2,1) etc in the output it should be labeled with the vertex labels as told in that order

And the edges straingth and good.

 

kind help  i apologize for the disturbance kind help please

First 403 404 405 406 407 408 409 Last Page 405 of 2097