function ex08b()
% Parameter
c  = 20;
m = 1;
L = 0.3;
g = 9.81;
L1 = 0.25;
L2 = 0.2;

% Definition der Funktion F = 2d-Funktion mit 2d-Argument
wurzp = @(x, y) sqrt((x + L).^2 + y.^2);
wurzm = @(x, y) sqrt((x - L).^2 + y.^2);
f1 = @(x,y) (1 - L1./wurzp(x, y)).*(x + L) ...
          + (1 - L2./wurzm(x, y)).*(x - L);
f2 = @(x,y) (1 - L1./wurzp(x, y)).*y ...
          + (1 - L2./wurzm(x, y)).*y + m*g/c; 
F = @(X) [f1(X(1), X(2)); f2(X(1), X(2))];

% Jacobimatrix DF = 2x2-Matrix-Funktion mit 2d-Argument
f1x = @(x,y) 2 - L1*y.^2./(wurzp(x, y).^3) ...
               - L2*y.^2./(wurzm(x, y).^3);
f1y = @(x,y) L1*(x + L).*y./(wurzp(x, y).^3) ...
           + L2*(x - L).*y./(wurzm(x, y).^3);
f2x = @(x,y) f1y(x,y);
f2y = @(x,y) 2 - L1*(x + L).^2./(wurzp(x, y).^3) ...
               - L2*(x - L).^2./(wurzm(x, y).^3);
DF = @(X) [f1x(X(1),X(2)), f1y(X(1),X(2));...
           f2x(X(1),X(2)), f2y(X(1),X(2))];

% Plotten
xvals = -0.7:0.05:0.7;
yvals = -0.7:0.05:0.3;
plotZeros(f1, f2, xvals, yvals)
frame = getframe(gcf);
imwrite(frame.cdata, "bild87.png");

% Startwert ablesen:
x0 = [0; -0.4];

% Lösen mit Newton-Verfahren
x = solveNewton(F, DF, x0, 1e-6)
