***,h2o !A title r=1.85,theta=104 !set geometry parameters geometry={O; !z-matrix geometry input H1,O,r; H2,O,r,H1,theta} basis=VTZ !use VTZ basis hf !closed-shell scf ccsd(t) !do ccsd(t) calculation