Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function vp = CalculezValProprii(d,s,m, tol)
- n = length(d);
- if(m > length(d))
- m = length(d)
- endif
- vp = zeros(1, m)
- r = IntervaleValProprii(d,s,m)
- for k = 1:m
- Mvechi = 0
- Mnou = (r(k)+r(k+1))/2
- while((abs(Mnou-Mvechi)/Mnou)<= tol)
- Pa = ValoriPolinoame(d,s,r(k));
- Pm = ValoriPolinoame(d,s,Mnou);
- if((Pa(n+1)*Pm(n+1)) < 0)
- r(k+1) = Mnou
- elseif((Pa(n+1)*Pm(n+1)) > 0)
- r(k) = Mnou
- else
- vp(k) = Mnou
- break;
- endif
- Mvechi = Mnou
- Mnou = (r(k)+r(k+1))/2
- endwhile
- endfor
- endfunction
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement