Guest User

Untitled

a guest
Jun 21st, 2018
69
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.36 KB | None | 0 0
  1. Solve[(47^(u + 1) + 3^(u + 1))/(u + 1) ==(50^1.5)/1.5,u]
  2.  
  3. Plot[(47^(u + 1) + 3^(u + 1))/(u + 1), {u, -50, 50},
  4. PlotRange -> {-1, 236}]
  5.  
  6. Solve[(47^(u+1)+3^(u+1))/(u+1)==(50^(3/2))/(3/2) && Abs[u]<1,u, Reals]
  7.  
  8. (* {{u->Root[{500 Sqrt[2]-3^(2+#1)-3 47^(1+#1)+500 Sqrt[2] #1&,-0.99133010386676537502}]},{u->Root[{500 Sqrt[2]-3^(2+#1)-3 47^(1+#1)+500 Sqrt[2] #1&,0.52441226301643480967}]}} *)
  9.  
  10. N[%]
  11.  
  12. (* {{u->-0.99133},{u->0.524412}} *)
  13.  
  14. Clear[f]
  15. f[u_] := (47^(u + 1) + 3^(u + 1))/(u + 1) - (50^1.5)/1.5
  16.  
  17. sol =
  18. Reap@
  19. NDSolve[{
  20. D[y[u], u] == D[f[u], u], y[0] == f[0],
  21. WhenEvent[y[u] == 0, Sow[u]]},
  22. y, {u, -2, 2}
  23. ];
  24.  
  25. sol[[2, 1]]
  26. (* Out: {-0.99133, 0.524412} *)
  27.  
  28. Show[
  29. Plot[f[u], {u, -2, 2}, PlotRange -> 1000],
  30. ListPlot[
  31. Callout[
  32. {#, 0}, Round[#, 0.001],
  33. LabelStyle -> Directive[Red, Medium]
  34. ] & /@ sol[[2, 1]],
  35. PlotStyle -> {Red, PointSize[0.015]}
  36. ]
  37. ]
  38.  
  39. With[
  40. {
  41. lhseqn = (47^(u + 1) + 3^(u + 1))/(u + 1),
  42. rhseqn = (50^1.5)/1.5
  43. },
  44. {u, lhseqn} /. Solve@Reduce[
  45. lhseqn == rhseqn
  46. , u
  47. , Reals
  48. ]
  49. ]
  50.  
  51. (* {{-0.99133, 235.702}, {0.524412, 235.702}} *)
  52.  
  53. Plot[
  54. {(3^(1 + u) + 47^(1 + u))/(1 + u), 235.70226039551585`}
  55. , {u, -2, 1}
  56. , Epilog -> {
  57. Red,
  58. PointSize[Large],
  59. Point[
  60. {
  61. {-0.9913301038667653`, 235.702260395515`},
  62. {0.5244122630164348`, 235.7022603955158`}
  63. }
  64. ]
  65. }]
Add Comment
Please, Sign In to add comment