Advertisement
Guest User

Untitled

a guest
Jun 26th, 2019
82
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.48 KB | None | 0 0
  1. ClearAll[y];
  2. DSolve[{y''[t] == -(1/(y[t] - 1)), y[0] == 0, y'[0] == 0}, y[t], t]
  3.  
  4. {y[t] -> E^-InverseErf[-Sqrt[(2/[Pi])] t]^2 (-1 + E^
  5. InverseErf[-Sqrt[(2/[Pi])] t]^2)}
  6.  
  7. y[t] := E^-InverseErf[-Sqrt[(2/[Pi])] t]^2 (-1 + E^
  8. InverseErf[-Sqrt[(2/[Pi])] t]^2);
  9. Plot[Evaluate[{y[t], 1/(y[t] - 1), D[y[t], t], D[y[t], {t, 2}]}], {t,
  10. 0, 2}]
  11.  
  12. Solve[y[t] == 1, t]
  13. {}
  14.  
  15. y[0.5]
  16. y[0.5]
  17.  
  18. Evaluate[y[0.5]]
  19. y[0.5]
  20.  
  21. y[t] /. t -> 0.5
  22. 0.12777
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement