gaussian formula for initial con in runme.c