SHARE
TWEET

Untitled

a guest Nov 14th, 2017 54 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  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.  ]
RAW Paste Data
Top