clear

function s = f(x)
    s(1,1) = -5.*x(1) + 2.*sin(x(1)) + 2.*cos(x(2));
    s(2,1) = 2.*cos(x(1)) + 2.*sin(x(2)) -5.*x(2); 
endfunction

function A = Df(x)
  //  A(1,1) = 
  //  A(1,2) =  
  //  A(2,1) = 
  //  A(2,2) = 
endfunction


function s = g(x)
  //  s(1,1) = 
  //  s(2,1) = 
endfunction


// itéré initial 

xinit = zeros(2,1);
//xinit = [10000;-10000];


///////////////////////////////////////////
/
// Resolution par point fixe contractant 
//
///////////////////////////////////////////


eps = 1.0E-14; // critere d'arret 
kmax = 100; // nb max d'iterations 

x = xinit 

// A COMPLETER 


while (nr/nr0 > eps)&(k<kmax)
    
   
   // A COMPLETER 
      
    
end

scf(1)
xtitle('residu')
//plot(res1,'-xb')  // a decommenter 
//print(%io(2),x,r)



////////////////////////////////
//
//  NEWTON 
//
///////////////////////////////




// Resolution par Newton 



eps = 1.0E-14;
kmax = 100;


x = xinit 
// A COMPLETER 

while (nr/nr0 > eps)&(k<kmax)



   // A COMPLETER 

     
   
 
end
scf(1)

//plot(res,'-xr')
//print(%io(2),x,r)



