“…We compute the values of 𝑦 n , 𝑦 n+1 , 𝑦 ′ n , and 𝑦 ′ n+1 and introduce them in the equations of the method ( 5)- (7). If we use Euler's identity, e iv = cos(𝜐) + i sin(𝜐), and compare the real and imaginary parts of 𝑦 n+1 , 𝑦 ′ n+1 and 𝑦(x n+1 ), 𝑦 ′ (x n+1 ), respectively, we get the following system (see Tsitouras 11 ):…”