Single integral (quad) in SciPy - Mini Project: Build & Apply

Calculate a Single Integral Using scipy.quad

📖 Scenario: Imagine you want to find the area under a curve for a simple function, like f(x) = x^2, between two points. This is a common task in science and engineering to understand total quantities.

🎯 Goal: You will write a Python program that uses the scipy library's quad function to calculate the definite integral of f(x) = x^2 from 0 to 3.

📋 What You'll Learn

Create a function f that returns the square of its input

Set the integration limits as 0 and 3

Use scipy.integrate.quad to calculate the integral

Print the result of the integral

💡 Why This Matters

🌍 Real World

Calculating areas under curves helps in physics for work done, in biology for growth models, and in finance for risk assessment.

💼 Career

Knowing how to compute integrals programmatically is useful for data scientists and engineers working with continuous data and models.

1

Define the function to integrate

Create a function called f that takes one input x and returns x**2.

SciPy

# Your code here
def f(x):
    # Return x squared
    pass

Hint

Use def to define the function and return x**2 to return the square.

2

Set the integration limits

Create two variables called lower_limit and upper_limit and set them to 0 and 3 respectively.

SciPy

def f(x):
    return x**2

# Your code here

Hint

Just assign the numbers 0 and 3 to the variables lower_limit and upper_limit.

3

Calculate the integral using scipy.quad

Import quad from scipy.integrate. Then use quad with the function f and limits lower_limit and upper_limit. Store the first returned value in a variable called result.

SciPy

def f(x):
    return x**2

lower_limit = 0
upper_limit = 3

# Your code here

Hint

Use from scipy.integrate import quad and call quad(f, lower_limit, upper_limit). Assign the first output to result.

4

Print the integral result

Write a print statement to display the value stored in result.

SciPy

from scipy.integrate import quad

def f(x):
    return x**2

lower_limit = 0
upper_limit = 3

result, error = quad(f, lower_limit, upper_limit)

# Your code here

Hint

Use print(result) to show the integral value.

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