Wielomian 3 stopnia – program wiesza się

Wielomian 3 stopnia – program wiesza się
0

Witam,
znajomy dostał od prowadzącego na uczelni zadanie do wykonania (załącznik), a ja staram się mu je rozwiązać. Polega ono na znalezieniu koncentracji składnika c dla podanego równania.
Metodą jaką wybrałem jest metoda Newtona. Napisałem kod w Pascalu niestety nie do końca działa on tak jak bym chciał. Przy wprowadzeniu wartości początkowej c dla której algorytm ma zacząć obliczać miejsce zerowe i dokładności program wiesza się jeśli dokładność jest większa niż 0.001. Nie wiem czy to wina mojego słabego laptopa czy to jakiś błąd w kodzie. Bardzo proszę o pomoc.

Kopiuj
program Untitled;

uses crt;

type
	tablica = array[1..100000] of real;

var
	c: tablica;
	Q,cin,k,K1,V,f,g,z,eps: real;
	n: integer;

begin
	clrscr;

WriteLn('Program oblicza koncentracje koncowa skladnika c w reaktorze z idealnym wymieszaniem metoda Newtona.');
Write('Podaj Q: ');
ReadLn(Q);
Write('Podaj cin: ');
ReadLn(cin);
Write('Podaj k: ');
ReadLn(k);
Write('Podaj K1: ');
ReadLn(K1);
Write('Podaj V: ');
ReadLn(V);
Write('Podaj x: ');
ReadLn(c[1]);
Write('Podaj dokladnosc obliczen: ');
ReadLn(eps);

n:=1;

repeat
	f:=Q*c[n]+2*Q*K1*c[n]*c[n]+Q*K1*K1*c[n]*c[n]*c[n]-Q*cin-2*Q*cin*K1*c[n]-Q*cin*K1*K1*c[n]*c[n]+V*k*c[n];
      	g:=Q+4*Q*K1*c[n]+3*Q*K1*K1*c[n]*c[n]-2*Q*cin*K1-2*Q*cin*K1*K1*c[n]+V*k;
      	c[n+1]:=c[n]+(f/g);
      	z:=c[n+1]-c[n];
until
      	z<eps;

WriteLn('Szukane stezenie c wynosi: ',c[n+1]);

ReadLn();

end.

edytowany 1x, ostatnio: flowCRANE
flowCRANE
Moderator Delphi/Pascal
  • Rejestracja:ponad 13 lat
  • Ostatnio:około 2 godziny
  • Lokalizacja:Tuchów
  • Postów:12166
0

Wartość zmiennej n nie zmienia się w tej pętli, a powinna.

Poza tym, ten kod jest dość nieczytleny i pasowałoby go sformatować.


Pracuję nad własną, arcade'ową, docelowo komercyjną grą z gatunku action/adventure w stylu retro (pixel art), programując silnik i powłokę gry od zupełnych podstaw, przy użyciu Free Pascala i SDL3. Więcej informacji znajdziesz na moim mikroblogu.
edytowany 1x, ostatnio: flowCRANE
0

Dodałem na końcu pętli inkrementację zmiennej n jednak wiele to nie pomogło. Teraz po uruchomieniu programu i wczytaniu zmiennych program wyrzuca błąd po którym konsola automatycznie się zamyka przez co nie jestem nawet w stanie zorientować się w czym rzecz. Co do formatowania kodu, to nie jestem programistą. Starałem się zastosować do instrukcji publikacji kodu, jeśli czegoś nie zrobiłem to powiedz proszę w czym rzecz, a następnym razem na pewno zrobię to lepiej.

Kopiuj
repeat
      f:=Q*c[n]+2*Q*K1*c[n]*c[n]+Q*K1*K1*c[n]*c[n]*c[n]-Q*cin-2*Q*cin*K1*c[n]-Q*cin*K1*K1*c[n]*c[n]+V*k*c[n];
      g:=Q+4*Q*K1*c[n]+3*Q*K1*K1*c[n]*c[n]-2*Q*cin*K1-2*Q*cin*K1*K1*c[n]+V*k;
      c[n+1]:=c[n]+(f/g);
      z:=c[n+1]-c[n];
      n:=n+1;
until
      z<eps;
