“…Many attempts have been made to formulate numerical algorithms capable of solving special problem of type (1) without reducing it to system of first order equations [7], [9], [11], [12] . Multiderivative methods for direct solution of problem (1), taking m = 4, have been considered using canonical polynomial as the basis functions [4], [5] . They produced methods with many function evaluations.…”