# ifs.ode
# simple iterated function system
# sierpinsky triangle
par c0=0,c1=.5,c2=0
par d0=0,d1=0,d2=.5
p=flr(ran(1)*3)
x'=.5*x+shift(c0,p)
y'=.5*y+shift(d0,p)
@ meth=discrete,total=20000,maxstor=100000
@ xp=x,yp=y,xlo=0,xhi=1,ylo=0,yhi=1,lt=0
aux pp=p
done