Using Simpson's Rule in Python

Question

I'm trying to get an array from an f(x) function this way:

array=list()

for i in range(x):
    parameter= z+(i*change)
    array=f(parameter)

Note that x is an integer, z and change are floats established in my code.

The next thing I want is to use simpson's rule using the simps function in scipy. I tried this:

Simpsons= integrate.simps(array, dx=change)

It says there is an error How can I solve this?

Christopher Paterson · Accepted Answer · 2019-10-30 18:27:13Z

1

The problem line is array=f(param). You're assigning array to the result of f, not appending it. You should do array.append(f(param)).

answered Oct 30, 2019 at 18:27

1861 silver badge9 bronze badges

Sign up to request clarification or add additional context in comments.

1 Answer 1