Fitting data to gaussian distributions

Hi!

I need to fit a data that shows up as two gaussians partially
superimposed to the corresponding gaussian distributions, i.e.

data=c(rnorm(100,5,2),rnorm(100,-6,1))

I figured it out how to do it with mle or fitdistr when only one
gaussian is necessary, but not with two or more. Is there a function in
R to do this?

Thank you very much in advance,

Luis