Advanced Search

Maple Questions and Posts

These are Posts and Questions associated with the product, Maple
View unanswered posts
Maple Questions and Posts Feed

The problem of saving photos in eps format with a ...

Asked by: delvin 35
February 13 2023
0  4

 Hi,

The problem of saving photos in eps format with a small size
When I draw a function diagram with Maple, when I save the photo in eps format, the size is high, but I want the size to be high and the volume to be low, what should I do? If anyone can and knows how to draw with MATLAB, write the code for me, I would be grateful

Why a restart is needed to get a helpfull warning ...

Asked by: C_R 1960
integration int warning
February 12 2023
2  2

Update: partly resolved in Maple 2023. int_warning.mw still does not show the warning

There is a very helpful warning about using assumptions from the int command that does not appear in the attached case:
Warning, unable to determine if -1 is between 0 and x0; try to use assumptions or use the AllSolutions option
Warning, unable to determine if 1 is between 0 and x0; try to use assumptions or use the AllSolutions option

Using infolevel[int]:=5 in Maple 2022.2,  one can see that the integrals listed under int/ellalg/elltype are not the same.

My preference would be to have this warning always working. It's a reminder that Maple cannot give a general answer but that there might be solutions for a restricted real domain. The warning also does not reproduce when execution the corresponding help page
(Update: the warning reproduces in Maple 2023)

https://www.maplesoft.com/support/help/errors/view.aspx?path=Warning,%20unable%20to%20determine%20if%201%20is%20between%200%20and%20x__0;%20try%20to%20use%20assumptions%20or%20use%20the%20AllSolutions%20option

Sidenote: This kind of warning is so useful that I wonder if and how such warnings could be extended beyond  integration with piecewise solutions. (Maybe also to other commands. Maple often returns input unevaluated because there is no general solution in the complex domain, without informing the inexperienced user that there might be solutions available in the real domain.)

int_warning.mw

int_warning_with_infolevel.mw

Can 'Maple 2022' Simplify to a constant using L'ho...

Asked by: paulmcquad 60
limit simplification tutor
February 12 2023
1  2

My question is can the last step be equal to 1/6.

I want to output to 1/6. it outputs to sqrt(9)/18 - (See Below).

NULL"5)Square root: undefined - DNE (Limit does not exist):"

NULL

NULL

NULL

limit((sqrt(x+1)-3)/(x-8), x = 8)"(=)"1/6 

 

 

limit((sqrt(x+1)-3)/(x-8), x = 8)Limit(((x+1)^(1/2)-3)/(x-8), x = 8) = (1/18)*9^(1/2)"(=)"Limit(((x+1)^(1/2)-3)/(x-8), x = 8) = 1/6

NULL

Download limit-sqrt-5.mw

How can I search multiple worksheets for a specifi...

Asked by: Earl 930
February 12 2023
1  13

I would like to know which of the more than 100 worksheets I have in a Windows 11 folder execute the events operation within dsolve numeric.

Is there a way to identify which of this large number of Maple worksheets contain a given phrase e.g. "events"?

strange behavior of simplify. ...

Asked by: nm 8552
simplify
February 12 2023
3  1

What is the logic behind this?

expr:=(c[2]+x)^3;
simplify(expr)

Which is what expected. But 

expr:=(c[2]+x)^3+a;
simplify(expr)

gives

Which does not look simpler to me. I expected it to be the same as before but with "a" added.  This is what Mathematica gives for comparison

I know I can use simplify with size option. But my question is, how did Maple decide that x^3 + 3*x^2*c[2] + 3*x*c[2]^2 + c[2]^3 + a is "simpler" than (c[2] + x)^3 + a ? It must use some logic which I am trying to understand.

Maple 2022.2 on windows 10

Maple unter Window 11...

Asked by: panke 140
windows crash incomplete-question
February 12 2023
0  7

My Maple Program run perfectly  with Window 10.  My new Laptop has Window 11 and I get system  crash wenn I applied

plot3d(...) .

So I think It is somethink with the Grapics Redering ?

Can somebody help ?

