f[x] = x * CDF[NormalDistribution[0,1], x] + PDF[NormalDistribution[0,1], x]