This script was working. For some reason I had to reinstall and not this
script no longer works.
load 'plot'
d=: 13 :'(0{x)-(-/x)*(i.>:y)%y' NB. domain intervals
T=:2p1 NB. Tau
rose=: 13 :'(1 o.x*y) (*"1) 2 1 o./y'
e=:(_1 1) d 5 NB. grid domain
INDEXES=:t=: ~.r=: <"1 |: b=: e I.(3 rose (0,T) d 9)*/0 1
INDEXES
plot (3 rose (0,T) d 5)*/0 1
plot ;/3 rose (0,T) d 385
plot ;/(3 rose (0,T) d 385)*/0 1
mask=: 13 :'1 INDEXES } 0 $~ 2 # >: y'
]shamrock=: ' X'{~ mask 25 NB. combine masks
Linda
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