Betriebsystemname    Microsoft Windows 11 Home
Version    10.0.22621 Build 22621
Weitere Betriebsystembeschreibung     Nicht verfügbar
Betriebsystemhersteller    Microsoft Corporation
Systemname    ELIPAN
Systemhersteller    SAMSUNG ELECTRONICS CO., LTD.
Systemmodell    750XED
Systemtyp    x64-basierter PC
System-SKU    SCAI-A5A5-A5A5-ADLP-PCFG
Prozessor    12th Gen Intel(R) Core(TM) i5-1235U, 1300 MHz, 10 Kern(e), 12 logische(r) Prozessor(en)
BIOS-Version/-Datum    American Megatrends International, LLC. P08CFG.033.220913.HQ, 13.09.2022
SMBIOS-Version    3.4
Version des eingebetteten Controllers    255.255
BIOS-Modus    UEFI
BaseBoard-Hersteller    SAMSUNG ELECTRONICS CO., LTD.
BaseBoard-Produkt    NP750XED-KC5DE
BaseBoard-Version    SAMSUNG_SW_REVISION_12345+10.0.22000
Plattformrolle    Mobil
Sicherer Startzustand    Ein
PCR7-Konfiguration    Erweiterung zum Anzeigen erforderlich
Windows-Verzeichnis    C:\WINDOWS
Systemverzeichnis    C:\WINDOWS\system32
Startgerät    \Device\HarddiskVolume1
Gebietsschema    Deutschland
Hardwareabstraktionsebene    Version = "10.0.22621.819"
Benutzername    ELIPAN\Pan
Zeitzone    Mitteleuropäische Zeit
Installierter physischer Speicher (RAM)    16,0 GB
Gesamter physischer Speicher    15,7 GB
Verfügbarer physischer Speicher    8,28 GB
Gesamter virtueller Speicher    16,7 GB
Verfügbarer virtueller Speicher    6,98 GB
Größe der Auslagerungsdatei    1,00 GB
Auslagerungsdatei    C:\pagefile.sys
Kernel-DMA-Schutz    Ein
Virtualisierungsbasierte Sicherheit    Wird ausgeführt...
Virtualisierungsbasierte Sicherheit – erforderliche Sicherheitseigenschaften    
Virtualisierungsbasierte Sicherheit – verfügbare Sicherheitseigenschaften    Allgemeine Virtualisierungsunterstützung, Sicherer Start, DMA-Schutz, UEFI-Code Readonly, SMM Security Mitigations 1.0, Modusbasierte Ausführungssteuerung, APIC-Virtualisierung
Virtualisierungsbasierte Sicherheit – konfigurierte Dienste    Durch Hypervisor erzwungene Codeintegrität
Virtualisierungsbasierte Sicherheit – ausgeführte Dienste    Durch Hypervisor erzwungene Codeintegrität
Windows Defender-Anwendungssteuerungsrichtlinie    Erzwungen
Windows Defender-Anwendungssteuerungs-Richtlinie für den Benutzermodus    Aus
Unterstützung der Geräteverschlüsselung    Erweiterung zum Anzeigen erforderlich
Es wurde ein Hypervisor erkannt. Features, die für Hyper-V erforderlich sind, werden nicht angezeigt.    

Save a plot as image directly via code...

Asked by: vs140580 490
February 12 2023
0  1

Save plot directly to a image file say as .jpg or .png 

Directly via maple code.

 

How do I create point structures in maple ...

Asked by: Liviu 5
February 11 2023
1  1

 The following code:

for i to np do
    pl_pp[i] := point([h_p[i], m_p[i]], color = red);
    pl_np[i] := point([h_n[i], m_n[i]], color = blue);
    pl_pn[i] := point([-h_n[i], -m_n[i]], color = red);
    pl_nn[i] := point([-h_p[i], -m_p[i]], color = blue);
end do

Gives following error:

Error, (in plottools:-point) incorrect arguments for creating points structure, try providing the dimension option
 

