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

Need a 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

Need a 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

Need a 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

Need a hint?

Use print(result) to show the integral value.