// --------------------------------------------------------------------------
// シューティング法を適用する境界値問題の例題の微分方程式を定義する関数
// --------------------------------------------------------------------------
function dz = f(x,z)
dz = [ z(2)                                                 // z(1)'= u'= z(2)
       2 * x./(x^2 + 1).*z(2) - 2/(x^2 + 1).*z(1) + x^2 + 1 // z(2)'= p(x)u'+q(x)u+r(x)
       z(4)                                                 // z(3)'= v'= z(4)
       2 * x./(x^2 + 1).*z(4) - 2/(x^2 + 1).*z(3) ];        // z(4)'= p(x)v'+q(x)v
endfunction