Guest User

Untitled

a guest
Nov 14th, 2017
82
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.31 KB | None | 0 0
  1. Block[{x0, y0, eqns, soln},
  2. x0 = -1;
  3. y0 = 1;
  4. eqns = Join[
  5. Thread[{x'[s], y'[s]} == ({Re[#], Im[#]}/Abs[#])&[((x[s] - I y[s])^2)^(1/3)]],
  6. Thread[{x[0], y[0]} == {x0, y0}]
  7. ];
  8. soln = First@NDSolve[eqns, {x, y}, {s, 0, 3}];
  9.  
  10. ParametricPlot[
  11. {x[s], y[s]} /. soln
  12. , {s, 0, 3}
  13. , Frame -> True
  14. ]
  15. ]
Add Comment
Please, Sign In to add comment