JMathLabTutorial:Numerical Integration

From jWork.ORG
Revision as of 20:31, 27 December 2016 by Jworkorg (talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search
Limitted access. Reguest membership or login to this link first if you are already a member
Contents


Numerical Integration

Two routines are supplied for numerical integration, which are both not very sophisticated. quad('expression',ll,ul) is modeled after the Octave/Matlab integration function, but much simpler. Simpson's method is applied with a fixed number of nodes. This function uses the "eval"-method rather than symbolic variables. The function has to be supplied as quoted string, and must be compatible with vector arguments. Finally, the variable name must be x. romberg uses symbolic function definitions, and a symbolic variable has to be supplied. The maximum number of iterations is set by the variable rombergit (default $11$) and accuracy by rombergtol (default: [math]10^{-4}[/math]).


syms x
y=quad('exp(-x.^2)',0,5)
printf('%f',y)


syms x
y=romberg(exp(-x^2),x,0,5)
printf('%f',y)