fix marbletrjbin.m
output rotation matrix in superimpose.m
add calcrg.m (calc radius of gyration)