function dy_dx = Fun_BVP( x, y )
  % Ableitung von z = y(1) ergibt v = y(2)
  dy_dx(1,1) = y(2);
  % Ableitung von v = y(2) über die DGL = Krümmung -2
  dy_dx(2,1) = -2;
end