#!/usr/bin/gnuplot set xl "Da [AU]" set yl "De []" set zl "DI [deg]" set ticslevel 0 set size 0.666,1 set view 90,0,1,1 set tit "ERROR in calculation of r" sp '<./genvel_const 1e3 10 -1 | ./v_gauss 2.5 0.5 45 60' u 1:2:3 w p,\ '<./genvel_const 1e3 10 -1 | ./v_gauss_ERROR 2.5 0.5 45 60' u 1:2:3 w p pa -1 set tit "argument of perihelion" sp '<./genvel_const 1e3 10 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w d,\ '<./genvel_const 1e3 10 -1 | ./v_gauss 2.5 0 0 60' u 1:2:3 w d pa -1 set tit "eccentric anomaly" sp '<./genvel_const 1e3 10 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w d,\ '<./genvel_const 1e3 10 -1 | ./v_gauss 2.5 0 45 0' u 1:2:3 w d pa -1 set tit "isotropic power-law ejection, eccentric anomaly (E)" sp '<./genvel3 1e2 3.25 10 1e3 50 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w p,\ '<./genvel3 1e2 3.25 10 1e3 50 -1 | ./v_gauss 2.5 0 30 0' u 1:2:3 w p,\ '<./genvel3 1e2 3.25 10 1e3 50 -1 | ./v_gauss 2.5 0 60 0' u 1:2:3 w p pa -1 set tit "isotropic power-law ejection, vmin (R)" sp '<./genvel3 1e2 3.25 1 1e3 50 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w p,\ '<./genvel3 1e2 3.25 100 1e3 50 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w p,\ '<./genvel3 1e2 3.25 500 1e3 50 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w p pa -1 set tit "isotropic power-law ejection, parent body size (R)" sp '<./genvel3 1e2 3.25 1 1e3 10 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w p,\ '<./genvel3 1e2 3.25 1 1e3 100 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w p,\ '<./genvel3 1e2 3.25 1 1e3 1000 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w p pa -1 set tit "isotropic power-law ejection, slope (alpha)" sp '<./genvel3 1e3 1.0 1 1e3 50 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w p,\ '<./genvel3 1e3 2.0 1 1e3 50 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w p,\ '<./genvel3 1e3 5.0 1 1e3 50 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w p pa -1 set tit "conic ejection position" sp '<./genvel_cone 1e3 10 10 0 0 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w d,\ '<./genvel_cone 1e3 10 30 45 0 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w d,\ '<./genvel_cone 1e3 10 50 90 0 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w d pa -1 set tit "conic ejection angle" sp '<./genvel_cone 1e3 10 10 0 0 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w d,\ '<./genvel_cone 1e3 10 30 0 0 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w d,\ '<./genvel_cone 1e3 10 50 0 0 -1 | ./v_gauss 2.5 0 0 0' u 1:2:3 w d pa -1 set tit "velocity components (R,T,W)" sp '