MaplePrimes Questions

If I have a list say 

L:=[1,"2","3",4,5,"6",7,8,9,"10"]

Now i would like to have fuction which takes L as input say    split(L)

and would return two list one with characters from that namely ["2","3","6","10"] and another of integers [1,4,5,7,8,9] like this 

The list will have characters and integers only

Kind help it will be acknowleged

Q1: clicking on a variable in the attached I get

What does it mean and how can I fix that when it happens?

Initialisation_diagnostics_error.msim

Q2: Why are initial conditions which are set to "Ignore" listed as "Guess"?

I intend to use the initialization diagnostics to filter a model for components that have initial conditions set. In the attached example (excerpt from a larger model), this does not work the way I want it to (i.e. I don't want P1 to be listed a guess).

Init_diagnostics_ignore_listed_as_guess.msim

Q3: The app seems to be a powerfull tool. Is it explained somewhere how to use the app to debug a model?

Given a diagonal metric of signature (++- -)  for example the diagonal elements are [1,1,-1,-1], then how to find the Dirac gamma matrix in the standard representation or any representation.

restart;
with(LinearAlgebra);
with(plots);
NULL;

EqBIS := proc(P, U, V) local a, eq1, M1, t, PU, PV, bissec1; a := (P - U)/LinearAlgebra:-Norm(P - U, 2) + (P - V)/LinearAlgebra:-Norm(P - V, 2); M1 := P + a*t; eq1 := op(eliminate({x = M1[1], y = M1[2]}, t)); RETURN(op(eq1[2])); end proc;
NULL;
A := <4, 8>;
B := <4, 2>;
C := <1, 4>;
EqBIS(A, B, C);
NULL;

Cen := proc(M, N, R) local eq1, eq2, sol; eq1 := EqBIS(M, N, R) = 0; eq2 := EqBIS(N, M, R) = 0; sol := simplify(solve({eq1, eq2}, {x, y})); RETURN([subs(sol, x), subs(sol, y)]); end proc;
Cen(A,B,C):

CircleParm := t -> [(-t^2 + 1)/(t^2 + 1), 2*t/(t^2 + 1)];
P1 := Transpose(convert(CircleParm(1/4), Vector));

P2 := Transpose(convert(CircleParm(5), Vector));
P3 := Transpose(convert(CircleParm(-1/10), Vector));
P4 := Transpose(convert(CircleParm(-3/2), Vector));

C1 := Transpose(convert(Cen(P1, P2, P3), Vector));
Pts := [P1, P2, P3, P4, C1];

display(implicitplot([x^2 + y^2 - 1], x = -2 .. 2, y = -4 .. 2, colour = [blue], scaling = constrained), pointplot(Pts, symbolsize = 16));
Error, (in plots:-pointplot) incorrect number of coordinates in points data
NULL;Why this error. Thank you.

How relabel only one particular vertex of G of my choice. others will remain as it is

Function I want to create 

V:=Vertices(G):

Relabel:=proc(Graph::G,S,K)

I want to relabel only that particalur S subset V of G where S is list with say the K is a list of labels to be assigned the corrsponding vertices i paased in S. 

end proc;

Where v:=V[i] the vertex , k can be a character like       a or 1 or [1,2,3,4] or anything

To Maple support,. 

fyi,

internal error generated on this ode when trying odetest
 

interface(version);

`Standard Worksheet Interface, Maple 2022.1, Windows 10, May 26 2022 Build ID 1619613`

restart;

ode:=x^3*diff(y(x),x$3)+x^2*diff(y(x),x$2)-3*x*diff(y(x),x)+(x-1)*y(x)=0;
sol:=dsolve(ode,y(x),'series',x=0):
odetest(sol,ode,'series','point'=0);

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

Error, (in odetest/series) complex argument to max/min: 1/4800*((-79*I+3*83^(1/2))*3^(1/2)+9*I*83^(1/2)-79)*(316+12*83^(1/2)*3^(1/2))^(2/3)+1/12*I*(316+12*83^(1/2)*3^(1/2))^(1/3)*3^(1/2)-1/12*(316+12*83^(1/2)*3^(1/2))^(1/3)+20/3


 

Download problem_odetest_july_23_2022.mw

He Seniors, How to increase the quality of graph by puting command in ploting command according to required desire

Graph_qaulity_help.mw

Hi,

I found a bug in convert, please help me if you know a solution.

When we use the code

convert(U(xi)^2*y(xi)^5,list);

we get 

[U(xi)^2, y(xi)^5]

but when we use

convert(U(xi)^2,list);

we obtain

[U(xi),2]

while I expect the following result:

[U(xi)^2]

* "convert" is one part of my code for solving a problem. Sometimes there are cases that just exist one function instead of >1.

Since I want to keep my codes general, hence I need [U(xi)^2] in response to convert(U(xi)^2,list).

Thanks in advance.

A polynomial divided by another polynomial of lower degree. I want to find the remainder and quotient. How do I get.

The question is what should I add to the numerator get perfect division when x-3 x  -12 x+19 is divided by x2+x-6.

I shd get the answer as 2x+5 to be added since the remainder is -2x-5.   (Note that x-3 x  -10 x+24 is excatly divided by 2+x-6.

Can any one state the code to get the answer.

Quotient ( ) and Remainder (  ) codes , i tried, donot know how to use them.

Thanks for help. Cheers. RK


I'm using Maple 2021, and my Euler-Lagrange equation is coming out with a R[0] term that should not be there. I want to know what is causing this error and how I can fix it.

restartNULL

Lagrangian

 

Leq := (1/2)*R^2*[-cosh(rho1(tau)^2)*(diff(t(tau), tau))^2+(diff(rho1(tau), tau))^2+sinh(rho1(tau)^2)*(diff(theta1(tau), tau))^2]+(1/2)*R^2*[-cosh(rho2(tau)^2)*(diff(t(tau), tau))^2+(diff(rho2(tau), tau))^2+sinh(rho2(tau)^2)*(diff(theta2(tau), tau))^2]-(1/2)*k*rho1(tau)^2-(1/2)*k*rho2(tau)^2-(1/2*(tanh(rho2(tau)-rho1(tau))+1))*(rho2(tau)-rho1(tau))^2; L := subs({diff(rho1(tau), tau) = var4, diff(rho2(tau), tau) = var6, diff(t(tau), tau) = var2, diff(theta1(tau), tau) = var8, diff(theta2(tau), tau) = var10, rho1(tau) = var3, rho2(tau) = var5, t(tau) = var1, theta1(tau) = var7, theta2(tau) = var9}, Leq)

(1/2)*R^2*[-cosh(rho1(tau)^2)*(diff(t(tau), tau))^2+(diff(rho1(tau), tau))^2+sinh(rho1(tau)^2)*(diff(theta1(tau), tau))^2]+(1/2)*R^2*[-cosh(rho2(tau)^2)*(diff(t(tau), tau))^2+(diff(rho2(tau), tau))^2+sinh(rho2(tau)^2)*(diff(theta2(tau), tau))^2]-(1/2)*k*rho1(tau)^2-(1/2)*k*rho2(tau)^2-(1/2)*(-tanh(-rho2(tau)+rho1(tau))+1)*(rho2(tau)-rho1(tau))^2

(1.1)

Time Equation

 

epr11 := diff(L, var2); epr12 := diff(L, var1); epr13 := subs({var1 = t(tau), var10 = diff(theta2(tau), tau), var2 = diff(t(tau), tau), var3 = rho1(tau), var4 = diff(rho1(tau), tau), var5 = rho2(tau), var6 = diff(rho2(tau), tau), var7 = theta1(tau), var8 = diff(theta1(tau), tau), var9 = theta2(tau)}, epr11); epr14 := subs({var1 = t(tau), var10 = diff(theta2(tau), tau), var2 = diff(t(tau), tau), var3 = rho1(tau), var4 = diff(rho1(tau), tau), var5 = rho2(tau), var6 = diff(rho2(tau), tau), var7 = theta1(tau), var8 = diff(theta1(tau), tau), var9 = theta2(tau)}, epr12); epr15 := diff(epr13, tau); teq := epr15-epr14 = 0

2*R*[0]+(1/2)*R^2*[-4*rho1(tau)*(diff(rho1(tau), tau))*sinh(rho1(tau)^2)*(diff(t(tau), tau))-2*cosh(rho1(tau)^2)*(diff(diff(t(tau), tau), tau))]+(1/2)*R^2*[-4*rho2(tau)*(diff(rho2(tau), tau))*sinh(rho2(tau)^2)*(diff(t(tau), tau))-2*cosh(rho2(tau)^2)*(diff(diff(t(tau), tau), tau))] = 0

(2.1)

NULL


 

Download Maple_Help_Example.mw

 

Dear All,

I'm wondering if there is a simple way to define, given a "module" that protects "e" as a global, to allow the user to write e[0], e[1], e[2], e[3], e[inf] such that I calculate with e[4] but I display e[inf]. And if I type in e[inf], it translates behind the curtain of the module e[inf] into e[4] such that I keep my indexes 0-4 for the handling of my entities.

I hope I made myself understood ;)

I was woundering if there is some way with the standard "TypeSetting Module" ?

Sincerely

Ivar 

How can I achieve something like this (i.e. an output in a worksheet)?

It is important that the variables of the lables match the columns.

In the attached file is a typical dataset. A display of row lables would be optional.

Context: Easier interpretation (of the structure) of a Jacobian and subsequent manual solving the equations of interest that correspond to the rows.

 

Labels_for_sparse.mw

Hi!

I would like Maple to show the full solution method to this pde.  How do I do this?

 

Show_Solution_to_NLPDE.mw

Hello my friends
Please help me solve the integration problem in the attached Maple code below (why was the integration not done in u[4]?) ... Thank you very much

``

u[1] :=  1/3*t^3*x + 2*t*x + (-1)*0.7522527780*t^(3/2)*x + (-1)*0.1719434922*t^(7/2)*x + (-1)*0.01587301587*t^7*x + (-1)*0.1333333333*t^5*x

(1/3)*t^3*x+2*t*x-.7522527780*t^(3/2)*x-.1719434922*t^(7/2)*x-0.1587301587e-1*t^7*x-.1333333333*t^5*x

(1)

u[2] := subs(t=Tau,u[1]);

(1/3)*Tau^3*x+2*Tau*x-.7522527780*Tau^(3/2)*x-.1719434922*Tau^(7/2)*x-0.1587301587e-1*Tau^7*x-.1333333333*Tau^5*x

(2)

u[3]:= int(u[2],Tau = 0 .. t);

t^2*x+0.8333333333e-1*t^4*x-0.2222222222e-1*t^6*x-0.1984126984e-2*t^8*x-.3009011112*t^(5/2)*x-0.3820966493e-1*x*t^(9/2)

(3)

u[4]:= int((t-Tau)^(-0.5)*u[2],Tau = 0 .. t);

int(((1/3)*Tau^3*x+2*Tau*x-.7522527780*Tau^(3/2)*x-.1719434922*Tau^(7/2)*x-0.1587301587e-1*Tau^7*x-.1333333333*Tau^5*x)/(t-Tau)^.5, Tau = 0 .. t)

(4)

``

Download help.mw

When displaying two tubeplots together, we may specify their colors at will, as long as they are different colors!  For instance, specifying red and green works correctly, but specifying red and red results in red and black!

See the attached worksheet.  Interestingly, when displaying the contents of the worksheet on this website, the colors are rendered correctly!  So don't go with what you see on this web page; look inside the worksheet instead.

restart;

kernelopts(version);

`Maple 2022.1, X86 64 LINUX, May 26 2022, Build ID 1619613`

with(plots):

Two intersecting tori colored red and green -- works as expected:

display(
        tubeplot([cos(t), 0, sin(t)], t=-Pi..Pi, radius=0.2),
        tubeplot([cos(t), sin(t), 0], t=-Pi..Pi, radius=0.2),
style=surface, color=[red,green]);

When we set both colors to red, one of the surfaces is painted black!  Why?

Please note: This website displays the colors corectly as red and red.  But

within the worksheet the colors are read and black.

display(
        tubeplot([cos(t), 0, sin(t)], t=-Pi..Pi, radius=0.2),
        tubeplot([cos(t), sin(t), 0], t=-Pi..Pi, radius=0.2),
style=surface, color=[red,red]);

Specifying colors as red/red within the tubeplots still produces red/black!

display(
        tubeplot([cos(t), 0, sin(t)], t=-Pi..Pi, radius=0.2, color=red),
        tubeplot([cos(t), sin(t), 0], t=-Pi..Pi, radius=0.2, color=red),
style=surface);

Download mw.mw

PS: As a workaround, we may replace the red & red specification with
COLOR(RGB, 1, 0, 0) and
COLOR(RGB, 1, 0, 0.01)
which are different enough to make Maple happy, but produce essentially the same red color.

First 173 174 175 176 177 178 179 Last Page 175 of 2308