How to Solve Error, (in dsolve/numeric/bvp/convert...

Asked by: Saha 70
February 11 2023
0  1

Please Help,

I am trying to plot some movement equations, but the following error keeps happening:

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

when I tryed to rewrite the equation the error became:

Can someone help me with those errors?

Here is the code

pp.mw

I Encountered Error Message trying to use dsolve t...

Asked by: Elisha 35
ode parameters dsolve
February 11 2023
0  6

I find it difficult to use dsolve to solve system of ordinary differential equations with assigned parameters and initial conditions. The error message "Error, (in dsolve/numeric) 'parameters' must be specified as a list of unique unassigned names" kept coming up.

Pls see the uploaded equation for more understanding

> 

restart
> 

interface(imaginaryunit = F)

I

 (1)
> 

I

I

 (2)
> 

sqrt(-4)

2*I

 (3)
> 

NULL
> 

Suscep := diff(S(t), t) = theta*epsilon+v__2*S__v(t)-S(t)*lambda-S(t)*(µ+v__1)

diff(S(t), t) = theta*varepsilon+v__2*S__v(t)-S(t)*lambda-S(t)*(µ+v__1)

 (4)
> 

Vacc := diff(S__v(t), t) = (1-theta)*epsilon+v__1*S(t)-(µ+alpha+v__2)*S__v(t)-(1-w)*S__v(t)*lambda
> 

Immun := diff(V(t), t) = alpha*S__v(t)+`ρ__A`*A(t)+(1-k)*`ρ__Q`*Q(t)+`ρ__I`*(I)(t)-µ*V(t)
> 

Exp := diff(E(t), t) = S(t)*lambda+(1-w)*S__v(t)*lambda-(q__E+delta+µ)*E(t)
> 

Asymp := diff(A(t), t) = delta*a*E(t)-(`ρ__A`+µ)*A(t)+k*`ρ__Q`*Q(t)
> 

Inf := diff((I)(t), t) = delta*(1-a)*E(t)-(`ρ__I`+q__I+`δ__I`+µ)*(I)(t)
> 

Quar := diff((I)(t), t) = q__E*E(t)+q__I*(I)(t)-(`ρ__Q`+`δ__Q`+µ)*Q(t)
> 

init_conds := S(0) = S_0, S__v(0) = S__v*_0, V(0) = V_0, E(0) = E_0, A(0) = A_0, (I)(0) = I_0, Q(0) = Q_0

S(0) = S_0, S__v(0) = S__v*_0, V(0) = V_0, E(0) = E_0, A(0) = A_0, I(0) = I_0, Q(0) = Q_0

 (5)
> 

sys := {Asymp, Exp, Immun, Inf, Quar, Suscep, Vacc, init_conds}
> 

``
> 

sol := dsolve(sys, numeric, parameters = [`δ__Q`, `δ__I`, a, k, epsilon, v[1], q[E], q[I], q[A], eta[A], eta[Q], rho[A], rho[Q], rho[I], v[2], alpha, mu, delta, alpha, beta, w, lambda, S_0, S__v*_0, V_0, E_0, A_0, I_0, Q_0], method = rkf45)

Error, (in dsolve/numeric) 'parameters' must be specified as a list of unique unassigned names

  
> 

sol(parameters = [delta = .125, `δ__Q` = 0.6847e-3, epsilon = .464360344, `δ__I` = 0.2230e-8, a = .6255, q[E] = 0.18113e-3, k = .15, v__1 = 0.5e-1, v__2 = 0.6e-1, `ρ__Q` = 0.815e-1, `ρ__A` = .1, `ρ__I` = 0.666666e-1, q__I = 0.1923e-2, q__A = 0.4013e-7, `η__A` = .1213, `η__Q` = 0.3808e-2*alpha and 0.3808e-2*alpha = .4, w = .5925, mu = 0.464360344e-4, lambda = 0.1598643e-7, S_0 = 1.0, S__v*_0 = 0.6e-4, V_0 = 0.35e-4, E_0 = 0.5e-4, I_0 = 0.32e-4, A_0 = 0.15e-4, Q_0 = 0.1e-4])

sol(parameters = [delta = .125, delta__Q = 0.6847e-3, varepsilon = .464360344, delta__I = 0.2230e-8, a = .6255, q[E] = 0.18113e-3, k = .15, v__1 = 0.5e-1, v__2 = 0.6e-1, rho__Q = 0.815e-1, rho__A = .1, rho__I = 0.666666e-1, q__I = 0.1923e-2, q__A = 0.4013e-7, eta__A = .1213, false, w = .5925, mu = 0.464360344e-4, lambda = 0.1598643e-7, S_0 = 1.0, S__v*_0 = 0.6e-4, V_0 = 0.35e-4, E_0 = 0.5e-4, I_0 = 0.32e-4, A_0 = 0.15e-4, Q_0 = 0.1e-4])

 (6)
> 

Evaluate*the*system*at*t = 2
> 

sol(2)

sol(2)

 (7)
> 

sol(1)

sol(1)

 (8)
> 

sol(.1)

sol(.1)

 (9)
> 

sol(.3)

sol(.3)

 (10)
> 

sol(.5)

sol(.5)

 (11)
> 

sol(.7)

sol(.7)

 (12)
> 

sol(.9)

sol(.9)

 (13)
> 

sol(1.1)

sol(1.1)

 (14)
> 

sol(1.3)

sol(1.3)

 (15)
> 

sol(1.5)

sol(1.5)

 (16)
> 

 
> 

 

Download Covid19_Simulation.mw

Substituting indexes by an operator...

Asked by: Kmohseni 5
multiplication tensor
February 10 2023
1  4

Dear Maple community,

I am trying to define an operator in Maple that multiplies to an expression and substitutes the indices.

The issue is that I don't want to add the expression as input of the operator/function, sth like:

p := proc(i, j, z) option operator; subs({i = j, j = i}, z); end proc

I want it to act on expression, like:

P(i,j)*X(i,j,k)=X(j,i,k)

and X can be any tensor or matrix.

Thanks in advance.

Is it possible to print temporary output (and then...

Asked by: sursumCorda 922
output mathematica
February 10 2023
2  2

The shortcut key Ctrl＋Delete is used to clear a math region or delete a cell within a group, so I can delete some output manually. However, I find that there exists some long-standing functions like PrintTemporary and Monitor in Mathematica, which can be used to print intermediate values of variables temporarily (and dynamically, if possible). The term "temporarily" means that the printed output disappears at the end of a computation / when a computation is complete / when the evaluation finishes / before returning the result of an evaluation. For instance: 

(Note that certain operations require a front end.)

Is it possible to implement such functionalities in Maple? 

A problem in selecting elements from set...

Asked by: mmcdara 6149
diff set
February 10 2023
1  8

I get a problem while selecting elements from sets.

I have a list of two sets apparently ordered in the same order.
When I pick this, for instance, first set from each of these two samples I get the first one from the first set, and the second one from the second set.

> 

restart:
> 

DF := [ {diff(Theta(x), x), diff(Theta(x), x, x)}, {diff(g(x), x), diff(g(x), x, x)} ]

[{diff(Theta(x), x), diff(diff(Theta(x), x), x)}, {diff(diff(g(x), x), x), diff(g(x), x)}]

 (1)
> 

# Why are the elements of DF[2] not returned in their correct order?

i := 1;
DF[i], DF[i][1], DF[i][2];
print():
i := 2;
DF[i], DF[i][1], DF[i][2];

1

  

{diff(Theta(x), x), diff(diff(Theta(x), x), x)}, diff(Theta(x), x), diff(diff(Theta(x), x), x)

  

 

2

  

{diff(diff(g(x), x), x), diff(g(x), x)}, diff(diff(g(x), x), x), diff(g(x), x)

 (2)
> 

# I suppose we must not trust what is displayed and that lprint is more
# trustworthy?

lprint(DF)

[{diff(Theta(x), x), diff(diff(Theta(x), x), x)}, {diff(diff(g(x), x), x), diff(g(x), x)}]

  
> 

# I expected that converting each element of DF into a list would fix this
# problem. but anothr does appear:
# "Why does the conversion into list change the order of the elements?"
# I had never seen that before

LDF := convert~(DF, list)

[[diff(Theta(x), x), diff(diff(Theta(x), x), x)], [diff(diff(g(x), x), x), diff(g(x), x)]]

 (3)
> 

 

Download Selection_from_a_set.mw

Running this command

map(sort, LDF)

displays LDF, which suggests that LDF is already sorted according to some order.
What is this order which makes diff(Theta(x), x$2) the successor of diff(Theta(x), x)  but diff(g(x), x) the successor of diff(g(x), x$2)  ?

To find the table values of Nusselt and Skin frict...

Asked by: KIRAN SAJJAN 35
February 10 2023
0  2

odeSys := {diff(Theta(x), x, x)+Pr*(R*(diff(Theta(x), x))*f(x)+Nb*(diff(Theta(x), x))*(diff(Phi(x), x))+Nt*(diff(Theta(x), x))^2), N2*(diff(G(x), x, x))-N1*(2*G(x)+diff(f(x), x, x))-N3*R*((diff(f(x), x))*G(x)-f(x)*(diff(G(x), x))), diff(Phi(x), x, x)+R*Sc*f(x)*(diff(Phi(x), x))+Nt*(diff(Theta(x), x, x))/Nb, (1+N1)*(diff(g(x), x, x))+R*((diff(g(x), x))*f(x)-g(x)*(diff(f(x), x)))-M*g(x)+2*Kr*(diff(f(x), x)), (1+N1)*(diff(f(x), x, x, x, x))-R*((diff(f(x), x))*(diff(f(x), x, x))-f(x)*(diff(f(x), x, x, x)))+N1*(diff(G(x), x, x))-M*(diff(f(x), x, x))-2*Kr*(diff(g(x), x))}; cond := f(0) = 0, (D(f))(0) = 1, g(0) = 0, Theta(0) = 1, G(0) = -n*((D@@2)(f))(0), Phi(0) = 1, f(1) = lambda, (D(f))(1) = 0, g(1) = 0, Theta(1) = 0, G(1) = n*((D@@2)(f))(1), Phi(1) = 0; ans := {};

n := .5; N1 := 0.; N2 := 1.0; N3 := .1; lambda := .1; M := .1; Kr := .1; Sc := 1.0; Nb := .1; Pr := 1.0; Nt := .1; R := .5;

ans := dsolve*{cond, eval*odeSys};

hello these are the pde and Boundary conditions  i want to calculate the value of f''(0) ,Theta(0) and  Phi(0)

what is the proper cammand to get the table values for the given equation
NBVs := [eval(ans(N1*G(x)+(1+N1)*(diff(f(x), x, x))), x = 0), eval(ans(-(diff(Theta(x), x))), x = 0), eval(ans(-(diff(Phi(x), x))), x = 0)];

Numerical value for some quadratic equation....

Asked by: Al86 145
rootfinding fsolve numeric
February 09 2023
0  2

I have this piece of command in maple which I want to get a numerical value out of it, but instead of that it shows the sines without their numerical values. What to add to these commands to get a numerical value for x? (I checked maple's help for "numerical value", but perhaps I am a bit impatient and didn't find such a thing):

with(Degrees);
solve(x^2 = (4/sind(10))^2 + ((4 + x)/sind(20))^2 - 8*(4 + x)/(sind(10)*sind(20)), x);

I get the following:

                      1                        /  /   
---------------------------------------------- \8 \16 
                   2                                  
4 Degrees:-sind(10)  - 4 Degrees:-sind(10) + 1        

                   2                                        2   
  Degrees:-sind(10)  Degrees:-cosd(10) - 8 Degrees:-sind(10)  - 

  /                     4                  2
  \256 Degrees:-sind(10)  Degrees:-cosd(10) 

                          4                  
   - 192 Degrees:-sind(10)  Degrees:-cosd(10)

                         4
   + 32 Degrees:-sind(10) 

                         3                  
   - 64 Degrees:-sind(10)  Degrees:-cosd(10)

                         3
   + 24 Degrees:-sind(10) 

                         2                  
   + 16 Degrees:-sind(10)  Degrees:-cosd(10)

                        \      \\  
   - 2 Degrees:-sind(10)/^(1/2)//, 

                        1                        /  /   
  ---------------------------------------------- \8 \16 
                     2                                  
  4 Degrees:-sind(10)  - 4 Degrees:-sind(10) + 1        

                   2                                        2   
  Degrees:-sind(10)  Degrees:-cosd(10) - 8 Degrees:-sind(10)  + 

  /                     4                  2
  \256 Degrees:-sind(10)  Degrees:-cosd(10) 

                          4                  
   - 192 Degrees:-sind(10)  Degrees:-cosd(10)

                         4
   + 32 Degrees:-sind(10) 

                         3                  
   - 64 Degrees:-sind(10)  Degrees:-cosd(10)

                         3
   + 24 Degrees:-sind(10) 

                         2                  
   + 16 Degrees:-sind(10)  Degrees:-cosd(10)

                        \      \\
   - 2 Degrees:-sind(10)/^(1/2)//

First 111 112 113 114 115 116 117 Last Page 113 of 2097