I have independent variables $ X_1, X_2,\ldots,X_n $ with normal distribution on range $ [0,1] $ . Next, variables $ Z_i $ are created according to this formula $ Z_i = - \frac{1}{\lambda} \ln(1-X_i) $ , where $ \lambda > 0 $. How to find probability density function of $ W = \max (Z_1 , \ldots , Z_n ) $ ?
After comments : So, I have a formula $ f_Z(z) = | \frac{d}{dz} (1-e^{-\lambda z}) | f_X(1-e^{-\lambda z}) $ , but $ f_X (x) = 1 $ , so $ f_Z(z) = 1 + \lambda e^{-\lambda z} $ And it represents all $ Z_i $
For every maximum $W$ of random variables, one has, for every $w$, $$ [W\lt w]=\bigcap_{i=1}^n[Z_i\lt w]. $$ If $(Z_i)_{1\leqslant i\leqslant n}$ is i.i.d. and distributed as $Z$, then, for every $w$, $$ P(W\lt w)=\prod_{i=1}^nP(Z_i\lt w)=P(Z\lt w)^n. $$ In your case, considering $X$ distributed as $X_1$, $$ [Z\lt w]=[\ln(1-X)\gt-\lambda w]=[X\lt1-\mathrm e^{-\lambda w}], $$ hence, for every $w\gt0$, $$ F_W(w)=F_Z(w)^n=F_X(1-\mathrm e^{-\lambda w})^n. $$ If every $X_i$ is uniform on $[0,1]$ then $F_X(x)=x$ for every $x$ in $(0,1)$ hence, for every $w\gt0$, $$ F_W(w)=(1-\mathrm e^{-\lambda w})^n. $$ Then the PDF $f_W$ is obtained by differentiation, that is, $$ f_W(w)=n\lambda\mathrm e^{-\lambda w}(1-\mathrm e^{-\lambda w})^{n-1}\mathbf 1_{w\gt0}. $$