Si no consideras que el cilindro es finito y desprecias los efectos de bordes, puedes plantear que la superficie gaussiana más adecuada es un cilindro coaxial al cilindro que contiene las cargas, y su área queda expresada: A = 2π*L*r.
Luego, la ecuación queda:
E*2π*L*r = (qp + qn)/ε0.
Haces pasaje de divisor como factor, y queda:
E*2πε0*L*r = qp + qn,
haces pasaje de término, y queda:
E*2πε0*L*r - qp = qn.
Luego, sustituyes la expresión: ε0 = 1/(4π*k), y queda:
qn = E*L*r/(2*k) - qp.
Luego, reemplazas datos (observa que empleamos el Sistema Internacional de Unidades de Medida,
y observa también que el campo es "atractivo" y tiene dirección y sentido hacia el eje del cilindro):
qn = -9876,9*0,12*0,12/(2*9*109) - 25*10-9 = -7,90152*10-9 - 25*10-9 = - 32,90152*10-9 c = -32,90152 nc.
Espero haberte ayudado.