Advertisement
Guest User

Untitled

a guest
Apr 16th, 2017
88
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Octave 0.71 KB | None | 0 0
  1. function vp = CalculezValProprii(d,s,m, tol)
  2.  
  3. n = length(d);
  4.  
  5. if(m > length(d))
  6.  
  7.   m = length(d);
  8.  
  9. endif
  10.  
  11. vp = zeros(1, m);
  12.  
  13.   r = IntervaleValProprii(d,s,m);
  14.  
  15.   for k = 1:m
  16.    
  17.     Mvechi = 0;
  18.     Mnou = (r(k)+r(k+1))/2;
  19.     a = r(k);
  20.     b = r(k+1);
  21.    
  22.     while(1)
  23.    
  24.    
  25.     Pa = ValoriPolinoame(d,s,a);
  26.    
  27.     Pm = ValoriPolinoame(d,s,Mnou);
  28.    
  29.       if((Pa(n+1)*Pm(n+1)) < 0)
  30.        b = Mnou;
  31.       elseif((Pa(n+1)*Pm(n+1)) > 0)
  32.         a = Mnou;
  33.        
  34.         endif
  35.         if((abs(Mnou-Mvechi)/Mnou)<= tol)
  36.         vp(k) = Mnou;
  37.         break;
  38.         endif
  39.          Mvechi = Mnou;
  40.          Mnou = (a+b)/2;
  41.       endwhile
  42.    
  43.   endfor
  44.  
  45. endfunction
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement