$#define FIGDIR ./Fig
$#define FIGLIST undamped damped subharm topharm nonlin nonlin2 nonlin3
#include <SEP.defs>

undamped.ascii: 
	duffing omega=1.0 epsilon=1.0 delta=0.0 gamma=0.0 y1=0.0 y2=1.0 pts=100 tmax=20.00
	cp phase undamped.ascii

FIGDIR/undamped.v: undamped.H 
	Graph <undamped.H title='Phase Space t=20' label1='u' label2='v' out=j1.v >/dev/null
	Real <undamped.H |Ftplot out=j2.v >/dev/null
	vp_OverUnderIso j1.v j2.v >FIGDIR/undamped.v
	
damped.ascii: 
	duffing omega=1.0 epsilon=1.0 delta=2.0 gamma=0.0 y1=0.0 y2=1.0 pts=100 tmax=20.00
	cp phase damped.ascii

FIGDIR/damped.v: damped.H
	Graph <damped.H title='Phase Space t=20' label1='u' label2='v' out=j1.v >/dev/null
	Real <damped.H |Ftplot out=j2.v >/dev/null
	vp_OverUnderIso j1.v j2.v >FIGDIR/damped.v


topharm.ascii: 
	duffing omega=4.67 epsilon=40.0 delta=0.0 gamma=10.0 y1=0.0 y2=1.0 pts=100 tmax=100.00
	cp phase topharm.ascii

FIGDIR/topharm.v: topharm.H 
	Graph <topharm.H title='Phase Space t=100' label1='u' label2='v' out=j1.v >/dev/null
	Real <topharm.H |Ftplot out=j2.v >/dev/null
	vp_OverUnderIso j1.v j2.v >FIGDIR/topharm.v

subharm.ascii: 
	duffing omega=4.67 epsilon=40.0 delta=0.0 gamma=0.67 y1=0.0 y2=1.0 pts=100 tmax=100.00
	cp phase subharm.ascii

FIGDIR/subharm.v: subharm.H 
	Graph <subharm.H title='Phase Space t=100' label1='u' label2='v' out=j1.v >/dev/null
	Real <subharm.H |Ftplot out=j2.v >/dev/null
	vp_OverUnderIso j1.v j2.v >FIGDIR/subharm.v
	
nonlin.ascii: 
	duffing omega=1 epsilon=10.0 delta=2.0 gamma=10 y1=0.0 y2=1.0 pts=100 tmax=100.00
	cp phase nonlin.ascii

FIGDIR/nonlin.v: nonlin.H 
	Graph <nonlin.H title='Phase Space' label1='u' label2='v' out=j1.v >/dev/null
	Real <nonlin.H |Ftplot out=j2.v >/dev/null
	vp_OverUnderIso j1.v j2.v >FIGDIR/nonlin.v
	
nonlin2.ascii: 
	duffing omega=10 epsilon=10.0 delta=2.0 gamma=10 y1=0.0 y2=1.0 pts=100 tmax=100.00
	cp phase nonlin2.ascii

FIGDIR/nonlin2.v: nonlin2.H 
	Graph <nonlin2.H title='Phase Space' label1='u' label2='v' out=j1.v >/dev/null
	Real <nonlin2.H |Ftplot out=j2.v >/dev/null
	vp_OverUnderIso j1.v j2.v >FIGDIR/nonlin2.v
	
nonlin3.ascii: 
	duffing omega=9.3 epsilon=40.0 delta=2.0 gamma=10 y1=0.0 y2=1.0 pts=100 tmax=100.00
	cp phase nonlin3.ascii

FIGDIR/nonlin3.v: nonlin3.H 
	Graph <nonlin3.H title='Phase Space' label1='u' label2='v' out=j1.v >/dev/null
	Real <nonlin3.H |Ftplot out=j2.v >/dev/null
	vp_OverUnderIso j1.v j2.v >FIGDIR/nonlin3.v
	

%.H: %.ascii
	atoF <%.ascii >%.H
	SETn1 %.H


clean: jclean
	-rm *.v *.ascii

$duffing:
$	make duffing

#include<SEP.idoc.rules>
#include<SEP.prog.rules>
#include<SEP.obj.rules>
