Unconstrained fitting of implied volatility surfaces using a mixture of normals