Guest User

Untitled

a guest
Jan 16th, 2018
79
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.89 KB | None | 0 0
  1. DSolve[(E^2 - L^2*((x - 1)^(2*d - 1)/((1 - y[x]^2)*(x + 1)^(2*d + 1))) - a*((x - 1)/(x + 1))^d)*Derivative[2][y][x] ==
  2. (E^2 - L^2*((x - 1)^(2*d - 1)/((1 - y[x]^2)*(x + 1)^(2*d + 1))) - a*((x - 1)/(x + 1))^d)*(((d^2 - 1)*y[x]*(1 - y[x]^2))/((x^2 - 1)*(x^2 - y[x]^2))) -
  3. L^2*((y[x]*(x - 1)^(2*d - 2))/((1 - y[x]^2)*(x + 1)^(2*d + 2))) +
  4. (a*((x - 1)^(d - 1)/(x + 1)^(d + 1)) - (E^2 - L^2*((x - 1)^(2*d - 1)/((1 - y[x]^2)*(x + 1)^(2*d + 1))) - a*((x - 1)/(x + 1))^d)*
  5. ((2*(x - d)*(x^2 - y[x]^2) + (d^2 - 1)*x*(1 - y[x]^2))/((x^2 - 1)*(x^2 - y[x]^2))) - L^2*(((x - 2*d)*(x - 1)^(2*d - 2))/((1 - y[x]^2)*(x + 1)^(2*d + 2))))*
  6. Derivative[1][y][x] - ((E^2 - L^2*((x - 1)^(2*d - 1)/((1 - y[x]^2)*(x + 1)^(2*d + 1))) - a*((x - 1)/(x + 1))^d)*
  7. (y[x]*((x^2 - y[x]^2 + (1 - d^2)*(1 - y[x]^2))/((1 - y[x]^2)*(x^2 - y[x]^2)))) + L^2*((y[x]*(x - 1)^(2*d - 1))/((1 - y[x]^2)^2*(x + 1)^(2*d + 1))))*
  8. Derivative[1][y][x]^2 + (((d*a)/(1 - y[x]^2))*((x - 1)/(x + 1))^d - (E^2 - L^2*((x - 1)^(2*d - 1)/((1 - y[x]^2)*(x + 1)^(2*d + 1))) - a*((x - 1)/(x + 1))^d)*
  9. ((x*(x^2 - 1) + d^2*x*(1 - y[x]^2) - 2*d*(x^2 - y[x]^2))/((1 - y[x]^2)*(x^2 - y[x]^2))) - L^2*(((x - 2*d)*(x - 1)^(2*d - 1))/((1 - y[x]^2)^2*(x + 1)^(2*d + 1))))*
  10. Derivative[1][y][x]^3, y[x], x]
  11.  
  12. Clear[y, d, L, a]
  13. y'[x] =.
  14. d = 1;
  15. L = 2;
  16. a = 3;
  17. NDSolve[{(E^2 -
  18. L^2*((x - 1)^(2*d - 1)/((1 - (y[x]^2))*(x + 1)^(2*d + 1))) -
  19. a*((x - 1)/(x + 1))^d)*
  20. y''[x] - (E^2 -
  21. L^2*((x - 1)^(2*d - 1)/((1 - (y[x]^2))*(x + 1)^(2*d + 1))) -
  22. a*((x - 1)/(x + 1))^d)*(((d^2 - 1)*
  23. y[x]*(1 - (y[x])^2))/((x^2 - 1)*(x^2 - (y[x]^2)))) -
  24. L^2*((y[x]*(x - 1)^(2*d - 2))/((1 - (y[x]^2))*(x + 1)^(2*d +
  25. 2))) + (a*(((x - 1)^(d - 1))/((x + 1)^(d + 1))) - (E^2 -
  26. L^2*((x - 1)^(2*d - 1)/((1 - (y[x]^2))*(x + 1)^(2*d + 1))) -
  27. a*((x - 1)/(x + 1))^
  28. d)*((2*(x - d)*(x^2 - (y[x]^2)) + (d^2 - 1)*
  29. x*(1 - (y[x]^2)))/((x^2 - 1)*(x^2 - (y[x]^2)))) -
  30. L^2*(((x -
  31. 2*d)*(x - 1)^(2*d - 2))/((1 - (y[x]^2))*(x + 1)^(2*d +
  32. 2))))*(y'[
  33. x]) - ((E^2 -
  34. L^2*((x - 1)^(2*d - 1)/((1 - (y[x]^2))*(x + 1)^(2*d + 1))) -
  35. a*((x - 1)/(x + 1))^d)*(y[
  36. x]*(x^2 - (y[x]^2) + (1 -
  37. d^2)*(1 - (y[x]^2)))/((1 - (y[x]^2))*(x^2 - (y[
  38. x]^2)))) +
  39. L^2*((y[x]*(x - 1)^(2*d - 1))/((1 - (y[x]^2))^2*(x + 1)^(2*d +
  40. 1))))*(y'[
  41. x])^2 + (((d*a)/(1 - (y[x]^2)))*((x - 1)/(x + 1))^
  42. d - (E^2 -
  43. L^2*((x - 1)^(2*d - 1)/((1 - (y[x]^2))*(x + 1)^(2*d + 1))) -
  44. a*((x - 1)/(x + 1))^d)*((x*(x^2 - 1) +
  45. d^2*x*(1 - y[x]^2) -
  46. 2*d*(x^2 - y[x]^2))/((1 - y[x]^2)*(x^2 - y[x]^2))) -
  47. L^2*(((x -
  48. 2*d)*(x - 1)^(2*d - 1))/((1 - y[x]^2)^2*(x + 1)^(2*d +
  49. 1))))*(y'[x])^3 == 0, y[0] == 2, y'[0] == 3},
  50. y[x], {x, 0, 30}]
Add Comment
Please, Sign In to add comment