function solve_BVP()
  % Vorgabe für Lösung: 6 Punkte für x von 0 bis 5
  %   konstanter Funktionswert 1, Ableitung 0
  solinit = bvpinit( [0 1 2 3 4 5], [1 0] );

  % Solver-Aufruf
  sol = bvp4c( @Fun_BVP, @BC_BVP, solinit );
  % Stützpunkte zum Zeichnen des x-Bereichs
  x = linspace(0,5,100);
  w = deval( sol, x );  % dazu Lösungswerte aus sol holen
  y = w(1,:);           % y ist 1. Komponente der Lösung
  % alles plotten
  plot( x,y )
end