Code Directory
 Visual Basic & VB.NET
New Code
dbForge Studio for PostgreSQL 2.3.212
HTMLPad 2020 16.2
WeBuilder 2020 16.2
Rapid CSS 2020 16.2
Rapid PHP 2020 16.2
C# HTML to PDF 2020.8.1
Vue Injector 3.3
Spectrum Analyzer pro Live 2019
Devart Excel Add-in for HubSpot 2.1
RentALLScript - Airbnb clone 2.2
SuiteCRM Theme Customization 7.11.6
iScripts NetMenus 3.1
iScripts EasyIndex 2.2
iScripts EasySnaps 2.0
Top Code
IcrediBB Bulletin Board System 1.0
Ruby FANN-Fast Artificial Neural Network 1.1.3
Java-2-Pseudo 1.0
Ticket Booking System 1.0
Simplified Fuzzy ARTMAP Neural Network 1.0
Block Matching Algorithms for Motion Estimation 1.0
dbForge Studio for PostgreSQL 2.3.212
Online Betting Exchange & Open API. 1.0
PHP reports generator 1.0
SkaDate Online Dating Software 4.0
Betting system 6.x-1.x-dev
Kalman Filter Tutorial 1.0
Gif4J Light 1.0
Solid File System OS edition 5.1
Auto Screen Capture and Word / PDFDumper 1.0
Top Rated
Uber Clone with Safety Measure Addons 2.0
Answers phpSoftPro 3.12
phpEnter 5.1.
Quick Maps For Dynamics CRM 3.1
Single Leg MLM 1.2.1
Azizi search engine script PHP 4.1.10
Paste phpSoftPro 1.4.1
Extreme Injector 3.7
Apphitect Airbnb Clone Script 1.0
Deals and Discounts Website Script 1.0.2
Solid File System OS edition 5.1
Classified Ad Lister 1.0
Aglowsoft SQL Query Tools 8.2
Invoice Manager by PHPJabbers 3.0
ICPennyBid Penny Auction Script 4.0
Predictor-corrector PECE method for fractional differential equations 1.0
File ID: 81018

Predictor-corrector PECE method for fractional differential equations 1.0
Download Predictor-corrector PECE method for fractional differential equations 1.0http://www.mathworks.comReport Error Link
License: Shareware
File Size: 10.0 KB
Downloads: 45
Submit Rating:
Predictor-corrector PECE method for fractional differential equations 1.0 Description
Description: FDE12 solves an initial value problem for a non-linear differential equation of fractional order (FDE).

This is an implementation of the predictor-corrector method of Adams-Bashforth-Moulton described in [1]. Convergence and accuracy of the method are studied in [2]. The implementation with multiple corrector iterations has been proposed and discussed for multiterm FDEs in [3]. In this implementation the discrete convolutions are evaluated by means of the FFT algorithm described in [4] allowing to keep the computational cost proportional to N*log(N)^2 instead of N^2 as in the classical implementation; N is the number of time-point in which the solution is evaluated, i.e. N = (TFINAL-T)/H. The stability properties of the method implemented by FDE12 have been studied in [5].


[T,Y] = FDE12(ALPHA,FDEFUN,T0,TFINAL,Y0,h) integrates the initial value problem for the FDE, or the system of FDEs, of order ALPHA > 0
Y^(k)(T0) = Y0(:,k+1), k=0,...,m-1
where m is the smallest integer grater than ALPHA and D^ALPHA is the fractional derivative according to the Caputo's definition. FDEFUN is a function handle corresponding to the vector field of the FDE and for a scalar T and a vector Y, FDEFUN(T,Y) must return a column vector. The set of initial conditions Y0 is a matrix with a number of rows equal to the size of the problem (hence equal to the number of rows of the output of FDEFUN) and a number of columns depending on ALPHA and given by m. The step-size H>0 is assumed constant throughout the integration.