flowCRANE
Moderator Delphi/Pascal
  • Rejestracja:ponad 13 lat
  • Ostatnio:około 2 godziny
  • Lokalizacja:Tuchów
  • Postów:12166
1
Atsutane napisał(a):

Teraz po uruchomieniu programu i wczytaniu zmiennych program wyrzuca błąd po którym konsola automatycznie się zamyka przez co nie jestem nawet w stanie zorientować się w czym rzecz.

Jak to nie – uruchom program z wiersza poleceń i zobacz co wyrzuciło. Albo skorzystaj z debuggera.

Co do formatowania kodu, to nie jestem programistą.

Nie trzeba być programistą, aby widzieć, że jeden ciąg znaków długaśnego wyrażenia bardzo źle się czyta. Przydałoby się chociaż rozsunąć fragmenty pomiędzy odejmowaniem – już łatwiej będzie.


Sprawdź ten kod pod debuggerem – postaw breakpoint na pierwszej linijce w pętli, wrzuć sobie zmienne do okna watches i klawiszem F8 wykonuj po jednej linijce, obserwując co ostatecznie wpisywane jest do zmiennych.

Poza tym, warunek końcowy pętli nie zabezpiecza przed wyjściem poza zakres tablicy, przez co n może być większe niż 100.000. Dlatego też na sam początek zmień ten warunek na taki:

Kopiuj
until (n = 100000) or (z < eps);

Końcowe wyświetlanie danych też trzeba zabezpieczyć, bo również możesz wyjść poza zakres tablicy tym c[n+1].


Pracuję nad własną, arcade'ową, docelowo komercyjną grą z gatunku action/adventure w stylu retro (pixel art), programując silnik i powłokę gry od zupełnych podstaw, przy użyciu Free Pascala i SDL3. Więcej informacji znajdziesz na moim mikroblogu.
edytowany 5x, ostatnio: flowCRANE
CL
  • Rejestracja:ponad 15 lat
  • Ostatnio:7 miesięcy
0

c[n+1]:=c[n]+(f/g);
A co to za wersja metody Newtona??
chyba miało być
c[n+1]:=c[n]-(f/g);

To samo z warunkiem zakończenia. Różnica z:=c[n+1]-c[n]; może wyjść -500 i wyjdzie z pętli a warunek nie będzie spełniony dla metody. Wartość różnicy ma być bezwzględna!

edytowany 1x, ostatnio: Clarc
GS
  • Rejestracja:ponad 14 lat
  • Ostatnio:około 3 godziny
1

W jakim celu zapisujesz do tablicy o stałym rozmiarze kolejne wyniki? Raz że ryzykujesz wyjście poza jej wymiar, dwa, te dane nie są do niczego potrzebne. Wystarczy że będziesz pamiętał w zmiennej wynik z poprzedniej iteracji.

0

W metodzie Newtona rzeczywiście miałem zły znak, wyłapałem to jak przeliczałem sobie algorytm w głowie. Co do wartości absolutnej zupełnie o tym zapomniałem wobec czego stokrotne dzięki za wskazówkę.
Co się zaś tyczy tablicy o stałym rozmiarze to nie wiem czemu ją stosuję. Zawsze jej używałem, ponieważ tak było mi najprościej. Obliczenia które wykonywał program zawsze sprawdzałem sobie w excel'u czy program wyjdzie poza skalę i odpowiednio dobierałem rozmiar tablicy. Nie bardzo wiedziałem jednak o co chodziło Panu powyżej więc zrobiłem to na czuja i program działa. Jednakże, jeśli nie o to chodziło to chciałbym prosić o informację aby w przyszłości nie popełnić podobnego błędu. Na sam koniec mam pytanie dotyczące tablicy o stałym rozmiarze. Znalazłem już informację, że tablice o zmiennym rozmiarze istnieją, jednak czy używanie ich jest na poziomie laika jakim jestem, jeśli chodzi o programowanie?

Kopiuj
program Untitled;

uses crt;

var

Q,cin,k,K1,V,f,g,z,eps,c,m: real;
l: integer;
p: char;

begin

