backward finite difference with fractional differential equat. (GL definition)by sedo_5 Tags: finite difference, fractional calculus, grunwaldletnikov, mathematica 

#1
Sep2111, 07:11 PM

P: 1

Hi everyone!
I'm having a problem with backward finite difference method. I'm dealing with this problem using Mathematica software (and I'm going to paste some code below, maybe you don't exactly know the syntax, but it's pretty easy I think).. I need to use the backward finite difference method with the GrunwaldLetnikov fractional derivative definition (probably it's an uncommon topic..). Long story short, GL definition for the fractional derivative is this(two links here): http://en.wikipedia.org/wiki/Gr%C3%B...kov_derivative and http://www.xuru.org/fc/GrunwaldLetnikov.asp and can be usually used for backw.FD in order to approximate some fractional differential equation. This is my problem. I have a first order massspring differetial equation and and I have to write a code in Mathematica to approximate it (with backwFD method) and then compare it with the fractional differential equation case. (ps: the same problem appears in the second order equation (massspringdamper) approximation, but let's just limit this on the first order). The problem is that, for the ordinary case I have no problem in using FD method, instead I have a big problem in using FD method for the fractional case: I'll write some code to be more precise.  first order  ordinary 
 first order  fractional  everything as before, but with this FD cicle
 The first code (ordinary case) works flawlessly, but the second doesn't: the curve is completely different depending on the "deltah" imposed (it gives no problems when deltah=1 , but smaller values are generally used in approximations...and that's where this code gives problems). Now, i've tried quite everything with this code but to no avail...the problem is for sure given by deltah but i cannot find a way. It would be awesome if you could give me some hints on this. I think you could give me some hints also if you don't know Mathematica. PS: Since it's my first post here (by the way...this forum is awesome and I say Hi to everybody :) ), please excuse me if I'm posting it in the wrong section or if I'm going against some netiquette. Thanks a lot in advance. R. 


Register to reply 
Related Discussions  
Finite Difference Discretization of a Fourth Order Partial Differential Term  Differential Equations  4  
Difference between Forward and Backward Fourier Transforms?  Calculus  1  
Numerical differentiation using forward, backward and central finite difference  Differential Equations  2  
System of ODE Boundary Value Problem with 2nd Order Backward Difference  Differential Equations  1  
Definition of Fractional Occupation  Classical Physics  0 