integral

问题 I am trying to Initialize a Non-Integral template Constant. Please find below the code: #ifndef _EXETENDED_CLASS_H #define _EXETENDED_CLASS_H template<class T> class BaseClass { public: BaseClass(); ~BaseClass(); }; template <class T> BaseClass<T>::BaseClass() {} template <class T> BaseClass<T>::~BaseClass() {} template<class T> class ExtendedClass:public BaseClass<T> { public: typedef ExtendedClass<T>* position; static const position NULLPOSITION; ExtendedClass(); ~ExtendedClass(); private:

问题 I have a 2D array (typical size about 400x100) as shown (it looks like a trapezium because elements in the lower right are nan). For each element in the array, I want to perform a numerical integral along the column for several elements (of the order of ~10 elements). In physics language think of the colour as the magnitude of the force, and I want to find the work done by calculating th integral of Fdz. I can use a double for-loop and use trap to do the job, but are there other more

问题 I have two positive-valued vectors x,y of the same length in R. Using plot(x, y, "l",...), gives me a continuous line plot in 2 dimensions out of my finite vectors x and y. Is there a way to compute a definite integral over some range of this line plot in R? edit1: I've looked into the integrate function in R. I'm not sure however how to make a function out of two vectors to pass to it, as my vectors are both finite. edit2: For some more background, The length of x and y ~ 10,000. I've

问题 I need to calculate this integral below, using R: The q_theta(x) function I managed to do in R with quantile regression (package: quantreg ). matrix=structure(c(0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.1, 0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18, 0.19, 0.2, 0.21, 0.22, 0.23, 0.24, 0.25, 0.26, 0.27, 0.28, 0.29, 0.3, 0.31, 0.32, 0.33, 0.34, 0.35, 0.36, 0.37, 0.38, 0.39, 0.4, 0.41, 0.42, 0.43, 0.44, 0.45, 0.46, 0.47, 0.48, 0.49, 0.5, 0.51, 0.52, 0.53, 0.54, 0.55, 0.56, 0.57,

问题 I want to compute a 2-fold repeated integral (not double integral) in R, for example, where In practice, both f(x) and g(x) are quite complicated, but to run an experiment, let's simplify g(x)=1 and f(x)=cos(x), in R I use integrate to compute: > phi = function(x){integrate(function(x){cos(x)},lower=x,upper=3)[["value"]]^2} > foldintegral = integrate(phi,lower=0,upper=3) And I got this error message: Error in integrate(phi, lower = 0, upper = 3) : evaluation of function gave a result of wrong

问题 This question was migrated from Cross Validated because it can be answered on Stack Overflow. Migrated 2 years ago . With the goal of turning the following into a function, I was wondering how I can write the following double integral in terms of R codes?: ($\bar{x} = \mu$): 回答1: Assuming pi0 and pi1 implement your functions $\pi_0$ and $\pi_1$ in a vectorized way, a possible solution is: integral <- function(n, mu, s, pi0, pi1) { C <- (2 * pi)^(-n/2) C * integrate(f = function(sigmavec)

问题 Here, I would like to have your helps on implementing calculation of integral on two vector. I checked pages on integral calculation relative to R. But, I have few training on mathematics, so I still can not do that by myself. My objectives is to implement the idea of this sentence "If you plot the rate estimates by position, the genetic map is just the integral of this plot." This mean I have variables (rates, positions), each position have a rate of its own. I want to calculate the integral

问题 There is a function which determine the intensity of the Fraunhofer diffraction pattern of a circular aperture... (more information) Integral of the function in distance x= [-3.8317 , 3.8317] must be about 83.8% ( If assume that I0 is 100) and when you increase the distance to [-13.33 , 13.33] it should be about 95%. But when I use integral in python, the answer is wrong.. I don't know what's going wrong in my code :( from scipy.integrate import quad from scipy import special as sp I0=100.0

问题 I was reading Setting an int to Infinity in C++. I understand that when one needs true infinity, one is supposed to use numeric_limits<float>::infinity() ; I guess the rationale behind it is that usually integral types have no values designated for representing special states like NaN , Inf , etc. like IEEE 754 floats do (again C++ doesn't mandate neither - int & float used are left to the implementation); but still it's misleading that max > infinity for a given type. I'm trying to

问题 The full mathematical problem is here. Briefly I want to integrate a function with a double integral. The inner integral has boundaries 20 and x-2 , while the outer has boundaries 22 and 30 . I know that with Scipy I can compute the double integral with scipy.integrate.nquad . I would like to do something like this: def f(x, y): return (x ** 2 + y ** 2) res = sp.integrate.nquad(f, [[22, 30], [20, x-2]]) Is it possible? Maybe using also sympy ? 回答1: I solved with sympy : from sympy import * x,