%aufgabe 3 Serie 5 Marco Weber
clear all
clc


F=@(x)[(atan(2*x(1))-(x(2)/15)^3) ; 1/30*(x(1)-x(2))];
Df=@(x) [2/(4*x(1)^2+1) -x(2)^2/1125 ;1/30 -1/30]; 
%a)

figure(1)
x1=[-15:0.25:15];
y1=[-15:0.25:15];

z1=[atan(2*x1)-(y1/15).^3];

plot3(x1,y1,z1)
grid on;
hold on;


figure(2)
x2=[-15:0.25:15];
y2=[-15:0.25:15];

z2=[1/30.*(x2-y2)];

plot3(x2,y2,z2)
grid on;
hold on;

%___________________________________________________________

maxit=5; %Nur die ersten 5 Iterationen
tol=10^(-10); %gleich wie bei aufgabe 2 

x0=[[0.5 0.7 ; 5 10]];
alpha=[0.5 0.8 1];

% beta=length(x0) Orientierung ob i.O.

for c=1:1:length(x0)


    for m=1:1:length(alpha)
        
    [x_n,x_hist,k]=my_damp_newton(F,Df,x0(c),tol,maxit,alpha(m));
    
    hist =[hist, x_hist'];
    end
hist
end