saysics

10 Reputation

2 Badges

11 years, 207 days

MaplePrimes Activity


These are questions asked by saysics

Hell

I write this code and didn't work , I have some erorrs as

Warning, The use of global variables in numerical ODE problems is deprecated, and will be removed in a future release. Use the 'parameters' argument instead (see ?dsolve,numeric,parameters)
Warning, cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

also I have question " How I can change the scale of plot"

parameters := [z = 0, Omega = 2.2758, tau = 13.8, T2 = 200, omega0 = 1, r = .7071, s = 2.2758, H = 1.05457173*10^(-34), omega = .5, k = 1666666.667, Delta = 1.7758]

 

sys1 := {diff(u(t), t) = s*v(t)-u(t)/T2, diff(v(t), t) = -s*u(t)-2*Omega*exp(-r^2/omega0^2-t^2*1.177^2/tau^2)*cos(kz-`ωt`)*w(t)-v(t)/T2, diff(w(t), t) = 2*Omega*exp(-r^2/omega0^2-t^2*1.177^2/tau^2)*cos(kz-`ωt`)*v(t)}; ICs1 := {u(-20) = 0, v(-20) = 0, w(-20) = -1}

 

ans1 := dsolve(`union`(eval(sys1, parameters), ICs1), numeric, output = listprocedure); plots:-odeplot(ans1, [[t, u(t)], [t, v(t)], [t, w(t)]], t = -20 .. 20, legend = [w, v, u])

 

U := eval(u(t), ans1); F := eval(((-2*10^33*Omega*H*r*U(t))*(1/omega0^2))*exp(-r^2/omega0^2-t^2*1.177^2/tau^2)*cos(kz-`ωt`), parameters)

 

plot(F, t = -20 .. 20)

I write this system but I have 2 error

 

restart; params := [z = 0,

Omega = 2.2758,

tau = 13.8, T2 = 200,

omega0 = 1,

r = .7071,

s = 2.2758,

omega = .5]

 

sys1 := {diff(q(t), t) = -2*Omega*v(t)-s*exp(-r^2/omega0^2-t^2*1.177^2/tau^2)*cos(k*z-omega*t)*(y(t)-x(t))-q(t)/T2,

diff(v(t), t) = Omega*q(t)-v(t)/T2,

diff(x(t), t) = 2*s*exp(-r^2/omega0^2-t^2*1.177^2/tau^2)*cos(k*z-omega*t)*q(t)+y(t)/T1,

diff(y(t), t) = -2*s*exp(-r^2/omega0^2-t^2*1.177^2/tau^2)*cos(k*z-omega*t)*q(t)-y(t)/T1};

ICs1 := {q(-20) = 0, v(-20) = 0, x(-20) = 1, y(-20) = 0}

 

 

ans1 := dsolve(`union`(eval(sys1, params), ICs1), numeric, output = listprocedure); plots:-odeplot(ans1, [[t, x(t)], [t, y(t)], [t, q(t)], [t, v(t)]], t = -20 .. 20, legend = [x, y, q, v])

 

Error, invalid input: eval received params, which is not valid for its 2nd argument, eqns
Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

Hello everybody

I'm new at using Maple

so what I'm trying to do is " solve system of differential equations numerically " and plot the result 

I use the floweing code

 

PDEtools[declare]((u, v, w)(t), prime = t)

> params := z = 0;

Omega= 2.2758;

tau = 13.8;

T2 = 200; s = 1;

r = 0.7071;

\[CapitalDelta] = 1.7758;

s = 2.2758;

Eta= 1.05457173*10^-34;

omega = 0.5; k = 1666666.667;

> sys1 := {diff(u(t), t) = Omega*v(t)-u(t)/T2,

diff(v(t), t) = -Omega*u*{t}-2*s*exp(-r^2/omega0^2-t^2*1.177^2/tau^2)*cos(k*z-omega*t)*w(t)-v(t)/T2,

diff(w(t), t) = 2*s*exp(-r^2/omega0^2-t^2*1.177^2/tau^2)*cos(k*z-omega*t)*v(t)};

Cs1 := {u(-20) = 0, v(-20) = 0, w(-20) = -1}

> ans1 := dsolve*RealRange(Open({ICs1, sys1}), {u(t), v(t), w(t)});
%;
Error, (in RealRange) invalid arguments

plot([u(t),t=-20..20])
plot([v(t),t=-20..20])
plot([w(t),t=-20..20])

 

 

:::::::::

also I need to use the result of v(t) in another equation as,

x=2*v(t)*cos(k*z-omega*t)

How I can do that ?

 

Page 1 of 1