$sigma = sqrt( ((@set_a-1)*$var_a + (@set_b-1)*$var_b) / $df );<br />

$p = ($avg_a - $avg_b) / ($sigma * sqrt( 1/@set_a + 1/@set_b));</p>