In order to assess the stability of DNA-DNA hybridizations-for example during PCR primer design or oligonucleotide selection for microarrays-one needs to predict the change in Gibbs free energy DeltaG during hybridization. The nearest neighbor model provides a good compromise between accuracy and computational simplicity for this task. To determine optimal combinations of reaction parameters (temperature, salt concentration, oligonucleotide length and GC-content), one would like to understand how DeltaG depends on all of these parameters simultaneously.