repeat
	clrscr;
	WriteLn('Program oblicza koncentracje koncowa skladnika c metoda Newtona.');
	Write('Podaj Q: ');
	ReadLn(Q);
	Write('Podaj cin: ');
	ReadLn(cin);
	Write('Podaj k: ');
	ReadLn(k);
	Write('Podaj K1: ');
	ReadLn(K1);
	Write('Podaj V: ');
	ReadLn(V);
	Write('Podaj x: ');
	ReadLn(m);
	Write('Podaj dokladnosc obliczen: ');
	ReadLn(eps);
	
	repeat
	      f:=Q*m + 2*Q*K1*m*m + Q*K1*K1*m*m*m - Q*cin - 2*Q*cin*K1*m - Q*cin*K1*K1*m*m + V*k*m;
	      g:=Q + 4*Q*K1*m + 3*Q*K1*K1*m*m - 2*Q*cin*K1 - 2*Q*cin*K1*K1*m + V*k;
	      c:=m - (f/g);
	      z:=abs(c - m);
	      m:=c;
	until
	      z<eps;
	
	WriteLn('Ile miejsc po przecinku program ma wyswietlic w wyniku?');
	ReadLn(l);
	clrscr;
	WriteLn('Szukane stezenie c wynosi: ',c:l+2:l);
	WriteLn();
	WriteLn('Czy chcesz zmienic warunki zadania? T/N');
	ReadLn(p);
until
	p='N';

ReadLn();


end.

KU
  • Rejestracja:około 7 lat
  • Ostatnio:ponad 6 lat
  • Postów:369
0
Kliknij, aby dodać treść...

Pomoc 1.18.8

Typografia

Edytor obsługuje składnie Markdown, w której pojedynczy akcent *kursywa* oraz _kursywa_ to pochylenie. Z kolei podwójny akcent **pogrubienie** oraz __pogrubienie__ to pogrubienie. Dodanie znaczników ~~strike~~ to przekreślenie.

Możesz dodać formatowanie komendami , , oraz .

Ponieważ dekoracja podkreślenia jest przeznaczona na linki, markdown nie zawiera specjalnej składni dla podkreślenia. Dlatego by dodać podkreślenie, użyj <u>underline</u>.

Komendy formatujące reagują na skróty klawiszowe: Ctrl+B, Ctrl+I, Ctrl+U oraz Ctrl+S.

Linki

By dodać link w edytorze użyj komendy lub użyj składni [title](link). URL umieszczony w linku lub nawet URL umieszczony bezpośrednio w tekście będzie aktywny i klikalny.

Jeżeli chcesz, możesz samodzielnie dodać link: <a href="link">title</a>.

Wewnętrzne odnośniki

Możesz umieścić odnośnik do wewnętrznej podstrony, używając następującej składni: [[Delphi/Kompendium]] lub [[Delphi/Kompendium|kliknij, aby przejść do kompendium]]. Odnośniki mogą prowadzić do Forum 4programmers.net lub np. do Kompendium.

Wspomnienia użytkowników

By wspomnieć użytkownika forum, wpisz w formularzu znak @. Zobaczysz okienko samouzupełniające nazwy użytkowników. Samouzupełnienie dobierze odpowiedni format wspomnienia, zależnie od tego czy w nazwie użytkownika znajduje się spacja.

Znaczniki HTML

Dozwolone jest używanie niektórych znaczników HTML: <a>, <b>, <i>, <kbd>, <del>, <strong>, <dfn>, <pre>, <blockquote>, <hr/>, <sub>, <sup> oraz <img/>.

Skróty klawiszowe

Dodaj kombinację klawiszy komendą notacji klawiszy lub skrótem klawiszowym Alt+K.

Reprezentuj kombinacje klawiszowe używając taga <kbd>. Oddziel od siebie klawisze znakiem plus, np <kbd>Alt+Tab</kbd>.

Indeks górny oraz dolny

Przykład: wpisując H<sub>2</sub>O i m<sup>2</sup> otrzymasz: H2O i m2.

Składnia Tex

By precyzyjnie wyrazić działanie matematyczne, użyj składni Tex.

<tex>arcctg(x) = argtan(\frac{1}{x}) = arcsin(\frac{1}{\sqrt{1+x^2}})</tex>

Kod źródłowy

Krótkie fragmenty kodu

