clear;
F=zeros(2,1);
e=zeros(2,1);
jak=zeros(2,2);
 a=[];
 %%%Maksimum İterasyon Sayısı%%%
iterasyon=50;
%%%ilk degerler%%%
e(1)=0.3;%nn
e(2)=1.5;%cc
%%%%%parametreler%%%
alpha1=3;
alpha2=2;
%%%Cizdirmek icin matris%%%
a(1,2)=e(1);
a(1,3)=e(2);
%%%Hata Tolerans%%%
epsilon=0.001;

for j=1:iterasyon
    %%%denklemler%%%
    F=[alpha1*(e(2)/1+e(2))*e(1)-e(1); 
    alpha2-(e(2)/1+e(2))*e(1)-e(2)];
    %%%denklemlerin jakobiyeni%%%
    jak=[alpha1*((e(2)/1+e(2))-1)  alpha1*e(1)/((1+e(2))^2); (-e(2)/(1+e(2)))  -1-(1/((1+e(2))^2))];
    %%%Newton Raphson%%% 
    e=e-(inv(jak)*F);
    %%%Epsilon ile norm karsilastirma%%%
    xnorm=norm(-inv(jak)*F);
    a(j+1,2)=e(1);
    a(j+1,3)=e(2);
    if xnorm<=epsilon
        break
    else j=j+1;
    end
    
end
%%%İterasyon sayisi=j%%%
j 
%%%Ulasilan Çözüm Değerlerini Goster%%%
xnorm
e
%%%Matlab Grafik Çizdirme%%%
plot(a)