Predictorcorrector PECE method for fractional differential equations 1.0 
Date Added: September 01, 2013  Visits: 431 


FDE12 solves an initial value problem for a nonlinear differential equation of fractional order (FDE).This is an implementation of the predictorcorrector method of AdamsBashforthMoulton 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 timepoint in which the solution is evaluated, i.e. N = (TFINALT)/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,...,m1where 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 stepsize 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.E6 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. 5771.[2] K. Diethelm, N.J. Ford, A.D. Freed, Detailed error analysis for a fractional Adams method, Numer. Algorithms 36 (1) (2004) 3152.[3] K. Diethelm, Efficient solution of multiterm fractional differential equations using P(EC)mE methods, Computing 71 (2003), pp. 305319.[4] E. Hairer, C. Lubich, M. Schlichte, Fast numerical solution of nonlinear Volterra convolution equations, SIAM J. Sci. Statist. Comput. 6 (3) (1985) 532541.[5] R. Garrappa, On linear stability of predictorcorrector algorithms for fractional differential equations, Internat. J. Comput. Math. 87 (10) (2010) 22812290.

License: Shareware  Size: 10 KB 
Development Tools

maroce 1.0
This mfile 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 mfile deals with the iterated principal factor method (principal axis factoring) thru the Rcorrelation matrix (without the matrix of data), the latent root criterion, iteration convergence criterion, and uses the varimax factor rotation. It... 
10 KB  
Network & Internet

djangomongokit 0.2.2
djangomongokit was built for Django 1.2 with the multidb support but you can use djangomongokit 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 WhatsOniPhone.com) For iPhone, iPod touch & iPad (HD) *** AROUND THE WORLD *** Being played in over 85 countries and... 
31.7 MB  
Games

Cricket Sudoku 1.0.0.0
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  
Spreadsheets

Portfolio Optimization 5.0
The Portfolio Optimization template identifies the optimal capital weightings for a portfolio of financial investments that gives the desired risk and return profile based on the correlation between individual investments. The design of the... 
544 KB 
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). Signup 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 ModelViewController (MVC) pattern. It provides the basic functionality developers need, without... 
51.2 KB  
Development Tools

7Zip for Script 4.42
7Zip 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 crossclassified according to two categorical variables, each of which has two levels or subcategories (2x2). It is a nonparametric... 
10 KB 