# New system of two second order differential equations

I know one more thing from physical system. If we can assume the solutions in form $$x=Ae^{jp_1t}, \quad y=Be^{jp_2t}, \quad j=-1^{1/2}$$ I know that $$p_1=2p_2$$

If someone can help me. It is need to find a analytic or numeric solutions where $D_i$ are known constants. If the system can describe by lower number of constants and lower order, how can I get a numerical solutions in function of this constants using some of methods (perturbation or some software - Mathematica)

$$D_1x''+D_2y''(x'-y')-D_2x'y'+D_3x=0$$

$$D_4y''+D_2x''(x'-y')+D_2x'y'+D_5=0$$

Initial conditions

$$x(0)=a, \quad y(0)=0, \quad x'(0)=0, \quad y'(0)=0$$

where $(')=d/dt, ('')=d^2/dt^2$.

-
"numeric solutions" - in Mathematica, you can use NDSolve[]; it doesn't look to me like an obviously symbolically solvable DE system... –  Ｊ. Ｍ. May 7 '11 at 17:48