Initial revision
[xgalaxy.git] / examples / three_generator.pl
1 use strict;
2 my $x = 0.97000436e0;
3 my $y = 0.24308753e0;
4 my $vx =0.93240737e0;
5 my $vy =0.86473146e0;
6
7 my $R=4.e7;
8 my $V=176.5e0;
9 my $M=7.4e22;
10
11 print<<EOT;
12 delta=3600
13 error=1E-10
14
15 EOT
16
17  
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";
21
22