www.imir12.fora.pl

Psulek

Gdyby ktoś potrzebował to tak wygląda program z ostatnich ćwiczeń który działał ;]

p=[6 3 0 2] % funkcja
x1=-2 % wstępna wartość x
b=1 % wstępna wartość błędu
k=0 % wstępna ilość przebiegu pętli
while b>0.0000000000001 % warunek dokładności błędu
k=k+1; % ilość przebiegu pętli +1
f=polyval(p,x1); % wart. funkcji w p-kcie
f1=polyder(p); % pochodna
f2=polyval(f1,x1); % wart. poch. w p-kcie
h=f/f2; % f/f'
x2=x1-h; % wzór newtona
o=x1-x2; % różnica wartości
b=abs(o); % wart. bezwzględna różnicy
x1=x2; % przepisanie wartości do następnego przebiegu pętli
end % koniec pętli
k % wypisanie ilości przebiegów pętli
x1 % wypisanie wartości obliczonej ze wzoru
roots(p) % wypisanie wartości obliczonej inną metodą przez program ( do porównania dokładności )

Miłego analizowania ;] Pozdrawiam

Skrzydlaty

Dzieki!