function [F, DF] =Fun(x)
% Kombi-Routine fuer das Newtonverfahren,
% berechnet F(x) und die Jacobimatrix DF(x)
% Zwei Ellipsen
F(1,1)=   x(1)^2/1+x(2)^2/4-1;
F(2,1)=   x(1)^2/2+x(2)^2/6-1;
DF(1,1)=2*x(1)/1;
DF(1,2)=2*x(2)/4;
DF(2,1)=2*x(1)/2;
DF(2,2)=2*x(2)/6;