MaplePrimes Questions

I am trying to find remainder of (2 + sqrt(3))^15 + (2 - sqrt(3))^15 with 2017. I tried
irem((2 + sqrt(3))^15 + (2 - sqrt(3))^15, 2017);

I do not get the result. How can I get the result?

I use like this 
a := expand((2 + sqrt(3))^15 + (2 - sqrt(3))^15);
irem(a, 2017);

get the result.

I am tyring to simulate the intensity autocorrelation of a femtosecond laser pulse.   The pulse intensity is set up as a 4096 point Vector of real values using a sech function for the pulse shape (see below). 

Being that it's the intensity, not the complex field, of the laser pulse, it is always real-valued and positive.    Therefore, the autocorrelation of the pulse intensity should be positive at all points in time.   But in fact, when I use the AutoCorrelation command on this Vector, it results many negative-valued points (see below).

What is going on?   Why am I getting negative values?

Thanks!

-Gregg-

I am new to Maple. I'm having trouble with simple caluclations such as the one in the image. The answer (after simplificatoin and cancellations) should be (|a|^2+|b|^2)(x^2 + y^2 + z^2)  but it's not happening! The result shown is correct; only the simplification is missing. What step am I missing?

The calculation is so simple, a look at the picture should suffice for the expert,

(sorry! I still don't know how to pick-off a single tab out from the whole workbook).

Question-2: judging from the missing overbar (complex conj) on z, the "assuming" seems to only apply to z and not x and y.

if that is correct, then why?

Question-3: How can I implement the constraint (|a|^2+|b|^2)=1 (preferably)before or after this  multiplication? I tried assume((|a|^2+|b|^2)=1), and though it compiled it did not help with the simplification.

TYVM

i have on ode with some parameters. it is ok and it is solved. but when i substitute the parameter itself, it is not solved, can i use any assumption to solve this? tnx for the help.

restart

ode := diff(T(x), x, x)+q/k = 0

diff(diff(T(x), x), x)+q/k = 0

(1)

dsolve(ode)

T(x) = -(1/2)*q*x^2/k+_C1*x+_C2

(2)

ics1 := -k*(D(T))(0) = h[1]*(T[inf1]-T(0)), -k*(D(T))(0.5e-1) = h[2]*(T(0.5e-1)-T[inf2])

-k*(D(T))(0) = h[1]*(T[inf1]-T(0)), -k*(D(T))(0.5e-1) = h[2]*(T(0.5e-1)-T[inf2])

(3)

dsolve({ics1, ode})

T(x) = -(1/2)*q*x^2/k-(1/40)*h[1]*(800*k*T[inf1]*h[2]-800*k*T[inf2]*h[2]-40*k*q-q*h[2])*x/((20*k*h[1]+20*k*h[2]+h[1]*h[2])*k)+(1/40)*(800*k*T[inf1]*h[1]+800*k*T[inf2]*h[2]+40*T[inf1]*h[1]*h[2]+40*k*q+q*h[2])/(20*k*h[1]+20*k*h[2]+h[1]*h[2])

(4)

ics2 := -k*(D(T))(0) = h[1]*(T[inf1]-T(0)), -k*(D(T))(L) = h[2]*(T(L)-T[inf2])

-k*(D(T))(0) = h[1]*(T[inf1]-T(0)), -k*(D(T))(L) = h[2]*(T(L)-T[inf2])

(5)

dsolve({ics2, ode})

Error, (in dsolve) found differentiated functions with same name but depending on different arguments in the given DE system: {T(L), T(x)}

 

``

Download ExactSol.mw

 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

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

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

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"?

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

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 plot directly to a image file say as .jpg or .png 

Directly via maple code.

 

 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
 

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 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

First 36 37 38 39 40 41 42 Last Page 38 of 2234