#!/usr/bin/gnuplot #!/usr/local/bin/gnuplot4 set title "RA/DE asteroidu z astorb.dat 15. 12. 2005; Hammerova stejnoplocha projekce" set parametric set trange [-180:180] set angles degrees set samples 1000 set size ratio -1 R = 1 X(lambda,phi) = 2.*R*sqrt(2.)*cos(phi)*sin(lambda/2.) / sqrt(1.+cos(phi)*cos(lambda/2.)) Y(lambda,phi) = R*sqrt(2.)*sin(phi) / sqrt(1.+cos(phi)*cos(lambda/2.)) set style line 1 lt 7 set style line 2 lt 2 lw 1 set style line 3 lt 3 pt 1 ps 0.5 set style line 4 lt 1 pt 7 ps 2 p "hammer.out" w d not ls 1,\ "Trojans.ham" w p t "Trojans" ls 3,\ "<./hammer.awk Jupiter.dat" w p t "Jupiter" ls 4,\ X(t,-60),Y(t,-60) not ls 2,\ X(t,-30),Y(t,-30) not ls 2,\ X(t, 0),Y(t, 0) not ls 2,\ X(t, 30),Y(t, 30) not ls 2,\ X(t, 60),Y(t, 60) not ls 2,\ X(-180, t/2),Y(-180, t/2) not ls 2,\ X(-150, t/2),Y(-150, t/2) not ls 2,\ X(-120, t/2),Y(-120, t/2) not ls 2,\ X( -90, t/2),Y( -90, t/2) not ls 2,\ X( -60, t/2),Y( -60, t/2) not ls 2,\ X( -30, t/2),Y( -30, t/2) not ls 2,\ X( 0, t/2),Y( 0, t/2) not ls 2,\ X( 30, t/2),Y( 30, t/2) not ls 2,\ X( 60, t/2),Y( 60, t/2) not ls 2,\ X( 90, t/2),Y( 90, t/2) not ls 2,\ X( 120, t/2),Y( 120, t/2) not ls 2,\ X( 150, t/2),Y( 150, t/2) not ls 2,\ X( 180, t/2),Y( 180, t/2) not ls 2 pa -1 call "png.plt" "Trojans.png"