$method=[hf,mp2,ccsd,qci,bccd,multi,mrci,acpf,rs3]              !some methods
basis=vdz                                                      !basis
geometry={o;h1,o,r;h2,o,r,h1,theta}                            !geometry
gdirect                                                        !direct option
r=1 ang,theta=104                                              !bond length and angle
do i=1,#method                                                 !loop over methods
$method(i)                                                     !run method(i)
e(i)=energy                                                    !save results in variables
dip(i)=dmz
enddo
table,method,e,dip                                             !print table of results