Single integral (quad) in SciPy - Practice Problems & Coding Challenges

Challenge - 5 Problems

❓ Predict Output

intermediate

Output of a definite integral using scipy.quad

What is the output of this code that calculates the integral of f(x) = x^2 from 0 to 2 using scipy.quad?

SciPy

from scipy.integrate import quad

def f(x):
    return x**2

result, error = quad(f, 0, 2)
print(round(result, 3))

A 4.000

B 2.667

C 1.333

D 8.000

💡 Hint

Recall the integral of x^2 from 0 to 2 is (2^3)/3.

The integral of x^2 from 0 to 2 is (2^3)/3 = 8/3 ≈ 2.667. The code uses scipy.quad to compute this numerically.

❓ data_output

intermediate

Number of points used by scipy.quad for integral approximation

What is the type and approximate number of points used internally by scipy.quad when integrating f(x) = sin(x) from 0 to π?

SciPy

from scipy.integrate import quad
import numpy as np

def f(x):
    return np.sin(x)

result, error = quad(f, 0, np.pi, full_output=1)
print(result)
print(error['neval'])

A Result: 0.0, Points used: around 10

B Result: 1.0, Points used: around 5

C Result: 2.0, Points used: around 50

D Result: 3.1415, Points used: around 100

💡 Hint

The integral of sin(x) from 0 to π is 2.

The integral of sin(x) from 0 to π is exactly 2. scipy.quad uses adaptive quadrature and typically evaluates around 50 points for this integral.

🧠 Conceptual

advanced

Understanding error estimate in scipy.quad

What does the second value returned by scipy.integrate.quad represent when computing an integral?

A The number of function evaluations performed

B The time taken to compute the integral

C The relative error between two integration methods

D The estimated absolute error of the integral approximation

💡 Hint

Think about what error means in numerical integration.

scipy.quad returns two values: the integral result and an estimate of the absolute error in the result.

🔧 Debug

advanced

Identify the error in this scipy.quad usage

What error will this code raise when trying to integrate f(x) = 1/x from 0 to 1 using scipy.quad?

SciPy

from scipy.integrate import quad

def f(x):
    return 1/x

result, error = quad(f, 0, 1)
print(result)

A No error, returns a finite result

B ValueError: Integration bounds must be finite

C RuntimeWarning: divide by zero encountered in true_divide

D ZeroDivisionError

💡 Hint

Consider what happens at x=0 for 1/x.

The function 1/x has an improper integral at x=0 but scipy.quad can handle improper integrals with singularities at the boundary and returns a finite result.

🚀 Application

expert

Calculate integral of a piecewise function using scipy.quad

Using scipy.quad, what is the value of the integral of the function f(x) defined as f(x) = x for x < 1 and f(x) = 2 - x for x ≥ 1, integrated from 0 to 2?

SciPy

from scipy.integrate import quad

def f(x):
    return x if x < 1 else 2 - x

result, error = quad(f, 0, 2)
print(round(result, 3))

A 1.0

B 2.0

C 0.5

D 1.5

💡 Hint

Split the integral at x=1 and sum the two parts.

The integral from 0 to 1 of x dx is 0.5, and from 1 to 2 of (2 - x) dx is also 0.5, total 1.0.

Practice

1. What does the scipy.integrate.quad function compute in Python?

easy

A. The derivative of a function at a point

B. The definite integral of a function over an interval

C. The sum of elements in a list

D. The maximum value of a function

Solution

1. Step 1: Understand the purpose of quad

   scipy.integrate.quad is designed to calculate the definite integral of a function over a given interval.

2. Step 2: Compare with other options

   Derivative, sum, and maximum are different operations and not what quad performs.

3. Final Answer:

   The definite integral of a function over an interval -> Option B

4. Quick Check:

   Integral calculation = The definite integral of a function over an interval [OK]

Hint: Remember: quad integrates functions over intervals [OK]

Common Mistakes:

• Confusing integration with differentiation
• Thinking quad sums list elements
• Assuming quad finds max values

2. Which of the following is the correct syntax to import the quad function from scipy.integrate?

easy

A. from scipy.integrate import quad

B. import quad from scipy.integrate

C. from scipy import quad.integrate

D. import scipy.integrate.quad

Solution

