SciPy Simpson

Among the numerical models that are used to calculate the integral includes “Simpson’s” rule. Typically, we employ the logic foundation theorem to get the scale factor which requires us to use the antiderivative integration methodologies. We can obtain the integration of function f(q) using values along the direction and synthesize the Simpson’s rule with the aid of the scipy.integrate.simps() function. Essentially, it’s a Python scientific workflow module that offers built-in utilities for plenty of well arithmetic operations. The SciPy integrate inter offers several integration methods albeit with integrated for common differential equations.


The method for the implementation of the “SciPy Simpson” function will be discussed and shown in this article. We must employ an odd range of grid units while having an even combination of spans. Half of the ordinates must be multiplied by a set of parameters known as “Simpson’s multipliers” to satisfy the “Simpson’s first rule” formula, considering the examples anywhere along the specified axis and the compounded Simpson’s rule to integrate “f(q)”. The interval of “dq” is expected if “q” is none. Since Simpson’s regulation calls for an even number of intervals and there must be an integer number of samples “n”, there is an odd integer of interims “(n-1)”. The way this is managed is controlled by the “even” argument. If the samples are not evenly spaced, the function must be a polynomial having the order of “2” or less than “2” to make the result exact.


$ scipy.integrate.simps(r,q)

We mentioned the syntax of SciPy Simpson’s function in the Python language. This function has two parameters of the “integrate” function and considers the two variables for storing or passing the value in the string as “r” and “q” within the simps function of Python.

Return Value:

Employing the samples would give the integrated value of q(r) in the output screen by finding out the integration of both variables and storing the numerical values in it. The result or return value is the integrated value of Simpson’s function of those variables.

Example 1:

Now, we are familiar with the syntax and the phenomenon of working with the SciPy Simpson function. Let us start implementing it into the Python code in different scenarios. We start by having the tool first. We install the “Spyder” tool. After the installation, we start writing our code into the console file. First of all, we need the “NumPy” library in the Python source file, so we import this library first as “np”. After that, we import another library of “integrate” from the source of “SciPy”. We add some comments in between to understand what we did in each step.

After importing both “integrate” and “NumPy” libraries, we require the variables for holding the numerical value to show us the integrated phenomenon. For this purpose, we create two variables of “s” and “c” where the variable “s” is assigned with “np” as the arranged value of range “2” to “12”. Whereas the variable “c” is similarly having the NumPy “np” extension along with the arrange value of relevant range from “2” to “12”, similar with the “s” variable assigned value.

After that, we now use our main function of “integrate.simps()” on the variables “c” and “s” and assign this result to a new user-defined function of the name, “Scipy_simpson”. Here, at this step, the result is stored and the values should be integrated and stored in the “Scipy_simpson” function. In the end, to display the result, we use the “print()” function and call the function value that is stored in “Scipy_simpson”.

#importing numpy library

import numpy as np

#importing integrate of scipy

from scipy import integrate

#declaring and assigning range to variables

s = np.arange(2, 12)

c = np.arange(2, 12)

#utlizing integrate.simps() module

Scipy_simpson = integrate.simps(c, s)

#printing the Scipy_simpson function


The output of our program code for the SciPy Simpson method that we used displays the integrated value of both variables “s” and “c” as the final result of the return value as “58.5”. This value varies differently for the different stored numerical values of variables according to their ranges.

Example 2:

Let’s examine how we can utilize the same SciPy Simpson method for the usage of only one numerical value by applying the “sqrt()” function on the provided variable used. Let’s perform the implememntation of the code on our tool where we import the first two libraries of “NumPy” as “np” and “integrate” from SciPy as we used in the previous example.

Now, we declare two variables of “q” and “r” where the variable “q” is the assigned “arrange()” function with the range value of “3” and “15” and the variable “r” uses the “sqrt()” function on the value of variable “q” and store it in the variable of “r”. After assigning the values to both variables, we come the “integrate.simps()”function. We apply it to our variables “r” and “q” by defining the new function of “Scipy_simp” and store it in this function. Then, we use the “print()” function on the last step and call the “Scipy_simp” function within the “print()” function. Then, it displays the integrated relation into the final return value.

#importing numpy and integrate library os scipy

import numpy as np

from scipy import integrate

#declaring variables

q = np.arange(3, 15)

r = np.sqrt(q)

# utilizing scipy.integrate.simps() method

Scipy_simp = integrate.simps(r, q)

#printing Scipy_simp function


After completing the code when the previous code is compiled, it displays the integrated value return result on the output screen which is nearly “31.46” or both variables of “q” and “r”.


The description and the theme implementation of the SciPy Simpson method is discussed in this article. Our article illustrated two examples of the Simpsons method of SciPy to find out the relation of integrated value between two variables which are defined in the program. The first one covers the value that range from minimal to “2” while the second one ranges to minimal of “3”. In the first example, the relation values were both defined separately. But in the second example, we defined the first value and the second value is derived by the “sqrt” function for the second variable value.

About the author

Kalsoom Bibi

Hello, I am a freelance writer and usually write for Linux and other technology related content