[T,Y] = FDE12(ALPHA,FDEFUN,T0,TFINAL,Y0,H,PARAM) solves as above with the additional set of parameters for the FDEFUN as FDEFUN(T,Y,PARAM).

[T,Y] = FDE12(ALPHA,FDEFUN,T0,TFINAL,Y0,H,PARAM,MU) solves the FDE with the selected number MU of multiple corrector iterations. The following values for MU are admissible:
MU = 0 : the corrector is not evaluated and the solution is provided just by the predictor method (the first order rectangular rule);
MU > 0 : the corrector is evaluated by the selected number MU of times; the classical PECE method is obtained for MU=1;
MU = Inf : the corrector is evaluated for a certain number of times until convergence of the iterations is reached (for convergence the difference between two consecutive iterates is tested).
The defalut value for MU is 1

[T,Y] = FDE12(ALPHA,FDEFUN,T0,TFINAL,Y0,H,PARAM,MU,MU_TOL) allows to specify the tolerance for testing convergence when MU = Inf. If not specified, the default value MU_TOL = 1.E-6 is used.


[1] K. Diethelm, A.D. Freed, The Frac PECE subroutine for the numerical solution of differential equations of fractional order, in: S. Heinzel, T. Plesser (Eds.), Forschung und Wissenschaftliches Rechnen 1998, Gessellschaft fur Wissenschaftliche Datenverarbeitung, Gottingen, 1999, pp. 57-71.

[2] K. Diethelm, N.J. Ford, A.D. Freed, Detailed error analysis for a fractional Adams method, Numer. Algorithms 36 (1) (2004) 31-52.

[3] K. Diethelm, Efficient solution of multi-term fractional differential equations using P(EC)mE methods, Computing 71 (2003), pp. 305-319.

[4] E. Hairer, C. Lubich, M. Schlichte, Fast numerical solution of nonlinear Volterra convolution equations, SIAM J. Sci. Statist. Comput. 6 (3) (1985) 532-541.

[5] R. Garrappa, On linear stability of predictor-corrector algorithms for fractional differential equations, Internat. J. Comput. Math. 87 (10) (2010) 2281-2290.

License: Shareware

Related: consecutive, iterates, Difference, reached, convergence, tested, defalut

O/S:BSD, Linux, Solaris, Mac OS X

File Size: 10.0 KB

Downloads: 45

More Similar Code

Durand-Kerner method for solving polynomial equations.

This submission contains the basic functions that are necessary for using the matrix approach to discretization of fractional differential equations, and demos.

The method is described in the following articles:
[1] I. Podlubny,...

the function is GUI that simulate double pendulum by using 4th-order Runge-Kutta algorithm for the differential equations.
the differential equation are very similar to eq. (12) and (13) here:

This demo shows examples of using the MATLAB Symbolic Math Toolbox for Ordinary Differential Equations (ODE). The examples range from relatively basic to rather advanced, and are placed according to the topics, as traditionally covered in...

cuLsoda is an implementation of the Livermore Solver for Ordinary Differential Equations with Automatic switching for stiff and non-stiff problems (LSODA), written in NVIDIA's C for CUDA GPGPU API.

fractPC.m provides the Predictor-corrector Method for fractional order Relaxation-Oscillation equation (u_t_alpha=-B u(t)+f(t), 0<alpha<=2).

variablerandomorder.m provides the Predictor-corrector Method for the variable-order...

Graphical user interface (GUI) is used to solve up to two ordinary differential equations (ODEs). Results can be plotted easily. Choose between MATLAB's ode45 (non-stiff solver) or ode15s (stiffer solver).

This is primarily a teaching...

Halley's method for solving equations of type f(x)=0.

Solves the linear least squares problem with nonnegative variables using the predictor-corrector algorithm in:
Portugal, Judice and Vicente, A comparison of block pivoting and interior point algorithms for linear least squares problems with...

Displaying the system of differential equations for the spatial model of the car and shows you how to make using the software MATLAB / Simulink.

User Review for Predictor-corrector PECE method for fractional differential equations
- required fields

Please enter text on the image