; bdry1.dat
restore expand.sav
def bound
bp = ball_head
loop while bp # null
section
cp = b_clist(bp)
loop while cp # null
if c_nforce(cp) # 0.0 then
b2 = c_gobj2(cp)
if pointer_type( b2 ) = 101 then
b_xfix(bp) = 1
b_yfix(bp) = 1
b_zfix(bp) = 1
b_color(bp) = 1
exit section
end_if
end_if
if c_gobj1(cp) = bp
cp = c_go1clist(cp)
else
cp = c_go2clist(cp)
end_if
end_loop
end_section
bp = b_next(bp)
end_loop
end
bound
ini xvel = 0.0 grad=(-0.1,0,0)
ini yvel = 0.0 zvel = 0.0
prop xdisp=0.0 ydisp=0.0 zdisp=0.0
del wall 1
del wall 2
del wall 3
del wall 4
del wall 5
del wall 6
plot create assembly
plot add ball red green
plot add axes black
plot add vel black
plot show
pause
plot add cf white
cycle 200
plot disp black