Mathematica 9 is now available
Student Support Forum
-----
Student Support Forum: 'ODE' topicStudent Support Forum > General > Archives > "ODE"

Next Comment >Help | Reply To Topic
Author Comment/Response
Yan
06/28/11 1:26pm

Hye, I have an ODE and I can't solve it with DSolve or NDSolve.

EDO = y''[x] + ((1/
x) - (Y^2*(e^2*ne0*mu*
BesselJ[1,
mu*x/a]/(me*Eo*omega*
a*(omega -
I*v)))/((Ko^2*(1 - (e^2*ne0*
BesselJ[0, mu*x/a]/(me*Eo*omega*(omega - I*v)))) -
Y^2)*(1 - (e^2*ne0*
BesselJ[0, mu*x/a]/(me*Eo*omega*(omega - I*v)))))))*
y'[x] + (Ko^2*(1 - (e^2*ne0*
BesselJ[0, mu*x/a]/(me*Eo*omega*(omega - I*v)))) - Y^2)*y[x]

and:
Ed = 4.52
a = 0.013
b = 0.015
f = 200000000
Eo = 8.85418782*10^-12
muo = 1.25663706*10^-6
t = 0.1
s = 0.1
omega = f*2.*Pi
v = t*omega
omegap = omega/s
Ko = Sqrt[(omega^2)*Eo*muo]
mu = 1
me = 9.10938188*10^(-31)
e = 1.602176487*10^(-19)
If[mu > 0, ne0 = mu*omega^2*me*Eo/(2*s^2*e^2*BesselJ[1, mu]),
ne0 = omega^2*me*Eo/(s^2*e^2)]

I try to solve with: NDSolve[{EDO == 0, y[0] == 1, y'[0] == 0}, y[x], {x, 0.0001, 0.013}]

But it doesn't work.
I'm only interest in the value of y'[0.013] and y[0.013].

URL: ,

Subject (listing for 'ODE')
Author Date Posted
ODE Yan 06/28/11 1:26pm
Re: ODE Yan 06/29/11 1:23pm
Next Comment >Help | Reply To Topic