Calculators and numerical differentiation

11012020, 05:39 PM
Post: #5




RE: Calculators and numerical differentiation
(11012020 05:39 AM)Wes Loewer Wrote: Instead of averaging the backwards and forwards differences, why not just check to see if f(a) exists and if it does then do the central difference method? Try taken derivative of f(x) = x, at x=0. f(0) = 0 = 0, thus exist. Central difference slope = (0+h  0h)/(2h) = 0 Forward difference slope = (h0)/h = 1 Backward difference slope = (0h)/h = 1 2 slopes does not match. f'(0) does not exist. Oh ... Casio does not check this, d/dx(√(X²), 0) → 0 Quote:The nonCAS Nspire apparently has a bit of CAS hidden under the hood because it does not use this approximation. It appears to evaluate the derivative symbolically and then evaluates that expression numerically, keeping the CAS carefully hidden from the user. How did you deduce there is hidden CAS under nonCAS Nspire ? May be an example ? Quote:Had anyone previously seen the del operator used for the backwards difference as shown in the Casio manual? I think many text use the same symbols (or equivalent): Ref: Fundamentals of Numerical Analysis, by Stephen Kellison Code: Operators: Or, operator notation: \(Δ=E1,\quad ∇=1E^{1}\) Combining µδ, we have the 3 points central difference slope formula \( hD ≈ µδ = \large\left({E^{½} + E^{½} \over 2}\right) \normalsize(E^{½}E^{½}) = \large{E\;\;E^{1} \over 2} = \large{Δ\;+\;∇ \over 2} \) For more accuracy, we can add more terms: (note, there is no even powers of δ ) \( hD ≡ µδ \large\left(1  {δ^2\over 6} + {δ^4\over 30}  {δ^6\over 140} + {δ^8\over 630}  {δ^{10}\over 2772} + {δ^{12}\over 12012}  {δ^{14}\over 51480} \;+\; ... \right) \) (10302020 11:41 PM)Paul Dale Wrote: The order 10 methods are a weighted summation of \( f(x \pm 1) \), \( f(x \pm 2) \), \( f(x \pm 3) \), \( f(x \pm 4) \) and \( f(x \pm 5) \). Lets build D (order 10), without using central difference table. XCas> c := E^(1/2)  E^(1/2); // central difference XCas> m := (E^(1/2) + E^(1/2))/2; // mean XCas> [mc, cc] := expand(simplify([m,c] .* c)) → [E/2  1/(E*2) , E  2 + 1/E] Note the symmetry of m*c, c*c. This suggested D also have similar symmetry XCas> simplify(mc/h * horner([1/630, 1/140, 1/30, 1/6, 1], cc)) \(→ D = \large \frac {2100(EE^{1})\; \;600(E^2E^{2})\; +\;150(E^3E^{3})\; \;25(E^4E^{4})\; +\;2(E^5E^{5})} {2520h}\) 

« Next Oldest  Next Newest »

Messages In This Thread 
Calculators and numerical differentiation  robve  10302020, 09:57 PM
RE: Calculators and numerical differentiation  Paul Dale  10302020, 11:41 PM
RE: Calculators and numerical differentiation  Albert Chan  10312020, 01:20 AM
RE: Calculators and numerical differentiation  Wes Loewer  11012020, 05:39 AM
RE: Calculators and numerical differentiation  Albert Chan  11012020 05:39 PM
RE: Calculators and numerical differentiation  Albert Chan  11012020, 11:43 PM
RE: Calculators and numerical differentiation  Wes Loewer  11032020, 06:09 PM
RE: Calculators and numerical differentiation  Albert Chan  11032020, 10:14 PM
RE: Calculators and numerical differentiation  Wes Loewer  11042020, 04:14 PM
RE: Calculators and numerical differentiation  CMarangon  11032020, 06:55 PM
RE: Calculators and numerical differentiation  Wes Loewer  11042020, 04:04 PM

User(s) browsing this thread: 1 Guest(s)