Questions tagged [quad]

Quad is a scipy package which is used for 1D integration in Python

135 questions
0
votes
1 answer

Can't integrate with dblquad

So I want to integrate a double integral with constants in it, like a, b, etc where the user can asign the value of this constants: The limits of the integral are x[0,1] and y[-1,2] import numpy as np import scipy.integrate as integrate def…
0
votes
0 answers

triple integration over arrays python

I have the following triple integral: the integral the data are: t 1D array(size 124) , pr 1D array(size 10), lat 1D array(241) and lon 1D array(size 480) V 4D array (124, 10, 241, 480) I want to integrate v over t, pr and lon at each lat So, I…
Mag
  • 11
  • 2
0
votes
0 answers

Bessel function using scipy.quad integration not working

I'm trying to plot a Bessel function in 2D using the scipy.quad integration method. My problem: the plot only has a tiny dot in the middle because my lowest value is far too low. Here is my code: import numpy as np import scipy.integrate as…
0
votes
1 answer

error in plotting a function, including integrate.quad

I want to plot a function which has integrate.quad in it. The function has true float value in the interval. import matplotlib.pyplot as plt import numpy as np import scipy.integrate as integrate def S(M): s = integrate.quad(lambda k:…
Smtl
  • 49
  • 6
0
votes
1 answer

How to propagate error when using scipy quad on a spline of data with measurement error?

I have a data set with N points which I fit a spline to and integrate using scipy.integrate.quad. I would like to use the N associated measurement errors to put an error estimate on the final integral value. I originally tried to use the…
handroski
  • 81
  • 2
  • 3
  • 15
0
votes
1 answer

Conversion of symbolic expression to numeric one for use in quad - use lambdify?

I want to convert an expression containing symbolic variables to a numeric one so that the expression may be subsequently used in an integration method 'quad'. import numpy import math as m import scipy import sympy #define constants …
CAF
  • 329
  • 3
  • 14
0
votes
0 answers

Numerical integration method 'quad' derivative error

I have defined a series of functions below with the end function fA being inserted for a numerical integration. The integration is with respect with one variable so the other arguments are passed as numeric so that the integration method (quad) may…
CAF
  • 329
  • 3
  • 14
0
votes
1 answer

"TypeError" message while performing quad integration method of scipy.integrate with mpmath functions

I am trying to calculate two integrals using scipy.integrate.quad. However, since gamma functions with negative first parameter are not defined in scipy, I had to choose the version from mpmath. After running the following code, from scipy.integrate…
Rebel
  • 472
  • 8
  • 25
0
votes
1 answer

IndexError: tuple index out of range, error when trying to integrate the fucntion

I am trying to calculate a value of function f_integ, which is a result of integration of function f from 0 to x_v. f = lambda x : x + 1 def f_integ(x_array): int_result = np.zeros_like(x_array) indexes =…
ani87
  • 17
  • 5
0
votes
2 answers

Minimizing a function using python for data fitting

I have a function as the following q = 1 / sqrt( ((1+z)**2 * (1+0.01*o_m*z) - z*(2+z)*(1-o_m)) ) h = 5 * log10( (1+z)*q ) + 43.1601 I have experimental answers of above equation and once I must to put some data into above function and solve…
Ethan
  • 180
  • 2
  • 15
0
votes
1 answer

How do you invert scipy quad function to pass in an arc length and receive a value for (t)?

If I were to run something like the code below where I pass in the minor and major arc of an ellipse as 3.05 and 2.23 with the 50 degree angle formed by the arc, how would I be able to take the output of 2.531432761012828 as the arc length and pass…
Fairly Factual
  • 161
  • 3
  • 13
0
votes
0 answers

Nonlinear magnetic model fit -with a definite integral- using scipy.integrate.quad and lmfit

I would like to fit an hysteresis curve, with a superparamagnetic behavior, using a magnetic model which includes a Langevin function and a pair distribution function [1]. To fit with this equation I must solve a definite integral. I was trying to…
Ernesto
  • 21
  • 4
0
votes
2 answers

python scientific programming why blank occur when plotting with right results?

I want to plot the integration results as described below , but it seems turned out to be a blank sheet, what' s the reason ? please help me! # -*- coding: utf-8 -*- import matplotlib.pylab as plt import numpy as np import scipy as sp from…
0
votes
3 answers

Python: Evaluating Integral for array

import numpy import matplotlib.pyplot as plt from scipy import integrate def f(x,y): return x*y + x**2 def integral(x,y): I = integrate.quad(f, 0, x, args=(y,))[0] return I def gau(x,y): return (1+x)*integral(x,y) xlist =…
user7248647
  • 67
  • 1
  • 2
  • 9
0
votes
0 answers

Scipy quad integrate with cdf as differential

I'm quite stuck at the moment with a little function I write. I want to embed the following function into my program: where is the convolution of of the cumulative distribution function of a poisson distribution. If I got it right, I'd calculate…
FriendlyGuy
  • 331
  • 5
  • 14
1 2 3
8
9