                                                 {
 //      Take a step using RK2
 
 //      Initialization
         Y0 = Y;
 
 //      Step 1
         (*F)(Y0,Yh); Yh *= h;     // Yh = h*F(Y0)
         Y0 += Yh;                 // Y0 = Y0 + h*Yh
         Yh *= 0.5; Y  += Yh;      // Y  = Y  + (h/2)*F(Y0)      
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
 //      Step 2
         (*F)(Y0,Yh); Yh *= 0.5*h; // Yh = (h/2)*F(Y0)
         Y += Yh;                  // Y = Y + (h/2)*F(Y0)
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
         return Y;
     }

◆ operator()() [2/2]

template<class T >

T & Algorithms::RK2< T >::operator()	(	T &	Y,
		double	h,
		AbstFunctor< T > *	F,
		size_t	dir
	)

Definition at line 524 of file lib-algorithms.h.

                                                             {
 //      Take a step using RK2
 
 //      Initialization
         Y0 = Y;
 
 //      Step 1
         (*F)(Y0,Yh,dir); Yh *= h;     // Yh = h*F(Y0)
         Y0 += Yh;                 // Y0 = Y0 + h*Yh
         Yh *= 0.5; Y  += Yh;      // Y  = Y  + (h/2)*F(Y0)      
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
 //      Step 2
         (*F)(Y0,Yh,dir); Yh *= 0.5*h; // Yh = (h/2)*F(Y0)
         Y += Yh;                  // Y = Y + (h/2)*F(Y0)
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
         return Y;
     }