Maple 13 Questions and Posts

how to get the output of the plots in rows and col...

Asked by: imparter 175

July 14 2018

0 3

any one can help me to get the output of the plot in 3 rows and 2 columns here is my codes. thanks in advance.

restart;
h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/Gamma2)-h(z)^2+Nb*h(z)^4):
lambda:=unapply(Int(K1,z=0..1), Gamma2):
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
with(plots):
display
( Vector[row]
( [ seq
( plot
( [ seq
( eval(lambda(Gamma2), Nb=j),
j in [0.1,0.2,0.3]
)
],
delta2=0.02..0.1,
legend=[Nb=0.1,Nb=0.2,Nb=0.3],
labels=[typeset(`δ1`), typeset(conjugate(`Δp`))],
title=typeset("Effect of ", ''alpha'', " when ", Gamma,"2=", Gamma2)
),
Gamma2 in [10,20,30,40,50,-10]
)
]
)
);

how to remove error in plot display...

Asked by: imparter 175

July 08 2018

0 4

Hellow ,

any body can help me to remove the error from the plot, here is my codes .

Note. 1. *Title {Effect of alpha when 'Gamma2'=Gamma2}

2. *on the y axis i want bar(deltap)

restart;

h:=z->1-(delta2/2)*(1 + cos(2*(Pi/L1)*(z - d1 - L1))):
K1:=((4/h(z)^4)-(sin(alpha)/Gamma2)-h(z)^2+Nb*h(z)^4):
lambda:=unapply(Int(K1,z=0..1), Gamma2):
L1:=0.2:
d1:=0.2:
alpha:=Pi/6:
plots:-display(Vector[row]([seq(plot([seq(eval(lambda(Gamma2), Nb=j), j=[0.1,0.2,0.3])], delta2=0.02..0.1, legend=[Nb=0.1,Nb=0.2,Nb=0.3],labels = ["d1 ", "Dp"],labelfont = ["SYMBOL", 12], title=(Effect of alpha when 'Gamma2'=Gamma2)), Gamma2=[10,20,30])]),axes=boxed);

height and width of plot in maple 13...

Asked by: imparter 175

July 05 2018

0 3

how to adjust the height and width of plot in maple 13

plot( x[3]^5, caption = typeset("\n A plot of %1.", x[3]^5), captionfont=[times, 20] );

caption size and font in plot...

Asked by: imparter 175

July 05 2018

0 3

hellow,

can any body help me to increase the size of the caption in maple plot

plot( x[3]^5, caption = typeset("A plot of %1.", x[3]^5) );

