Solving the Telegraph and Oscillatory Differential Equations by a Block Hybrid Trigonometrically Fitted Algorithm

We propose a block hybrid trigonometrically fitted (BHT) method, whose coefficients are functions of the frequency and the step-size for directly solving general second-order initial value problems (IVPs), including systems arising from the semidiscretization of hyperbolic Partial Differential Equations (PDEs), such as the Telegraph equation. The BHT is formulated from eight discrete hybrid formulas which are provided by a continuous two-step hybrid trigonometrically fitted method with two off-grid points. The BHT is implemented in a block-by-block fashion; in this way, the method does not suffer from the disadvantages of requiring starting values and predictors which are inherent in predictor-corrector methods. The stability property of the BHT is discussed and the performance of the method is demonstrated on some numerical examples to show accuracy and efficiency advantages.