                                                 {
 //      Take a step using RK3
 
 //      Initialization
         Y0 = Y;
 
 //      Step 1
         (*F)(Y0,Yh); Yh *= h;                         // Yh = h*F(Y0)
         Y0 += Yh;
 //      Y0 = Y0 + h*Yh
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
 //      Step 2
         (*F)(Y0,Yh); Yh *= h;                      // Yh = h*F(Y0)
         Y0 += Yh; Y *= 3.0; Y0 += Y;  Y0 *= 0.25;  // Changed Y to 3*Y!
 //      Y0 = 1/4 * ( 3*Y + (Y0 + h*Yh) )
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
 //      Step 3
         (*F)(Y0,Yh); Yh *= h;                        // Yh = h*F(Y0)
         Y0 += Yh; Y0 *= 6.0; Y += Y0; Y *= (1.0/9.0);
 //      Y  = 1/3 * ( Y + 2 * (Y0+h*Yh) )
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 
 
         return Y;
     }

◆ operator()() [2/2]

template<class T >

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

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

                                                             {
 //      Take a step using RK3
 
 //      Initialization
         Y0 = Y;
 
 //      Step 1
         (*F)(Y0,Yh,dir); Yh *= h;                         // Yh = h*F(Y0)
         Y0 += Yh;
 //      Y0 = Y0 + h*Yh
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
 //      Step 2
         (*F)(Y0,Yh,dir); Yh *= h;                      // Yh = h*F(Y0)
         Y0 += Yh; Y *= 3.0; Y0 += Y;  Y0 *= 0.25;  // Changed Y to 3*Y!
 //      Y0 = 1/4 * ( 3*Y + (Y0 + h*Yh) )
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
 //      Step 3
         (*F)(Y0,Yh,dir); Yh *= h;                        // Yh = h*F(Y0)
         Y0 += Yh; Y0 *= 6.0; Y += Y0; Y *= (1.0/9.0);
 //      Y  = 1/3 * ( Y + 2 * (Y0+h*Yh) )
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 
 
         return Y;
     }