# generate the RMS differences ;
# first invert;
mov i 1; mov k numatm; add k i;
loopi end: i k;
index j i; mul j.x -1;
end: ;
# then superimpose;
mov i 120;
add i 1; mul i 100;
use none tether;
gdock 0 10 100 1 10 0 i;
gdock 0 10 100 .1  1 0 i;
analyze 0 i;
use none bond angle hybrid torsion nonbon;