![$z(x)=c_1 y(x)+ c_2 y(x) \int_0^x \frac{ds}{y^2(s)} + y(x)\int_0^x \frac{\int_0^s y(t) f(t) dt \, ds}{y^2(s)}$ $z(x)=c_1 y(x)+ c_2 y(x) \int_0^x \frac{ds}{y^2(s)} + y(x)\int_0^x \frac{\int_0^s y(t) f(t) dt \, ds}{y^2(s)}$](https://dxdy-04.korotkov.co.uk/f/b/c/5/bc54b03f7c6877517954535da0fe261882.png)
Формула не точная (вспомнил по памяти)
Формула не то чтобы неточная, но довольно невнятно выписанная: последнее слагаемое на самом деле должно выглядеть как
![$y(x)\int\limits_0^x \frac{ds}{y^2(s)}\int\limits_0^s y(t) f(t) dt$ $y(x)\int\limits_0^x \frac{ds}{y^2(s)}\int\limits_0^s y(t) f(t) dt$](https://dxdy-03.korotkov.co.uk/f/a/0/1/a01349195c2f095e5aaff4ee8a4fc0e182.png)
. Она получается из стандартной формулы метода вариации произвольных постоянных:
![$z(x)=C_1y_1(x)+C_2y_2(x)-y_1(x)\int\limits_0^x\frac{f(t)y_2(t)}{w}\,dt+y_2(x)\int\limits_0^x\frac{f(t)y_1(t)}{w}\,dt,$ $z(x)=C_1y_1(x)+C_2y_2(x)-y_1(x)\int\limits_0^x\frac{f(t)y_2(t)}{w}\,dt+y_2(x)\int\limits_0^x\frac{f(t)y_1(t)}{w}\,dt,$](https://dxdy-01.korotkov.co.uk/f/0/a/2/0a2d98aafc6d50e11b0d4b6e957a112882.png)
где
![$w=y_1(x)y_2'(x)-y_1'(x)y_2(x)$ $w=y_1(x)y_2'(x)-y_1'(x)y_2(x)$](https://dxdy-03.korotkov.co.uk/f/6/8/a/68a7f46afc9524ae147762050e56b51282.png)
-- вронскиан двух независимых решений однородного уравнения. И если теперь выразить второе решение через первое из условия, например,
![$w\equiv1$ $w\equiv1$](https://dxdy-02.korotkov.co.uk/f/d/d/a/dda0b222ec5805d3a5dc8d690a625dda82.png)
, т.е. взять
![$y_2(x)=y_1(x)\int\limits_0^x\frac{ds}{y_1^2(s)}$ $y_2(x)=y_1(x)\int\limits_0^x\frac{ds}{y_1^2(s)}$](https://dxdy-04.korotkov.co.uk/f/7/6/c/76c8e929a16591ac8aaf77caf46cef6e82.png)
, то после подстановки там появятся два двойных интеграла, причём в третьем слагаемом интеграл будет браться по треугольнику, а в четвёртом -- по квадрату. Интегралы по общему треугольнику сократятся, и останется то, что нужно.
Это если по существу; а ссылку на литературу дать, к сожалению, не могу.