dslove trigger time with rhs(SOL[1](last)...

Asked by: panke 145

June 29 2018

0 6

I want to fetch the trigger time.
I use rhs function.

>

restart; ode := diff(y(x), x) = 1/(y(x)+1); event := [[abs(y(x)+1)-0.1e-1, halt]]

>

Warning, cannot evaluate the solution further left of -.49995006, event #1 triggered a halt

>

(1)

>

(2)

>

Download RHS_TEST.mw

dsolve numeric error ...

Asked by: panke 145

June 28 2018

1 0

>

`odeλ` := diff(lambda(tau), tau, tau)+2*a*Ep(tau)*(diff(lambda(tau), tau))/(Xi(tau)*a+1) = 0:

>

th := algsubs((diff(lambda(tau), tau))^2 = (-1-(diff(xi(tau), tau))^2)/(-(xi(tau)*a+1)^2), `odeξ`):

a*(diff(xi(tau), tau))^2/(xi(tau)*a+1)+diff(diff(xi(tau), tau), tau)+a/(xi(tau)*a+1) = 0

a*(diff(xi(tau), tau))^2/(xi(tau)*a+1)+diff(diff(xi(tau), tau), tau)-a/(xi(tau)*a+1) = 0

(1)

>

ds2 := -(xi(tau)*a+1)^2*(diff(lambda(tau), tau))^2+(diff(xi(tau), tau))^2:

>

`λp1` := 1:

>

`icsλ1` := lambda(0) = 0, (D(lambda))(0) = `λp1`:

(2)

>

`lpξ1` := dsolve([`odeξ1`, `icsξ1`], numeric, output = listprocedure, range = teu .. te, events = event1); Xi := eval(xi(tau), `lpξ1`); `Ξp` := eval(diff(xi(tau), tau), `lpξ1`); `odeλ`; `lpλ1` := dsolve([`odeλ`, `icsλ1`], numeric, output = listprocedure, range = teu .. te, events = event1)

diff(diff(lambda(tau), tau), tau)+.6*Ep(tau)*(diff(lambda(tau), tau))/(.3*Xi(tau)+1) = 0

Error, (in dsolve/numeric/DAE/make_proc) number of unknown functions and equations must match, got 3 functions {Ep, Xi, lambda}, and 1 equations

>

(3)

>

`λep1` := eval(diff(lambda(tau), tau), lp1); `ξep1` := eval(diff(xi(tau), tau), lp1); `λe1` := eval(lambda(tau), lp1); `ξe1` := eval(xi(tau), lp1)

(4)

>

`λep2` := eval(diff(lambda(tau), tau), lp2); `ξep2` := eval(diff(xi(tau), tau), lp2); `λe2` := eval(lambda(tau), lp2); `ξe2` := eval(xi(tau), lp2)

(5)

>

ds2h1 := subs([lambda(tau) = `λe1`(tau), xi(tau) = `ξe1`(tau), diff(lambda(tau), tau) = `λep1`(tau), diff(xi(tau), tau) = `ξep1`(tau)], ds2); dse1 := unapply(ds2h1, tau); ds2h2 := subs([lambda(tau) = `λe2`(tau), xi(tau) = `ξe2`(tau), diff(lambda(tau), tau) = `λep2`(tau), diff(xi(tau), tau) = `ξep2`(tau)], ds2); dse2 := unapply(ds2h2, tau)

>

(6)

>

t := proc (xi, lambda) options operator, arrow; (xi+1/a)*sinh(a*lambda) end proc; x := proc (xi, lambda) options operator, arrow; (xi+1/a)*cosh(a*lambda) end proc; xi0 := 1; teg := 7

(7)

>

p1 := plot([[x(xi0, `λe1`(tau)), t(xi0, `λe1`(tau)), tau = teu .. te]], legend = ['ds' = -1], color = ["red"]); p2 := plot([[t(xi0, `λe2`(tau)), x(xi0, `λe2`(tau)), tau = teu .. te]], legend = ['ds' = 1], color = ["blue"]); p3 := plot([[tau, tau, tau = -teg .. teg], [tau, -tau, tau = -teg .. teg]], color = ["black", "black"], legend = ['ds' = 0, 'ds' = 0])

>

Download Rindler_simulation_v2.mw

extracting coefficients of exponential functions...

Asked by: guras 100

June 20 2018

1 2

Consider for instance the following equation:

Eq:=(a-4)*exp(4*x)+(b+1)*exp(2*x)+(c-2)=0

How can I list the coefficients of the exponential functions and also solve the equation for the constant parameters

a, b, and c?

I tried

[coeffs(collect(lhs(EQ), exp), exp)] =~ 0;

but it did not work. Thank you for your help.

fetch the dsolve event time ...

Asked by: panke 145

June 19 2018

1 1

1) How can i get the event time. This example 0.69313.

2) I don't want to see the Warning.How can i do it?

Warning, cannot evaluate the solution further right of .69314729, event #1 triggered a halt

dsolve with event...

Asked by: panke 145

June 12 2018

0 10

Why this event dose'nt work ?
S(t) is a state , a parameter
dsolve(...,numeric,events = [[[s(t), a*arcsinh(2/a) < s(t)], halt]]);

Having problems with expansion series (mtaylor)...

Asked by: Mark1993 0

June 11 2018

0 5

Good evening!!!

Let me briefly describe the problem I've faced recently.

