Iterated residues and multiple Bernoulli polynomials On the representation of the Kummer surface of fourth order with six nodal points with the Göpel biquadratic relation between four theta functiions with two variables. 