Advertisement
Guest User

Untitled

a guest
Apr 16th, 2017
79
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Octave 0.66 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.    
  20.     while((abs(Mnou-Mvechi)/Mnou)<= tol)
  21.    
  22.     Pa = ValoriPolinoame(d,s,r(k));
  23.    
  24.     Pm = ValoriPolinoame(d,s,Mnou);
  25.    
  26.       if((Pa(n+1)*Pm(n+1)) < 0)
  27.        r(k+1) = Mnou
  28.       elseif((Pa(n+1)*Pm(n+1)) > 0)
  29.         r(k) = Mnou
  30.        else
  31.        vp(k) = Mnou
  32.           break;
  33.         endif
  34.          Mvechi = Mnou
  35.          Mnou = (r(k)+r(k+1))/2
  36.       endwhile
  37.    
  38.   endfor
  39.  
  40. endfunction
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement