Help required to plot the real part only in  multiple plots. i have written some codes  please help me to rectify the errors.Thanks in advance

restart:
with(plots):
n:=0.75:Eh:=100:mn:=1:t0:=0.2:
a1:=(mn/t0)^((n-1)/(n))*(tb)^(1/n):a3:=Eh/tb:a4:=(Eh)^2:
U1:=(a1/Eh)*(-1+(a3*r/Eh))^(1/n)*(n*a4/(1+n))*(1/a3-(r/Eh))-a1*(-1+(1/tb))^(1/n)*(n/(n+1)*(tb-1)):
plot([Re(seq(eval(U1,tb=j),j in[0.8,0.9,1.0]))],r=0..1,legend = [tb =0.8, tb=0.9,tb =1.0],  labels = ["z ", "U"], labeldirections = ["horizontal", "vertical"],  linestyle = [solid,dash,dot],color = [black, red,green]);

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

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. 
 

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 ?