Single integral (quad) in SciPy

Introduction

We use single integrals to find the total amount or area under a curve. It helps us understand how things add up continuously.

Calculating the total distance traveled when speed changes over time.

Finding the area under a curve in a graph, like total sales over a period.

Estimating the accumulated quantity, such as total rainfall over a day.

Computing probabilities in continuous distributions.

Determining work done when force varies along a path.

Syntax

SciPy

from scipy.integrate import quad

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

function is the function you want to integrate.

lower_limit and upper_limit define the range of integration.

Examples

Integrate x squared from 0 to 1.

SciPy

from scipy.integrate import quad

def f(x):
    return x**2

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

Integrate x cubed from 1 to 2 using a lambda function.

SciPy

result, error = quad(lambda x: x**3, 1, 2)

Integrate sine function from 0 to π.

SciPy

import math
result, error = quad(math.sin, 0, math.pi)

Sample Program

This program calculates the integral of e^-x² from 0 to 1, which is a common function in probability and statistics.

SciPy

from scipy.integrate import quad
import math

def integrand(x):
    return math.exp(-x**2)

result, error = quad(integrand, 0, 1)
print(f"Integral result: {result}")
print(f"Estimated error: {error}")

Important Notes

The quad function returns two values: the integral result and an estimate of the error.

Always check the error to understand how accurate the result is.

You can integrate any function that takes a single variable as input.

Summary

Single integrals help find total amounts or areas under curves.

Use scipy.integrate.quad to calculate integrals easily in Python.

The function returns both the integral value and an error estimate.

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