Add numpy π calculation