function solve_BVP1a()
    % Vorgabe: 10 Punkte für x von 0 bis xInt
    xInt = 5;
    xinit = linspace( 0, xInt, 10 );
    % Übergabe der Funktion initBVP1a
    solinit = bvpinit( xinit, @initBVP1a );

    % Solver-Aufruf
    sol = bvp4c( @Fun_BVP1, @BC_BVP1, solinit );

    % Stützpunkte zum Zeichnen des x-Bereichs
    x = linspace(0,xInt,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, [0,xInt], [0,0], 'k:' )
end
	