The program (attached) deals with a rather complicated function f depending on parametrs eps1, eps2, eps3, eps4 and variable w. The aim is to expand the function f(w1) into Taylor series with respect to all parametrs (eps1, eps2, eps3, eps4) in order to study its asymptotic behavior as function depending only on k; 0<k<1.

I decided to use mtaylor-function for that problem, which (as I've understood) is the only one to be applied in such cases, but the result was rather unsatisfactory, an error:

Error, (in gcd/LinZip) input must be polynomials over the integers

Programm code: (1)-(12) only announcing functions....(((, see below

>

(1)

>

(2)

>

(1/3)*(B1-D1+(B1^2-2*B1*D1+D1^2+3*A1^2-3*A1*C1)^(1/2))/A1

(3)

>

-B1+(1/3)*(A1-C1)*(B1-D1+(B1^2-2*B1*D1+D1^2+3*A1^2-3*A1*C1)^(1/2))/A1+(1/9)*(B1-D1)*(B1-D1+(B1^2-2*B1*D1+D1^2+3*A1^2-3*A1*C1)^(1/2))^2/A1^2-(1/27)*(B1-D1+(B1^2-2*B1*D1+D1^2+3*A1^2-3*A1*C1)^(1/2))^3/A1^2

(4)

>

s := eps4*sin(l*tau)+(4*(l*sqrt(k/(1-k))+l*eps3)+2*l*((1-2*k)/sqrt(k*(1-k))+eps1))/l^2

eps4*sin(l*tau)+(4*l*(k/(1-k))^(1/2)+4*l*eps3+2*l*((1-2*k)/(k*(1-k))^(1/2)+eps1))/l^2

(5)

>

A1 := (2*(l*sqrt(k/(1-k))+l*eps3)+l*((1-2*k)/sqrt(k*(1-k))+eps1))/s

(2*l*(k/(1-k))^(1/2)+2*l*eps3+l*((1-2*k)/(k*(1-k))^(1/2)+eps1))/(eps4*sin(l*tau)+(4*l*(k/(1-k))^(1/2)+4*l*eps3+2*l*((1-2*k)/(k*(1-k))^(1/2)+eps1))/l^2)

(6)

>

A1 := (2*(l*sqrt(k/(1-k))+l*eps3)+l*((1-2*k)/sqrt(k*(1-k))+eps1))/s

(2*l*(k/(1-k))^(1/2)+2*l*eps3+l*((1-2*k)/(k*(1-k))^(1/2)+eps1))/(eps4*sin(l*tau)+(4*l*(k/(1-k))^(1/2)+4*l*eps3+2*l*((1-2*k)/(k*(1-k))^(1/2)+eps1))/l^2)

(7)

>

4/(eps4*sin(l*tau)+(4*l*(k/(1-k))^(1/2)+4*l*eps3+2*l*((1-2*k)/(k*(1-k))^(1/2)+eps1))/l^2)^2

(8)

>

C1 := (((1-2*k)/sqrt(k*(1-k))+eps1)^2+(-(1-2*k)/sqrt(k*(1-k))+eps2)^2)/s^2

(((1-2*k)/(k*(1-k))^(1/2)+eps1)^2+(-(1-2*k)/(k*(1-k))^(1/2)+eps2)^2)/(eps4*sin(l*tau)+(4*l*(k/(1-k))^(1/2)+4*l*eps3+2*l*((1-2*k)/(k*(1-k))^(1/2)+eps1))/l^2)^2

(9)

>

D1 := (2*((1-2*k)/sqrt(k*(1-k))+eps1))*(-(1-2*k)/sqrt(k*(1-k))+eps2)/s^2

2*((1-2*k)/(k*(1-k))^(1/2)+eps1)*(-(1-2*k)/(k*(1-k))^(1/2)+eps2)/(eps4*sin(l*tau)+(4*l*(k/(1-k))^(1/2)+4*l*eps3+2*l*((1-2*k)/(k*(1-k))^(1/2)+eps1))/l^2)^2