load(draw);
draw2d(terminal=aquaterm,points_joined=true,
    grid=true,line_width=1.5,
    xlabel="Planet",ylabel="Mean distance from Sun (AU)",
    xtics={["Mercury",1],["Venus",2],["Earth",3],["Mars",4],["(Ceres)",5],
        ["Jupiter",6],["Saturn",7],["Uranus",8],["Neptune",9],["(Pluto)",10]},
    key="T-B Rule prediction",color=red,points([1,2,3,4,5,6,7,8,9,10],
        [0.4,0.7,1.0,1.6,2.8,5.2,10.0,19.6,38.8,77.2]),
    key="Actual distance",color=blue,points([1,2,3,4,5,6,7,8,9,10],
        [.39,.72,1.,1.52,2.77,5.2,9.54,19.2,30.06,39.44])
    );
draw_file(terminal=svg, file_name="Titus-Bode law.svg");