18 print $M,"\t",sprintf("%.13E",-$x*$R),"\t",sprintf("%.13E", $y*$R),"\t0\t",sprintf("%.13E", $vx*$V),"\t",sprintf("%.13E", $vy*$V), "\t0\n";
19 print $M,"\t", 0,"\t", 0,"\t0\t",sprintf("%.13E", -2*$vx*$V),"\t",sprintf("%.13E", -2*$vy*$V), "\t0\n";
20 print $M,"\t",sprintf("%.13E", $x*$R),"\t",sprintf("%.13E", -$y*$R),"\t0\t",sprintf("%.13E", $vx*$V),"\t",sprintf("%.13E", $vy*$V), "\t0\n";