function vdth3(hamso,x)
    pp_euler(hamso,x,1);
    midpoint(hamso,x,1);
    rk2(hamso,x,1);