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