1. Step 1: Recall Python import syntax

   The correct way to import a function is from module import function.

2. Step 2: Match with options

   from scipy.integrate import quad matches this syntax exactly for quad from scipy.integrate.

3. Final Answer:

   from scipy.integrate import quad -> Option A

4. Quick Check:

   Correct import syntax = from scipy.integrate import quad [OK]

Hint: Use 'from module import function' for specific imports [OK]

Common Mistakes:

• Using 'import function from module' which is invalid
• Trying to import nested attributes incorrectly
• Using dot notation in import statements wrongly

3. What is the output of this code?

from scipy.integrate import quad
result, error = quad(lambda x: x**2, 0, 2)
print(round(result, 2))

medium

A. 2.67e-1

B. 2.00

C. 2.67

D. 1.33

Solution

1. Step 1: Understand the integral being calculated

   The code integrates the function x² from 0 to 2.

2. Step 2: Calculate the definite integral

   Integral of x² dx from 0 to 2 is (x³/3) evaluated from 0 to 2 = (8/3) - 0 = 2.6667 approximately.

3. Step 3: Check the printed output

   The code rounds the result to 2 decimal places: 2.67.

4. Final Answer:

   2.67 -> Option C

5. Quick Check:

   Integral of x² from 0 to 2 = 2.67 [OK]

Hint: Integrate x² from 0 to 2 = 8/3 ≈ 2.67 [OK]

Common Mistakes:

• Confusing integral result with function value at 2
• Rounding incorrectly
• Mixing error estimate with result

4. Identify the error in this code snippet:

from scipy.integrate import quad
result = quad(x**3, 1, 3)
print(result)

medium

A. quad does not return any value

B. The integration limits should be reversed

C. The print statement is missing parentheses

D. The function to integrate must be passed as a callable, not an expression

Solution

1. Step 1: Check the function argument to quad

   The argument x**3 is an expression, not a function. quad expects a callable function like lambda x: x**3.

2. Step 2: Confirm other parts are correct

   Integration limits and print statement are correct syntax. quad returns a tuple.

3. Final Answer:

   The function to integrate must be passed as a callable, not an expression -> Option D

4. Quick Check:

   Function argument must be callable [OK]

Hint: Pass a function or lambda, not a direct expression [OK]

Common Mistakes:

• Passing expressions instead of functions
• Swapping integration limits unnecessarily
• Forgetting quad returns a tuple

5. You want to calculate the integral of sin(x) from 0 to π using scipy.integrate.quad. Which code snippet correctly computes this and prints the integral value rounded to 3 decimals?

hard

A. from scipy.integrate import quad
import math
result, error = quad(math.sin, 0, math.pi)
print(round(result, 3))

B. from scipy.integrate import quad
result, error = quad('sin(x)', 0, 3.14)
print(result)

C. import scipy
result = scipy.integrate.quad(lambda x: sin(x), 0, 3.14)
print(round(result, 3))

D. from scipy.integrate import quad
result = quad(lambda x: math.sin(x), 0, 3.14)
print(round(result, 3))

Solution

1. Step 1: Check function import and usage

   from scipy.integrate import quad import math result, error = quad(math.sin, 0, math.pi) print(round(result, 3)) imports math and passes math.sin as a callable, which is correct.

2. Step 2: Verify integration limits and print

   Limits 0 to π are correct using math.pi. Result is rounded to 3 decimals before printing.

3. Step 3: Identify errors in other options

   from scipy.integrate import quad result, error = quad('sin(x)', 0, 3.14) print(result) passes a string instead of a function. import scipy result = scipy.integrate.quad(lambda x: sin(x), 0, 3.14) print(round(result, 3)) misses math.sin import and uses undefined sin. from scipy.integrate import quad result = quad(lambda x: math.sin(x), 0, 3.14) print(round(result, 3)) does not unpack the tuple returned by quad.

4. Final Answer:

   from scipy.integrate import quad import math result, error = quad(math.sin, 0, math.pi) print(round(result, 3)) -> Option A

5. Quick Check:

   Use callable math.sin and unpack quad result [OK]

Hint: Import math and pass math.sin as function, unpack result [OK]

Common Mistakes:

• Passing string instead of function
• Not unpacking quad's tuple output
• Using undefined sin without import