Hola. Yo veo bien tu resolución. Solo ojo que si pones vector de un lado, lo del otro lado también debe ser vector (linea 3 al final).
Hay dos formas de implementar la ley de gauss, podes hacerla integrando o podes plantear la int [(D2-D1).nda] =QL = int(sigma da) - > Int[(D2-D1).n -sigma]da =0, y como el Área de la pastilla la tomo como quiero lo de adentro es nulo, entonces la condición se ve como (D2-D1).n = sigma, por eso mantengo los versores.
eso te resuelve la duda?