This is a follow-up question on this one: Expected value for maximum of n normal random variable
@RobertIsrael states the following:
Presumably the $X_i$ are independent. If $\Phi$ is the standard normal cdf, $$P(\max_i X_i < \mu + t \sigma) = \prod_i P(X_i < \mu + t \sigma) = \Phi(t)^n$$ so $$ E[\max_i X_i] = \mu + \sigma \int_{-\infty}^\infty t \dfrac{d}{dt} \Phi(t)^n\ dt $$
I can see some of the ideas that lead to this formula (e.g. the multiplication because of independence) but I don't see all of the details. It would be very helpful if somebody could explain the derivation in detail... Thank you!