# -*- coding: utf-8 -*-
"""
Created on Fri Jun  4 15:05:18 2021

@author: Daniel Powell
"""

import numpy as np
import scipy.integrate


def f(x):
    return 1.0 / (np.sqrt(1 - (x ** 2.5)))

(ans, abserror, infodict) = scipy.integrate.quad(f, 0, 1, full_output=1,
                                                 epsabs=6e-14, epsrel=6e-14)

print(ans)