>> From: C W "sylvester7"
Hi !
The following is an attempt to plot the lines crossig the point x =
21/50, y = 3/10 in the system [dx,dy] = [-1/20*x*(10*y-3),
-1/10*y*(10*y^2-7*y+10*x-3)]. How could I obtain the exact equations
and plots ?
S:=map(unapply,[-1/20*x*(10*y-3), -1/10*y*(10*y^2-7*y+10*x-3)],t,x,y);
solve(convert(S(0,x,y),set),{x,y});
with(share);
readshare(ODE,plots):
Digits:=16:smp_inits:=seq(evalf([0,cos(i*Pi*2/20)*0.01+0.42,
0.3+sin(i*Pi*2/20)*0.01]),i=0..20):
phaseplot(S,0..1,0..1,[25,25],{smp_inits},numsteps=600,
stepsize=0.05,iterations=5);
Chris
|
(NAG)
