%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)
%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],

for c=1:1:length(x0)
x_hist=[x_hist,x]

    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