added formula for cpp-netlib