The common Numerical Integration techniques.
πŸ‘︎ 3k
πŸ“°︎ r/mathmemes
πŸ’¬︎
πŸ‘€︎ u/12_Semitones
πŸ“…︎ Apr 30 2020
🚨︎ report
[F#] Numerical Integration in 30 lines, plus 100 lines of tests /r/CritiqueMyCode/comment…
πŸ‘︎ 2
πŸ“°︎ r/codereview
πŸ’¬︎
πŸ‘€︎ u/CorrettoSambuca
πŸ“…︎ Sep 25 2020
🚨︎ report
[F#] Numerical Integration in 30 lines, plus 100 lines of tests

I am a mathematician by profession, and am trying my hand at F#.

This is not my first exercise, but I would like some feedback especially on the testing part - I'll soon be building a much bigger system for more sophisticated numerics.

Pastebin: https://pastebin.com/R6z08kER

Gist: https://gist.github.com/OrlandoRiccardo/543fb187e6f00589a7f99cbd8884623b

πŸ‘︎ 3
πŸ“°︎ r/CritiqueMyCode
πŸ’¬︎
πŸ‘€︎ u/CorrettoSambuca
πŸ“…︎ Sep 25 2020
🚨︎ report
|Numerical Integration| What is Simpson's Rule? [Intuition] :) v.redd.it/m572spvr18r41
πŸ‘︎ 548
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/Aravindh_Vasu
πŸ“…︎ Apr 06 2020
🚨︎ report
[College math: Numerical integration] How can I remove singularities from integral?

I've got an integral from 2 to infinity of 1/[x*ln^2(x+sqrt(x))] to integrate numerically using matlab. First I removed singularities and made it into a proper integral by substitution y=1/x, which gave me integral from 0 to 1/2 of 1/[y*ln^2(1/y+1/sqrt(y))], which still have a singularity and I don't know how to remove it in order to integrate this numerically. Can you help me with this?

πŸ‘︎ 3
πŸ“°︎ r/HomeworkHelp
πŸ’¬︎
πŸ‘€︎ u/rikoslav
πŸ“…︎ Jul 09 2020
🚨︎ report
AP Calc: Numerical Integration with Tables

Self studying BC here. For some rather inexplicable reason, I just suck at numerical integration (Riemann sums, midpoint, trapezoidal, etc). This is specifically the case when the info is in tables with unequal intervals. It always costs me two points when it comes up. I've looked for resources but there's nothing specifically geared to applying those techniques with tabular data. The FRQ solution guides don't make much sense to me.

Could someone give me some resources to work with or guide me through it with an example? That would be much appreciated.

πŸ‘︎ 4
πŸ“°︎ r/APStudents
πŸ’¬︎
πŸ‘€︎ u/BlueMonkeys090
πŸ“…︎ Apr 30 2020
🚨︎ report
A Video Lecture on Numerical Double Integration
πŸ‘︎ 14
πŸ“°︎ r/mathematics
πŸ’¬︎
πŸ‘€︎ u/dishantpandya777
πŸ“…︎ Apr 08 2020
🚨︎ report
|Numerical Integration| What is Simpson's Rule? [Intuition] youtube.com/watch?v=DdNAc…
πŸ‘︎ 3
πŸ“°︎ r/VisualMath
πŸ’¬︎
πŸ‘€︎ u/FunVisualMath
πŸ“…︎ May 08 2020
🚨︎ report
Ideas for parameterizing this curve in the complex plane and calculating its length by (numerical) contour integration?
πŸ‘︎ 2
πŸ“°︎ r/mathematics
πŸ’¬︎
πŸ‘€︎ u/cl214
πŸ“…︎ May 11 2020
🚨︎ report
Numerical integration is better in certain situations.
πŸ‘︎ 57
πŸ“°︎ r/mathmemes
πŸ’¬︎
πŸ‘€︎ u/12_Semitones
πŸ“…︎ Mar 18 2020
🚨︎ report
Numerical Integration of Differential Equations

When can a differential equation be solved through analytical integration, and when can it be solved ONLY through numerical integration. Or in other words: What boxes must a differential equation need to "tick" before I know it can be analytically integrated?

πŸ‘︎ 27
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/GuyNew1234
πŸ“…︎ Dec 06 2019
🚨︎ report
Error bounds of numerical integration

I am trying to understand how to use the error bound of a trapezoidal numerical integration.

If the error bound is

abs(E) <= (K(b-a)^3 ) / (12n^2)

and

f''(x) <= K

How do I find (or reasonably estimate) K if I don't have the exact f(x) ?

Thanks.

Edit: To elaborate, I'm trying to propagate errors through a rather complex equation that includes a trapezoidal integration. I have the error due to my 6 measured variables accounted for, but I haven't been able to figure out how to account for the error to to the numerical integration. Any help on this would be greatly appreciated.

πŸ‘︎ 2
πŸ“°︎ r/learnmath
πŸ’¬︎
πŸ‘€︎ u/Proton_Driver
πŸ“…︎ Jan 22 2020
🚨︎ report
Numerical Methods in C++ Part 1: Newton-Cotes Integration thoughts-on-cpp.com/2019/…
πŸ‘︎ 38
πŸ“°︎ r/cpp
πŸ’¬︎
πŸ‘€︎ u/bmahr
πŸ“…︎ Apr 17 2019
🚨︎ report
Very weird results with numerical integration using trapezoid method - Online Technical Discussion Groupsβ€”Wolfram Community

The graph I got using the code in the link is VERY peculiar to say the least. I have no idea what's causing it.

https://community.wolfram.com/groups/-/m/t/1840568

πŸ‘︎ 2
πŸ“°︎ r/learnmath
πŸ’¬︎
πŸ‘€︎ u/MatheMaticaHelp_
πŸ“…︎ Dec 12 2019
🚨︎ report
Numerical integration and optimization

I want to solve something like min_x \int f(x,y)dy automatically. Not sure if there's a way to do this in JuMP with a numerical integration method or if I have to write out the integral approximation.

πŸ‘︎ 10
πŸ“°︎ r/Julia
πŸ’¬︎
πŸ“…︎ Jun 18 2019
🚨︎ report
Numerical Integrators Galore: Harmony of the Spheres, an n-body gravity simulator featuring 12 different numerical integration routines ranging from euler to the 12'th order Runge Kutta Nystrom routine thehappykoala.github.io/H…
πŸ‘︎ 44
πŸ“°︎ r/Physics
πŸ’¬︎
πŸ‘€︎ u/TheVeryHappyKoala
πŸ“…︎ Mar 28 2019
🚨︎ report
Numerical Methods in C++ Part 2: Gauss-Legendre Integration thoughts-on-cpp.com/2019/…
πŸ‘︎ 26
πŸ“°︎ r/cpp
πŸ’¬︎
πŸ‘€︎ u/bmahr
πŸ“…︎ Apr 25 2019
🚨︎ report
Been learning numerical integration recently
πŸ‘︎ 396
πŸ“°︎ r/MathJokes
πŸ’¬︎
πŸ‘€︎ u/esquared722
πŸ“…︎ Nov 01 2018
🚨︎ report
[1902.00648] Numerical Integration Method for Training Neural Network arxiv.org/abs/1902.00648
πŸ‘︎ 4
πŸ“°︎ r/MachineLearning
πŸ’¬︎
πŸ‘€︎ u/ihaphleas
πŸ“…︎ Feb 05 2019
🚨︎ report
Is there a way to express the position of a celestial body in orbit in r(t) and ΞΈ(t) without numerical integration?

As far as I know, there's no mathematical way to analytically express the position of a celestial body orbiting another body in a two-body problem at time t, just like calculating precisely the circumference of an ellipse.

Am I wrong? If I'm right is there an equation that approximates r(t) and ΞΈ(t)?

If it helps, I'm trying to plot the density wave theory for my students so they can see where the individual stars go.

Edit: I'm talking about elliptical orbits because there are no perfectly circular and parabolic orbits in real life.

πŸ‘︎ 2
πŸ“°︎ r/astrophysics
πŸ’¬︎
πŸ‘€︎ u/CelestialSegfault
πŸ“…︎ Aug 07 2019
🚨︎ report
New rigorous numerical integration in Arb fredrikj.net/blog/2017/11…
πŸ‘︎ 79
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/fredrikj
πŸ“…︎ Nov 25 2017
🚨︎ report
So We're Calculating Pi by Numerical Integration? [OC] imgur.com/a/t4gFbsF
πŸ‘︎ 16
πŸ“°︎ r/dataisbeautiful
πŸ’¬︎
πŸ‘€︎ u/equationsofmotion
πŸ“…︎ May 19 2018
🚨︎ report
Could you suggest a numerical integration package that can work on a function sampled on a grid?

I'm trying to look for something similar to scipy.integrate. There are several packages out there performing adaptive and other advanced integrations provided a generic function, but surprisingly I can't find one that integrates a function sampled on a grid.

I understand that I can write a Simpson's rule function myself, but if there's a package that has verified tests etc that'll be great. For some reason I'm not being able to import scipy correctly on the computer cluster, so that's not an option.

πŸ‘︎ 3
πŸ“°︎ r/Julia
πŸ’¬︎
πŸ‘€︎ u/slipnips
πŸ“…︎ Dec 03 2018
🚨︎ report
Numerical Integration with pymc3?

I am building a model of random variables in pymc3 that involves a numerical integration of some of my variables (and some data arrays), for which there is not an analytical solution. In the original pymc, I can use numpy.trapz or scipy.quad, but these do not work with the theano variables of pymc3. Surely there is a way to do non analytic integrals? Does anyone know how to accomplish this? For an example, see my unanswered question on StackOverflow.

πŸ‘︎ 3
πŸ“°︎ r/learnpython
πŸ’¬︎
πŸ‘€︎ u/stellar_guy
πŸ“…︎ Apr 06 2019
🚨︎ report
Need help in doing double integration of a function f(x,y) using two different numerical methods

For my numerical methods project I was given a task to do double integration of a function f(x,y), where I need to use the Simpson's 3/8 rule in respect to x, and then use the Composite Trapezoidal rule in respect to y.

I've implemented the functions for both two methods but I've been struggling how to use those methods to integrate x while treating y as a constant, and to integrate y while treating x as a constant. For example, I would like my Trapezoidal function to take in the function f(x,y) and do the method for the y variable, treat x as a constant, then have an output value having the form of a function f(x,y).

I tried searching for a long time and couldnt really come up with any example anywhere so far, and I would appreciate the help.

this is the Trapezoidal Rule function. I tried putting f(x(i),y) in place of f(x(i)) for example, etc. but so far Im stuck.

EDIT: figured it out! Newton's 3/8 Method

Composite Trapezoidal Method

And this is the function to run them both.

πŸ‘︎ 4
πŸ“°︎ r/matlab
πŸ’¬︎
πŸ‘€︎ u/krsuma
πŸ“…︎ Jan 29 2019
🚨︎ report
Numerical Integration & Stability in Structural Dynamics

Hi All,

I'm working on a simple structural dynamics FEA program in MATLAB. For now I am using a thin rectangular plate as the object to be analyzed, and I'm meshing it with rectangular plate elements with three degrees of freedom per node (one translation, two rotation). I have successfully implemented a normal modes solution that works seamlessly, however I am now working on a transient base excitation solver and hitting a hard wall with the numerical integration. My solutions are blowing up except for the case of very small timesteps (~1e-7 sec), even when I am using only ~100 elements (~350 DOFs).

Questions:

  1. How do I calculate the stable time increment for an MDOF system using a time-marching scheme, e.g. when using the Newmark-beta method to simulate a dynamic transient response?
  2. If there is a better integration method for transient response in structural dynamics, in terms of greater stability and quicker computation (likely from allowing larger stable timesteps), what would you suggest?

Some background on my project: The test/control plate for my program is 6" x 4" x 0.125". I'm integrating the response in modal coordinates and only using contributions from the 10 lowest modes, with the highest having a frequency of 23,142 rad/s (3,863 Hz). The excitation is a 1 lb, 5 ms half-sine shock pulse on an interior node ((x, y) = (1.75 in, 1.75 in)) with the four corner nodes pinned. I am currently using the Newmark-beta method for integrating the response with beta = 1/6 and gamma = 1/2. These are all mostly arbitrary, but this is my starting point.

I've run a bit of a numerical stability case study using this plate, and I can easily see that the stable time increment is a function of (modal) damping. It follows the behavior of the function f(x) = x*e^(-x) or something very similar. At any rate, I cannot seem to find an answer on Google on stable time increments and don't have a formal academic background in numerical methods.

Any and all help is greatly appreciated!

πŸ‘︎ 2
πŸ“°︎ r/numerical
πŸ’¬︎
πŸ‘€︎ u/GeeFLEXX
πŸ“…︎ Mar 04 2019
🚨︎ report
2D Numerical Integration on a sphere

I have to do some surface integrals over a sphere as part of my thesis. What methods do people generally use for that? I've been using trapezium rule up to now (quite easy to extend from a line to a rectangle) but I don't know if that works on a sphere.

πŸ‘︎ 3
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/jm_13
πŸ“…︎ Aug 18 2018
🚨︎ report
Numerical Order of Integration

Hello all! I had a question about numerically calculating a 2D integral over a rectangular region. The TLDR is, if I use nested 1D methods, does the order of these methods matter, either in terms of the value of the result or the speed at which the integral is calculated?

The integral is of the form \int f(x,y) dx dy, where the bounds of y are from 0 to a, and the bounds of x are from 0 to infinity. The specific form of f is rather complicated; it takes a page of text for me to define all the terms in it. But key properties is that f decays exponentially in x and y (as in both e^-x^2 and e^-y^2 ), is very well behaved in x (holding y constant), and has an oscillatory region in y (holding x constant). I can quickly calculate the location of these regions ahead of time, and their location is independent of x. Speed is an important factor here. f is also even in both x and y.

These behaviors lead me to think I should use nested 1D integration schemes: Gaussian-Hermite to integrate over x, and a recursive adaptive scheme over y (first trying Simpson’s, then going from there), while calculating the oscillatory region separately from the rest of the domain. The Gaussian-Hermite quadrature would also be adaptive; the basic plan is compute it using more and more abscissa until the relative error is good enough (same idea for adaptive Simpson’s obviously).

Given my strategy, there is a choice to make: I can numerically integrate over x given a value of y, and then integrate that over y, or do the reverse. Or, put another way, I can calculate

G(x) = \int f(x,y) dy using adaptive Simpson’s, and then \int G(x) dx using Gauss-Hermite

Or

H(y) = \int f(x,y) dx using Gauss-Hermite, and then \int H(y) dy using adaptive Simpson’s.

Of course, changing the order of integration wouldn’t change the result if I were to analytically calculate the integral. But to carry out the integration numerically, is there an obvious right choice, either for calculating a more accurate value or for converging on that value quickly ? My gut feeling is that one may be more efficient than the other, but it’s not obvious to me which one it would be.

πŸ‘︎ 8
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/cdstephens
πŸ“…︎ Jul 24 2018
🚨︎ report
Help with numerical integration in matlab.

Hi I want to use matlab to perform numerical integration on a function that looks something like

https://preview.redd.it/04cwipndery11.png?width=769&format=png&auto=webp&s=18ceebd1b2d6b374438c1f15a58df72c1af4bd54

where a,b and c are constants between the limits of (c,a) but I want to iterate over various values of c.

The following is the for loop that I currently have written to do this

for Ed=0.1:0.1:10
    fun = @(T) (1 - b^2 .*T/a)*(1./T.^2)*(exp(-0.5.*((T-6.1875)/2.11612).^2))*(0.8.*T/Ed);
    value = integral(fun,Ed,T_hat);
    dpa_rate_linear  = [dpa_rate_linear;value];
end

but I keep getting the following error

Error using *

Inner matrix dimensions must agree.

Error in @(T) (1 - b^2 .*T/a)*(1./T.^2)*(exp(-0.5.*((T-6.1875)/2.11612).^2))*(0.8.*T/Ed)

Error in integralCalc/iterateScalarValued (line 314)

fx = FUN(t);

Error in integralCalc/vadapt (line 132)

[q,errbnd] = iterateScalarValued(u,tinterval,pathlen);

Error in integralCalc (line 75)

[q,errbnd] = vadapt(@AtoBInvTransform,interval);

Error in integral (line 88)

Q = integralCalc(fun,a,b,opstruct);

Could someone help me with my mistake?

πŸ‘︎ 2
πŸ“°︎ r/matlab
πŸ’¬︎
πŸ‘€︎ u/wonder_bro
πŸ“…︎ Nov 16 2018
🚨︎ report
Numerical surface integration for the magnetic field outside of a bar magnet

I'm attempting to write a FORTRAN program that calculates the magnetic field, B, at any point outside of a bar magnet. I have asked about this already but I figured my questions weren't clear enough.

I'm going to use a first order euler scheme, where each side of the bar magnet is split into small cells, each with centres at (xi,yi,zi). I know I can ignore all of the sides that have any z values, and just focus on the top and bottom sides that are orientated in the x-y plane. So the method says this:

[;\int f(x,y,z)dS;] = [;\Delta S;] [;\cdot;] [;\sum;] f( [;x_{i};] [;y_{i};] [;z_{i};] )

where the integral is over the surface S, and the summation is over i.

Delta S, each area, is given by [;\hat{n};] [;\cdot;] [;d\vec{S};] , so if the cell is oriented in the x-y plane it's just [;\Delta x;] [;\cdot;] [;\Delta y;] .

Here is a screenshot of the specific method instructions with a figure that demonstrates it: http://prntscr.com/idnaik

The function for the magnetic field is this:

[;\vec{B};] [;(\vec{r});] = [;\frac{\mu _{0}}{4\pi };] [;\cdot;] [;\int;] [;\frac{(r-{r}'M(r)[;\cdot[;\hat{n}];}{|(r-{r}')|^{3}};]

Where the integral is over the surface S

The ultimate aim is to use python to then trace a field line in 3d, and hopefully get the shape you'd expect for a bar magnet. I'm struggling to understand this method. I've tried to construct a flow chart, but can't get very far. Any help to understand it would be appreciated, and also any help with the flow chart would be great.

Here's my flow chart: http://prntscr.com/idnbn9

Thanks (I posted this a minute ago, but the LaTeX wasn't showing up- i'm sure it's formatted correctly)

πŸ‘︎ 3
πŸ“°︎ r/AskPhysics
πŸ’¬︎
πŸ‘€︎ u/EnvironmentalDish
πŸ“…︎ Feb 12 2018
🚨︎ report
[Undergraduate, Numerical Integration] Richardson method for error control (for integral approximation using composite Simpson's 1/3 rule)

Hello, everyone. :)

There's this handout I was given for my Numerical Methods class, and it's asking me to write a program (which is not what I need help with) to calculate the integral of e^(-x^2) dx from 0 to 1 (using the composite Simpson's 1/3 rule) with a given tolerance, where this error control is achieved using the "Richardson method" (which is what I need help with), but when I research that online, I either just find what to seems the be the "differentiation version", with the central divided difference, [f(a+h) - f(a-h)] / (2h), or a bunch of stuff that doesn't seem like what I am looking for (but maybe it is, and I'm just not noticing).

Could someone please point me to a resource or resources that are as close as possible to what I need or provide me with an explanation themselves?

Any input would be GREATLY appreciated!

P.S.

Do let me know if you need more information from me!

πŸ‘︎ 2
πŸ“°︎ r/learnmath
πŸ’¬︎
πŸ‘€︎ u/alfredkaminski
πŸ“…︎ Dec 03 2018
🚨︎ report
When can Simpson's 3/8 rule be applied? (Numerical Integration)

I sat a test today and was asked to numerically integrate a data set with had 5 or 6 data points (can't remember which) using Simpsons 3/8 rule. How is this possible as the rule has place for only 4 data points?

πŸ‘︎ 4
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/glover261
πŸ“…︎ Jun 09 2018
🚨︎ report
[HELP] Numerical integration over a discrete dataset - trapz doesn't cut it

I am calculating a gas diffusional profile in a solid sphere for which the first and second derivatives of the curve are both negative. This causes a systematic bias such that linear interpolation/trapz causes volume calculation errors to rise asymptotically under certain circumstances. I have been told that Romberg integration is a good antidote to the problem, however, it appears that a requirement of the Romberg method is to have a smooth, continuous function. Unfortunately, this curve can not be described by such a function.

I am wondering if someone more skilled than myself could help me come up with an interpolation scheme (piecewise?) that would give me the continuous function that could be fed into a Romberg script.

I have also tried using Simpson's rule integration, but it seems my implementation is causing significant errors as well. I am probably doing something stupid wrong, but I have been staring at this stuff too long to really even tell. The approach here is to calculate and add the volume contributed by the concentration in each sphere between two spatial nodes, where v is the height of the concentration curve at radial position r, j is the index for spatial nodes and n is the number of spatial nodes:

vol = simpsons([v(1) v(1)],0,r(1),[]) * (4/3 * pi * r(1)^3);
for j = 2:n
    vol = vol + (4/3 * pi * r(j)^3 * simpsons(v(1:j),0,r(j),[]) - vol);
end

I am probably making a stupid mistake here, and I'm not too familiar with interpolation schemes in MATLAB, so any advice or help on either front is greatly appreciated.

Here is an example of one of the curves I am trying to integrate over. v is on the y-axis and r (fractional) is on the x-axis.

πŸ‘︎ 3
πŸ“°︎ r/matlab
πŸ’¬︎
πŸ‘€︎ u/ohshitgorillas
πŸ“…︎ Sep 22 2017
🚨︎ report
A Cute Little Analysis I did of Numerical Integration In Python. Thought You all Might Enjoy. [X-post r/DataIsBeautiful] github.com/Yurlungur/nume…
πŸ‘︎ 8
πŸ“°︎ r/Python
πŸ’¬︎
πŸ‘€︎ u/equationsofmotion
πŸ“…︎ May 19 2018
🚨︎ report
Numerical Integration for Rotational Dynamics

Article: http://cwzx.wordpress.com/2013/12/16/numerical-integration-for-rotational-dynamics/

This article may be of interest to game developers, particularly those that implement rotational motion, such as rigid bodies.

The article explains how the rotational equations of motion are more complicated than the translational versions (position, velocity, acceleration). The rotational equations are then solved under the assumption of constant angular acceleration. This allows for rotations and angular velocities to be advanced forwards in time by a time step. Unlike a generic integration scheme, this directly takes the rotational nature of the system into account, and the approximations that are made are explicitly stated.

Both rotation matrices and unit quaternion versions are included.

πŸ‘︎ 26
πŸ“°︎ r/gamedev
πŸ’¬︎
πŸ‘€︎ u/injx
πŸ“…︎ Dec 16 2013
🚨︎ report
Numerical integration of a perturbed ODE

So, I have an ODE (Actually there's an additional noise term but I don't think that's relevant), of the form

dx/dt = Ξ΅ f(x) + h(x) where Ξ΅ is not that large, but where f(x) is computationally hard to evaluate. If we take a timestep Ξ”, then an n-th order Runge Kutta method will give a local error of order Ξ”^(n+1) .

Since Ξ΅f(x) is small, an error in this term will give a smaller contribution to the total error than an error in the g(x) term. Will it give good results if we give the Ξ΅f(x) a larger time step, either directly (making it contribute only once every so often), or indirectly by using, say, 4th order Runge-Kutta for the g(x) but 1st order Runge-Kutta for f(x)?

I hope this is clear, thanks for reading!

πŸ‘︎ 3
πŸ“°︎ r/askmath
πŸ’¬︎
πŸ‘€︎ u/freemath
πŸ“…︎ Jan 30 2018
🚨︎ report
Programming Numerical Methods | The Trapezoidal Rule of Integration youtube.com/attribution_l…
πŸ‘︎ 2
πŸ“°︎ r/numerical
πŸ’¬︎
πŸ‘€︎ u/MechTutor
πŸ“…︎ Jan 17 2018
🚨︎ report
Anchor Physics - simulation of anchor + chain physics as an example of Verlet integration (numerical method), rendered on HTML5 Canvas pjama.github.io/projects/…
πŸ‘︎ 8
πŸ“°︎ r/javascript
πŸ’¬︎
πŸ‘€︎ u/magenta_placenta
πŸ“…︎ May 17 2017
🚨︎ report
Mathematics/Numerical Analysis: Estimating/Approximating an unknown function from a set of known integration results?

So here is the problem I've been thinking about. I apologize in advance, I don't know how to make mathematical markup on reddit. Please feel free to ask for clarifications if needed.

Statement of the Problem

  1. Let f(x) be some function with the following properties:
  • f(x) has finite support over some known interval, S.
  • f(x) is everywhere non negative.
  • f(x) is continuous with continuous derivatives, except maybe at the boundaries of the support.

&nbsp;

  1. Let P(x,y) be some known function which is integrable on the interval [0,infinity).

&nbsp;

  1. Let {Y_i} be a set of N values for which the result of the following operation is known:

&nbsp;

R(Y_i) = R_i = Integral[dx * f(x) * P(x,Y_i),{from 0 to +infinity}]

&nbsp;

  1. From the knowledge of {R_i} and P(x,y), find f(x) or some approximation of f(x).

&nbsp;

My solution:
I assume f(x) doesn't change too much over small enough intervals dx. I can then approximate f(x) as a series of N constant values over small intervals. Formally,

&nbsp;

f(x) = Sum[w_j * g_j(x),{j = 0 .. N}]

&nbsp;

Where w_j is the set of approximating constant values and g_j(x) is a set of functions which take the value 1 for the appropriate interval, j, and 0 everywhere else. The j intervals are simply equally spaced subdivisions of the finite support of f(x), S (though they don't have to be).

&nbsp;

I then have,

&nbsp;

R_i = Sum[w_j * P_i,j, {j = 1 .. N}]

&nbsp;

Where P_i,j is the result of the integration of the known function P(x,Y_i) over the j'th interval.

&nbsp;

The problem is now reduced to inverting the matrix P_i,j and calculating inv(P_i,j)*R_i to solve for the coefficients w_j. The resulting solution is adequate for my current needs but I do have some questions with an eye to future improvements:

&nbsp;

  1. I'm assuming this problem (or something similar/more general) has already been studied before. I'd appreciate any links to relevant material.
  2. Is there a way to choose the intervals g_j optimally?
  3. My solution results in a discontinuous approximation of f(x) which is known to be continuous. Is there a method find an approximation that is continuous?

&nbsp;

Any input would be much appreciated!
(edit: formatting)

πŸ‘︎ 14
πŸ“°︎ r/askscience
πŸ’¬︎
πŸ‘€︎ u/prblynot
πŸ“…︎ Mar 21 2015
🚨︎ report
ti-84 and numerical integration

I'm wondering if there isn't a way to add up sums in the ti-84 using a syntax similar to say where I have a function f(x) and can just enter f(0)+f(0.5)+f(1), etc. I know there is the summation key and there are also programs you can add to the calculator but specifically wondering if there's something like this already built in.

πŸ‘︎ 2
πŸ“°︎ r/calculus
πŸ’¬︎
πŸ‘€︎ u/alos87
πŸ“…︎ Oct 10 2017
🚨︎ report
Numerical Integration/derivation ON FPGA

I'm working on a project where I will be reading in a TDM data stream from an ADC attached to a sensor and doing some simple mathematical operations (Calculating RMS, peak detection, integration, derivatives, etc) on the data and displaying it to a user. This will all be done using the microzed eval kit.

I have two possible options:

  1. Read the TDM data stream into the PL, place each raw channel into a shared register, and have the PS access each raw sample, and do all of the processing in software...

  2. Read the TDM data stream into the PL, and have the FPGA do the numerical integration/derivative/etc on the fly, and store the raw sample along with the processed values into registers, leaving the PS open to do other things and access the latest data when it needs to...

I'm leaning towards the second option, which I think is a slicker solution. This way the PS can just peak at a register pertaining to the processed data when ever it needs it and won't have to waste time calculating the latest value. The latest processed value will be updated on the fly by the FPGA. I've never done any actual mathematical operations in VHDL, I'm a software guy. I wanted to know if anyone has any experience performing numerical integration, derivatives, etc in the fabric itself. If so, any helpful links/hints/example VHDL code? Are there ready made open source IPs for things like this? Is this more trouble than its worth?

πŸ‘︎ 6
πŸ“°︎ r/FPGA
πŸ’¬︎
πŸ‘€︎ u/l33tsk33t
πŸ“…︎ Mar 05 2015
🚨︎ report
Homework Help Numerical Integration

Need help with the second exercise of my intro to Matlab course here are the specific requirements: http://imgur.com/NWyxhpu

πŸ‘︎ 4
πŸ“°︎ r/matlab
πŸ’¬︎
πŸ‘€︎ u/bfitz2
πŸ“…︎ Apr 18 2016
🚨︎ report
Collected Algorithms of the ACM - From: Feb 1960, "obsolete, numerical integration, quadrature", title: QuadI, by: R.J. Herbold (); To: Nov 2012, Parallel solver for semidefinite programming problem having sparse Schur complement matrix, by: Yamashita, Fujisawa, Fukuda, Nakata & Nakata calgo.acm.org/
πŸ‘︎ 35
πŸ“°︎ r/programming
πŸ’¬︎
πŸ‘€︎ u/NevilleDNZ
πŸ“…︎ Feb 22 2013
🚨︎ report
Can someone help me debug my program that performs Numerical Integration? scratch.mit.edu/projects/…
πŸ‘︎ 3
πŸ“°︎ r/CasualMath
πŸ’¬︎
πŸ‘€︎ u/tmojeed
πŸ“…︎ Mar 10 2017
🚨︎ report
Need help with numerical integration (3 body problem)

This summer I am working on my matlab programming skills and I need some help with numeric integration. I am currently modeling a 3 body problem with the earth, moon and spacecraft. I have a script written with a function for the forces, accelerations, velocities and position. My issue is that I am not sure how to totally implement these equations correctly even though I have the written out. I will gladly post the code I have written so far so I can get some feedback on the next step to take.

πŸ‘︎ 2
πŸ“°︎ r/matlab
πŸ’¬︎
πŸ‘€︎ u/darthtumtum
πŸ“…︎ Jun 06 2017
🚨︎ report
Issue with numerical integration

I'm trying to do some numerical integration with numpy's trapz function, but I'm getting some peculiar behavior: the measured area under the curve seems to decrease when the slope is negative. I've got a symmetrical shape made up of ordered pairs and a function that plots and integrates it, and oddly enough it's reporting that the area is negative. Does anyone see the issue?

Pastebin to code: http://pastebin.com/38yEv39Y

System: Linux Mint, Python 3.5.2

πŸ‘︎ 2
πŸ“°︎ r/learnpython
πŸ’¬︎
πŸ‘€︎ u/jdr525
πŸ“…︎ Oct 07 2016
🚨︎ report
really need advice from people who use SciPy for numerical methods/ integration

Okay, so, I'm trying to use odeint correctly. I've gotten to the point where I can get solutions over time intervals consistently. I'm specifically using Python(x,y) with Cantera to model chemical kinetics and combustion processes.

My problem is that odeint often hangs my Python Interpreter regardless of my relative/ absolute tolerances, and regardless of my time steps.

Can anyone point me in the right direction on how to solve partial differential equations using odeint efficiently, and how to properly set my parameters?

πŸ‘︎ 5
πŸ“°︎ r/learnpython
πŸ’¬︎
πŸ‘€︎ u/casketballer
πŸ“…︎ Nov 13 2014
🚨︎ report
Quadrature – Numerical Integration Comparison austingwalters.com/quadra…
πŸ‘︎ 3
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/austingwalters
πŸ“…︎ Apr 02 2014
🚨︎ report
Good sources from which to learn about numerical integration?

A lot of you probably know a lot about numerical integration. If you do, could you provide the names of books or papers that one should read to get up to date on the techniques in the field?

πŸ‘︎ 3
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/zomglings
πŸ“…︎ May 26 2014
🚨︎ report
Lush (Lisp Universal Shell) for large-scale numerical and graphic applications with the easy integration of code written in C, C++, or other languages. lush.sourceforge.net/
πŸ‘︎ 10
πŸ“°︎ r/programming
πŸ’¬︎
πŸ‘€︎ u/dlyund
πŸ“…︎ Oct 20 2014
🚨︎ report
An introduction to numerical integration for beginners [PDF] ece.utah.edu/~ece6340/LEC…
πŸ‘︎ 25
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/roger_
πŸ“…︎ Feb 07 2011
🚨︎ report
[University Numerical Analysis] Error term within Numerical integration

http://i.imgur.com/SHnlHCM.png

^That is the question I'm trying to finish, I've done parts (a) and (b) and the first part of (c). The bit that I'm struggling with is the highlighted part.

So I've assumed the error term is of the form that was given. So I need to get a function f and take the fourth derivative. I'm thinking this will have to end up being a constant so that epsilon (sorry dont know the actual name of the greek letter) can disappear.

Would I just need to let f(x) = x^4 ? That just doesn't seem right to me..

Also have one more question, how would you use the trapezium rule to approximate an integral? E.g. integrate sin(x) from 0 to pi

When I apply the trapezium rule I get 0..

πŸ‘︎ 9
πŸ“°︎ r/learnmath
πŸ’¬︎
πŸ‘€︎ u/miserablestoic
πŸ“…︎ Nov 25 2015
🚨︎ report
Numerical methods for integration; log graph of the truncation error

I will use example to make my question clear:

My understand is that, this means:

  • error &lt;= C|h^3|

  • log error = 3 log(h) + C

  • if I draw the graph it should have slope 3

However, it has slope 2 instead... incidentally, there is something called the "order of the method" which in this case is equal to 2 (read beginning of page 44) as well.

  • Apparently, the slope of the graph is equal to "order of the method"

My question is:

  • Why is the slope of the log graph equal to (3 - 1) = 2? why do we deduct 1?
  • Why is this same number called the "order of the method"?

I posted this question on math.stackexchange and I have no idea what is wrong with that website.. after 24 hours I had only 7 views. This happened before and reddit helped me instead, so I hope this will happen again :)


I made a mistake in the measurement, I wasn't measuring the local error correctly. Now that it is corrected, it does give "3" and not "2".

  • still not sure if "order of the method" simply means the slope of the graph of the log of the global error.
  • still don't understand why the order of global error is +1 compared to the order of local error.
  • measuring the local error for multi-step method gives strange results, I will have to look into that. Currently what I'm doing is simply abs(true value - method_approximation_based_on_previous_true_values)
πŸ‘︎ 2
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/incraved
πŸ“…︎ Aug 30 2014
🚨︎ report
"Cumulative trapezoidal numerical integration," also known as... mathworks.com/access/help…
πŸ‘︎ 23
πŸ“°︎ r/funny
πŸ’¬︎
πŸ‘€︎ u/voidoid
πŸ“…︎ Jun 10 2009
🚨︎ report
Help with a numerical integration problem

Hi, I have a problem we have been given in an assignment for a class. A graph of output over time is given to us, along with the exact values, 20 values a second for 30 seconds. I understand the concept of the numerical integration - finding the area for each split and finding the sum of all, but i cant figure out how to add together the areas.

How can I take the values i need for one calculation, and then add it to the area for the second one and so on in VBA?

Sorry for what I imagine is a very basic thing, we have only been doing this a couple weeks and my notes are no help to me, as all they do is explain the theory of adding the areas together, and dont actually show us any of the code. Been googling for quite a while now and cant find anything that makes sense to me either.

Thanks

πŸ‘︎ 2
πŸ“°︎ r/vba
πŸ’¬︎
πŸ‘€︎ u/thelazyfool
πŸ“…︎ Feb 21 2015
🚨︎ report
Numerical integration?

I have an vector of position values. If I take the diff(x) to get dx/dt, I see that there are a few outliers (The data is supposed to have a constant velocity). If I correct these values, is it possible to integrate to get the approximately correct vector x?

πŸ‘︎ 4
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/j_lyf
πŸ“…︎ May 14 2011
🚨︎ report
Numerical Integration In Excel Using The Trapezoidal Rule myengineeringworld.net/20…
πŸ‘︎ 3
πŸ“°︎ r/excel
πŸ’¬︎
πŸ‘€︎ u/Jordna33
πŸ“…︎ Jun 18 2013
🚨︎ report
Quadrature – Numerical Integration Comparison austingwalters.com/quadra…
πŸ‘︎ 2
πŸ“°︎ r/statistics
πŸ’¬︎
πŸ‘€︎ u/austingwalters
πŸ“…︎ May 16 2014
🚨︎ report
newbie help: understanding curried functions in the context of numerical integration using hmatrix

Greetings,

I am working on my first Haskell project, and in this project I need to numerically evaluate some 2D integrals over rectangular regions. I am using hmatrix which interfaces with some numerical integration routines from GSL.

To start, I would like to replicate the double integral example from Wikipedia.

let f x y = x^2 + 4 * y

Now, my understanding is that I want the inner integral to be evaluated first, which should return me a function that I can then integrate over the outer bounds of integration (y = 7 to 10). However, I am having trouble translating this into Haskell. There is an example on p.19 of the hmatrix pdf, but I am not seeing how to apply that to this problem.

Help would be very much appreciated!

πŸ‘︎ 2
πŸ“°︎ r/haskell
πŸ’¬︎
πŸ‘€︎ u/stochasticMath
πŸ“…︎ Jun 28 2013
🚨︎ report
Cumulative trapezoidal numerical integration mathworks.com/help/techdo…
πŸ‘︎ 9
πŸ“°︎ r/funny
πŸ’¬︎
πŸ‘€︎ u/jweebo
πŸ“…︎ Dec 10 2010
🚨︎ report
What's generally more accurate for numerical integration, the Mid-point rule or the Trapezium rule?
πŸ‘︎ 3
πŸ“°︎ r/math
πŸ’¬︎
πŸ‘€︎ u/Fierytoadfriend
πŸ“…︎ May 11 2018
🚨︎ report

Please note that this site uses cookies to personalise content and adverts, to provide social media features, and to analyse web traffic. Click here for more information.