Download Shareware and Freeware Software for Windows, Linux, Macintosh, PDA

line Home  |  About Us  |  Link To Us  |  FAQ  |  Contact

Serving Software Downloads in 976 Categories, Downloaded 29.541.181 Times

Predictor-corrector PECE method for fractional differential equations 1.0

  Date Added: September 01, 2013  |  Visits: 400

Predictor-corrector PECE method for fractional differential equations

Report Broken Link
Printer Friendly Version

Product Homepage
Download (29 downloads)

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].USAGE:[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 D^ALPHA Y(t) = FDEFUN(T,Y(T)) Y^(k)(T0) = Y0(:,k+1), k=0,...,m-1where 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.REFERENCES:[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.

Requirements: No special requirements
Platforms: Matlab
Keyword: Consecutive Convergence Defalut Difference Iterates Reached Tested
Users rating: 0/10

License: Shareware Size: 10 KB
More Reviews or Write Review

Development Tools  -  maroce 1.0
This m-file deals with the iterated principal factor method (principal axis factoring) thru the complete data matrix, the latent root criterion, iteration convergence criterion, and uses the varimax factor rotation. It works with an iterative...
10 KB  
3D Graphic Tools  -  HDR Visual Difference Predictor 1.6
Visual difference metrics can predict whether differences between two images are visible to the human observer or not. Such metrics are used for testing either visibility of information (whether we can see important visual information) or...
133.12 KB  
Communication  -  Binary Conversion Methods 1.0
Five methods for converting a signal to a binary representation. Methods outlined in Arthur Petrosian, "Kolmogorov Complexity of finite sequences and recognition of different preictal eeg patterns".Methods are:1. average method: sample assigned 1...
10 KB  
Development Tools  -  fminsearch2 1.0
This is a (small) change to the current fminsearch routine to improve convergence.The current fminsearch selects n additional simplex points based on a small change in each variable. This implementation uses an adaptive procedure to select these n...
20.48 KB  
Development Tools  -  MCSpreadOption.m 1.0
The MCSpreadOption function can be used for valuation of spread options Spread options are actively traded on New York Mercantile Exchange (NYMEX) on the difference/spread between different oil qualities. Spread options are also popular in many...
10 KB  
Development Tools  -  marocewod 1.0
This m-file deals with the iterated principal factor method (principal axis factoring) thru the R-correlation matrix (without the matrix of data), the latent root criterion, iteration convergence criterion, and uses the varimax factor rotation. It...
10 KB  
Network & Internet  -  django-mongokit 0.2.2
django-mongokit was built for Django 1.2 with the multi-db support but you can use django-mongokit in Django 1.1 (tested in Django 1.1.1) as a secondary database. For example, you might want to continue running your application in...
10.24 KB  
Games  -  Avoidance Dynamite 1.4.0
Avoidance Dynamite (iOS 5 tested) Reached No.1 in New & Noteworthy! 4 **** review (by For iPhone, iPod touch & iPad (HD) *** AROUND THE WORLD *** Being played in over 85 countries and...
31.7 MB  
Games  -  Cricket Sudoku
The time tested Sudoku game with a difference. Enjoy the challenge of a general Sudoku game but this time not with numbers. Requirements: Windows Phone 8.1, Windows Phone 8, Windows Phone 7.5
3 MB  
WYSIWYG HTML Editors  -  Author.CMS | Usability-tested Content Management System 3.3
Author CMS is a usability-tested content management system. The package offers a host of remarkable features such as: · Usability-tested administration interface · Special installer that easily installs the system on a server · High-level...
Scripts  -  Freelancer Script 5.05
Main Features: 100% Secured. Email Support (3 Years). FREE Updates (3 Years). Post projects. Featured projects. Private projects. Sealed projects. Edit/delete projects. Select freelancers....
5.49 MB  
Scripts  -  B2B Script 4.20
Main Features: 100% Secured. Email Support (3 Years). FREE Updates (3 Years). Sign-up Account (Registration of account). Lead generation tools (for the sellers). Email verification to...
5.49 MB  
Scripts  -  B2C Script 5.06
Main Features: 100% Secured. Email Support (3 Years). FREE Updates (3 Years). The script comes with totally editable site colors, icons and graphics Multilevel categories allows extensive browsing Admin can change Category ordering or...
5.49 MB  
Scripts  -  Social Networking Script 2.86
Main Features: 100% Secured. Email Support (3 Years). FREE Updates (3 Years). Registration with name, email, password, date of birth etc. User can add multiple school, college, university with start...
5.49 MB  
Scripts  -  Business Networking Script 8.04
Main Features: 100% Secured. Email Support (3 Years). FREE Updates (3 Years). Ajax based interface. Profile creation. Different types of profile. Profile for jobseekers, employers and employed...
5.49 MB  
Development Tools  -  VMP Viewer 1.0
This is a very rudimentary tool to visualize the VMP files generated by BrainVoyager. Useful to share files with people who do not have BV.
10 KB  
Development Tools  -  Sending reports and timestamped file by emailing 1.0
main executing reference usage:[1] usage_send_mail.mIllustrates email sending with multiple separate files or single timestamped tar file. Attachment failure is properly handled, with continuation of report emailing without the attachment.[2]...
768 KB  
Development Tools  -  IrisMVC 2.0 rc1
IrisMVC is an OOP PHP framework that developers can use as a strong and secure foundation to build on various web applications following the Model-View-Controller (MVC) pattern. It provides the basic functionality developers need, without...
51.2 KB  
Development Tools  -  7-Zip for Script 4.42
7-Zip is a file archiver with a high compression ratio.Features:- High compression ratio in new 7z format with LZMA compression- Supported formats:- Packing / unpacking: 7z, ZIP, GZIP, BZIP2 and TAR- Unpacking only: RAR, CAB, ISO, ARJ, LZH, CHM,...
624.64 KB  
Development Tools  -  Barnardextest 1.0
This file, as the Fisher's exact test, performs the exact probability test for a table of frequency data cross-classified according to two categorical variables, each of which has two levels or subcategories (2x2). It is a non-parametric...
10 KB