function solve_BVP1b()
    
    % Vorgabe: 10 Punkte für x von 0 bis xInt
    xInt = 5;
    xinit = linspace( 0, xInt, 10 );
    solinit = bvpinit( xinit, @initBVP1b );
    
    % Solver-Aufruf
    sol = bvp4c( @Fun_BVP1, @BC_BVP1, solinit );
    
    % Gitter für x-Bereich
    x = linspace(0,xInt,100);
    % dazu Lösungswerte holen
    w = deval( sol,x );
    % y ist 1. Komponente der Lösung
    y = w(1,:);
    % alles plotten
    plot( x,y, [0,xInt], [0,0], 'k:' )
end
    
function yinit = initBVP1b( x )
  yinit = [-1 0];
end