Turbo Pascal

sepek

Program polega na tym ze ma rozwiazac uklad rownan nie liniowych metoda Newtona

F(x,y)=2x(do potegi 2)y+3xy-2xy(do potegi 2)-3=0
G(x,y)=x+3y(do potegi 3)x(do potegi 2)-2yx(do potegi 2)-1=0

dla x0=1.1 i y0=1.05

czy ktos zna dobrze ta metode. pomozcie
ja napisalem go ale cos jest w nim nie tak bo nie kalkuluje. sprawdzcie go
jak mozecie. pozdro

{$N+}
program newton;
uses crt; {wczytujemy potrzebne biblioteki}
var
x, x1, y, y1,s,E:extended; {deklarujemy zmienne}
wynik :text;
z: string;
label 1,2,st,kon; {deklarujemy kotwice}
begin {rozpoczynamy program}
st:
clrscr; {czyscimy ekran}
writeln('Warunki poczatkowe: ');
writeln('x0=1.1; y0:=1.05');

x1:=1.1;
y1:=1.05;

write('Podaj dokladnosc z jaka okreslic wynik ');
readln(E); {wczytujemy zmienna}
assign(wynik,'wynik.txt'); {przypisujemy sciezke do pliku}
rewrite(wynik); {otwieramy i kasujemy zawartosc pliku}
x:=0;
s:=0;
repeat {rozpoczynamy petle repeat}
if abs(x1-x)<=E then goto 1;
x:=x1;
y:=y1;

x1:=4*y*x+3*y-2*y*y; {pochodna pierwszego rownania po x}

write(wynik,x1:1:5,' '); {wpisujemy dane do pliku}
1:
if abs(y1-s)<=E then goto 2;
y:=y1;
s:=y;
y1:=9*y*y*x*x-2*x*x; {pochodna drugiego rownania po y}
write(wynik,y1:1:5); {wpisujemy dane do pliku}
2:
writeln(wynik)
until (abs(x1-x)<=E) and (abs(y1-y)<=E);
writeln('x= ',x1:1:5);
writeln('y= ',y1:1:5);
close(wynik); write('rozwiazac ponownie? T/N ');
readln(z);
if z='t'then goto st;
if z='n'then goto kon;
kon:
end.

sepek

blad w nim polega na tym ze poprostu krzyczy ze nie moze kalkulowac bo zaduza liczba, nie wiem moze nie takim kodem i nie takich jak trzeba komend uzywam. prosze was o pomoc musze oddac ten program na zaliczenie przedmiotu a nie wiem jak go ugryzc. Pozdrawiam wszystkich

				fora.pl - załóż własne forum dyskusyjne za darmo Powered by phpBB © 2001, 2002 phpBB Group Avalanche style by What Is Real © 2004