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


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 
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  
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  
Educational Games

Find the Difference Game 1.00.60
In this game, children find 12 differences between two pictures. Something may be a different color, something may be missing, something may be added, something may be a different size, or something may be changed in some other way, but with good... 
3.71 MB  
Games

Harmonic Convergence for Mac OS 1.0
Enter a world of beautiful graphics, hypnotic music, and addictive gameplay.Harmonic Convergence Deluxe blends these elements together into an unforgettable experience that will have you losing track of time trying to outdo your last top score,... 
18 MB  
Utilities

Difference Between Two Dates Software 7.0
This software offers a solution for users who want to find the difference between two dates. The user simply selects a date from the top calendar and one from the lower calendar and with one click the difference will be displayed in four lines:... 
396 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

Finite Difference Laplace Equation Solver using unequal square grid xy grids 1.0
The attachment contains:1. Related matlab files.2. Picture files of possible outputs.3. A read me text file.4. A report containing detailed explanations about the basics and about coding algorithm used herein.Finite difference (central) method is... 
15.8 MB  
Development Tools

temperature analysis using finite difference method 1.0
Study of heat transfer and temperature of a 1x1 metal plateheat is dissipated through the left right and bottom sides and emp at infinity is t n(n1) points in consideration,Temperature at top end is 500*sin(((i1)*pi)/(n1)A*Temp=U where A is... 
10 KB  
Development Tools

Babbage Difference Engine Emulator 1.0
Emulate Babbage Difference Engine 2 at Mountain View Computer History Museum Program to emulate evaluation of 7th order polynomial with 31 digit accuracy for Difference Engine designed by Charles Babbage around 1830. Two versions of the engine... 
1.03 MB  
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 
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 