integrate √(1 -(y/b)²)*dy