Wszelkie jednolinijkowe instrukcje języka programowania powinny być zawarte pomiędzy obróconymi apostrofami: `kod instrukcji` lub ``console.log(`string`);``.

Kod wielolinijkowy

Dodaj fragment kodu komendą . Fragmenty kodu zajmujące całą lub więcej linijek powinny być umieszczone w wielolinijkowym fragmencie kodu. Znaczniki ``` lub ~~~ umożliwiają kolorowanie różnych języków programowania. Możemy nadać nazwę języka programowania używając auto-uzupełnienia, kod został pokolorowany używając konkretnych ustawień kolorowania składni:

```javascript
document.write('Hello World');
```

Możesz zaznaczyć również już wklejony kod w edytorze, i użyć komendy  by zamienić go w kod. Użyj kombinacji Ctrl+`, by dodać fragment kodu bez oznaczników języka.

Tabelki

Dodaj przykładową tabelkę używając komendy . Przykładowa tabelka składa się z dwóch kolumn, nagłówka i jednego wiersza.

Wygeneruj tabelkę na podstawie szablonu. Oddziel komórki separatorem ; lub |, a następnie zaznacz szablonu.

nazwisko;dziedzina;odkrycie
Pitagoras;mathematics;Pythagorean Theorem
Albert Einstein;physics;General Relativity
Marie Curie, Pierre Curie;chemistry;Radium, Polonium

Użyj komendy by zamienić zaznaczony szablon na tabelkę Markdown.

Lista uporządkowana i nieuporządkowana

Możliwe jest tworzenie listy numerowanych oraz wypunktowanych. Wystarczy, że pierwszym znakiem linii będzie * lub - dla listy nieuporządkowanej oraz 1. dla listy uporządkowanej.

Użyj komendy by dodać listę uporządkowaną.

1. Lista numerowana
2. Lista numerowana

Użyj komendy by dodać listę nieuporządkowaną.

* Lista wypunktowana
* Lista wypunktowana
** Lista wypunktowana (drugi poziom)

Składnia Markdown

Edytor obsługuje składnię Markdown, która składa się ze znaków specjalnych. Dostępne komendy, jak formatowanie , dodanie tabelki lub fragmentu kodu są w pewnym sensie świadome otaczającej jej składni, i postarają się unikać uszkodzenia jej.

Dla przykładu, używając tylko dostępnych komend, nie możemy dodać formatowania pogrubienia do kodu wielolinijkowego, albo dodać listy do tabelki - mogłoby to doprowadzić do uszkodzenia składni.

W pewnych odosobnionych przypadkach brak nowej linii przed elementami markdown również mógłby uszkodzić składnie, dlatego edytor dodaje brakujące nowe linie. Dla przykładu, dodanie formatowania pochylenia zaraz po tabelce, mogłoby zostać błędne zinterpretowane, więc edytor doda oddzielającą nową linię pomiędzy tabelką, a pochyleniem.

Skróty klawiszowe

Skróty formatujące, kiedy w edytorze znajduje się pojedynczy kursor, wstawiają sformatowany tekst przykładowy. Jeśli w edytorze znajduje się zaznaczenie (słowo, linijka, paragraf), wtedy zaznaczenie zostaje sformatowane.

  • Ctrl+B - dodaj pogrubienie lub pogrub zaznaczenie
  • Ctrl+I - dodaj pochylenie lub pochyl zaznaczenie
  • Ctrl+U - dodaj podkreślenie lub podkreśl zaznaczenie
  • Ctrl+S - dodaj przekreślenie lub przekreśl zaznaczenie

Notacja Klawiszy

  • Alt+K - dodaj notację klawiszy

Fragment kodu bez oznacznika

  • Alt+C - dodaj pusty fragment kodu

Skróty operujące na kodzie i linijkach:

  • Alt+L - zaznaczenie całej linii
  • Alt+, Alt+ - przeniesienie linijki w której znajduje się kursor w górę/dół.
  • Tab/⌘+] - dodaj wcięcie (wcięcie w prawo)
  • Shit+Tab/⌘+[ - usunięcie wcięcia (wycięcie w lewo)

Dodawanie postów:

  • Ctrl+Enter - dodaj post
  • ⌘+Enter - dodaj post (MacOS)