***, Optimization of 2 d functions geometry={Ne} dexp=[2.0,1.0] basis={ sp,Ne,vdz;c; d,Ne,dexp(1),dexp(2) } hf mp2 eval=energy minimize,eval,dexp(1),dexp(2) ---