Add formula for simpleitk