Advertisement
Guest User

Untitled

a guest
Apr 16th, 2017
71
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Octave 0.54 KB | None | 0 0
  1. function vp = CalculezValProprii(d,s,m, tol)
  2.  
  3. if(m > length(d))
  4.   m = length(d)
  5. endif
  6. vp = zeros(1, m)
  7.   r = IntervaleValProprii(d,s,m)
  8.  
  9.   for k = 1:m
  10.    
  11.     Mvechi = 0
  12.     Mnou = (r(k)+r(k+1))/2
  13.    
  14.     %while((norm(Mnou-Mvechi)/Mnou)<= tol)
  15.    
  16.    %   if((r(k)*Mnou) < 0)
  17.     %    r(k+1) = Mnou
  18.     %  elseif((r(k)*Mnou) > 0)
  19.     %    r(k) = Mnou
  20.     %   else
  21.     %   vp(k) = Mnou
  22.     %      break;
  23.     %     endif
  24.     %     Mvechi = Mnou
  25.     %     Mnou = (r(k)+r(k+1))/2
  26.     %  endwhile
  27.    
  28.   endfor
  29.  
  30. endfunction
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement