Shape of molecular bonds

I’d like to model some simple molecules using VSEPR. With what functions should I approximate the sigma and pi bond orbitals?
A simple Gauss would do it?
I’d like to have a good precision.
Thanks in advance