“…Q pqr (t) y (p) (t) y (q) (t) y (r) is considered, where P k (t) , Q pq (t) , Q pqr (t) and g (t) are the given analytic functions defined on the interval a ≤ t ≤ b; λ j , a kj and b kj are the known rael coefficients. In order to solve the nonlinear problem (1.1)-(1.2), we utilize the matrix-collocation method, which have been developed by Sezer and Coworkers [7,10,[12][13][14], and research the numerical solution in the truncated Morgan-Voyce series form…”