M a t l a b   f i l e s   i m p l e m e n t i n g   t h e   a d a p t i v e   d e n s i t y   p r o p a g a t i o n   p r o c e d u r e   ( A D P )   a s   d e s c r i b e d   i n 
 
 A n d r e a   Y .   W e i s s e   &   W i l h e l m   H u i s i n g a ,   " E r r o r - c o n t r o l l e d   g l o b a l   s e n s i t i v i t y   a n a l y s i s   o f   o r d i n a r y   d i f f e r e n t i a l   e q u a t i o n s " ,   J o u r n a l   o f   C o m p u t a t i o n a l   P h y s i c s ,   2 0 1 1 
   
 w r i t t e n   b y 
 A n d r e a   Y .   W e i s s e 
 C e n t r e   f o r   S y s t e m s   B i o l o g y   a t   E d i n b u r g h 
 U n i v e r s i t y   o f   E d i n b u r g h 
 
 e m a i l :   a n d r e a . w e i s s e @ e d . a c . u k 
 
 C o p y r i g h t   ( C )   2 0 1 1 ,   U n i v e r s i t y   o f   E d i n b u r g h 
 
 F O R   A C A D E M I C   U S E   t h i s   p r o g r a m   i s   f r e e   s o f t w a r e ;   y o u   c a n   r e d i s t r i b u t e   
 i t   a n d / o r   m o d i f y   i t   u n d e r   t h e   t e r m s   o f   t h e   G N U   G e n e r a l   P u b l i c   L i c e n s e 
 a s   p u b l i s h e d   b y   t h e   F r e e   S o f t w a r e   F o u n d a t i o n ;   e i t h e r   v e r s i o n   2 
 o f   t h e   L i c e n s e ,   o r   ( a t   y o u r   o p t i o n )   a n y   l a t e r   v e r s i o n .   
 S e e   h t t p : / / w w w . g n u . o r g / l i c e n s e s / g p l . h t m l   f o r   d e t a i l s . 
 
 T h i s   p r o g r a m   i s   d i s t r i b u t e d   i n   t h e   h o p e   t h a t   i t   w i l l   b e   u s e f u l , 
 b u t   W I T H O U T   A N Y   W A R R A N T Y ;   w i t h o u t   e v e n   t h e   i m p l i e d   w a r r a n t y   o f 
 M E R C H A N T A B I L I T Y   o r   F I T N E S S   F O R   A   P A R T I C U L A R   P U R P O S E .     S e e   t h e 
 G N U   G e n e r a l   P u b l i c   L i c e n s e   f o r   m o r e   d e t a i l s .   
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 F i l e s   i n c l u d e d : 
 A D P . m 	 	 	 c t a u A . m 	 	 	 h i l l 1 d _ F . m 	 	 h i l l 2 d _ D F . m 	 	 h i l l 2 d _ t r a c e D F . m 	 s c a l a r P r o d . m 
 a d j a c e n c y M a t r i x . m 	 h i l l 1 d _ D F . m 	 	 h i l l 1 d _ t r a c e D F . m 	 h i l l 2 d _ F . m 	 	 p r e c o m p u t a t i o n . m 	 t i m e I n t e g r a t i o n . m 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 U S A G E : 
 
 T o   c o m p u t e   t h e   e x a m p l e s   s h o w n   i n   t h e   a r t i c l e   c a l l   t h e   s c r i p t   A D P . m ,   w h i c h   i n i t i a l i s e s   a l l   m o d e l   s p e c i f i c   a n d   n u m e r i c a l   p a r a m e t e r s   a n d   c a l l s   t h e   t i m e   i n t e g r a t i o n   r o u t i n e   i m p l e m e n t e d   i n   t h e   s c r i p t   t i m e I n t e g r a t i o n . m . 
 
 T h e   m o d e l   c a n   b e   s w i t c h e d   b e t w e e n   t w o   e x a m p l e s   ' h i l l 1 d '   a n d   ' h i l l 2 d '   b y   c h a n g i n g   t h e   v a r i a b l e   M O D E L   i n   A D P . m 
 
 T h e   f i l e s   
 h i l l 1 d _ F . m 	 h i l l 1 d _ D F . m 	 h i l l 1 d _ t r a c e D F . m 
 h i l l 2 d _ F . m 	 h i l l 2 d _ D F . m 	 h i l l 2 d _ t r a c e D F . m 
 i m p l e m e n t   t h e   R H S   F ,   t h e   J a c o b i a n   o f   F   a n d   t h e   t r a c e   o f   t h e   J a c o b i a n ,   r e s p e c t i v e l y ,   f o r   t h e   ' h i l l 1 d '   a n d   ' h i l l 2 d '   m o d e l s . 
 
   
 
 