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;
- a = r(k);
- b = r(k+1);
- while(1)
- Pa = ValoriPolinoame(d,s,a);
- Pm = ValoriPolinoame(d,s,Mnou);
- if((Pa(n+1)*Pm(n+1)) < 0)
- b = Mnou;
- elseif((Pa(n+1)*Pm(n+1)) > 0)
- a = Mnou;
- endif
- if((abs(Mnou-Mvechi)/Mnou)<= tol)
- vp(k) = Mnou;
- break;
- endif
- Mvechi = Mnou;
- Mnou = (a+b)/2;
- endwhile
- endfor
- endfunction
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement