se = sqrt(np.diag(pcov))