%Aufgabe 2b mainfile marco Weber

clear all 
clc

tol=10^(-10);

x0=[0.5;0.4];

maxit=1000;

f=@(x)[exp(x(1)*x(2))+x(1)^2+x(2)-1.4;x(1)^2+x(2)^2+x(1)-0.46];
Df=@(x)[x(2)*exp(x(1)*x(2))+2*x(1),x(1)*exp(x(1)*x(2))+1 ;2*x(1)+1, 2*x(2)];


[x_n,k]=my_newton(Df,f,x0,tol,maxit)

[x_n,k]=my_quasi(Df,f,x0,tol,maxit)

%Lösung Aufgabe 2b.

% x_n =
%    0.310757421718111
%    0.229504690864918
% k =
%      4

%quasi newton
% x_n =
%    0.310757421736084
%    0.229504690889023
% k =
%     16

