Program rozwiązujący funkcje matematyczną w pascalu

[jak781 0]

Witam. Napisałem taki oto program rozwiązujący funkcje ze zdjęcia. Jednak nie jestem przekonany, czy jest on do końca poprawny. Mianownik zmienia się o 2! i nie wiem jak mam to zapisać w pascalu.

var x,eps,suma,wyraz,i:real;
begin
   writeln('Podaj x');
   readln(x);
   i:=1;
   eps:=0.001;
   wyraz:=1;
   suma:=wyraz;
   while wyraz>eps do
   begin
       wyraz:=wyraz*(-1/i)*x*x;
       suma:=suma+wyraz;
       i:=i*(i*1);
   end;
   writeln('Wartosc funkcji wynosi',suma);
end.

 

 

Bardzo proszę o pomoc! Pozdrawiam.

[_Piotr_Z_ 0]

Od dobrych 15 lat nie bawiłem się takimi rzeczami, więc cos mogło mi się pomieszać, ale:

1. element zerowy powinien być raczej równy x, a nie 1
   
2. (-1/i) - a to skąd Ci się wzięło?
   
3. dodaj kolejną zmienna, w której będziesz przechowywał kolejne elementy silni (a(i) = a(i-1)*(4*i*i + 2*i) - chyba, policz sam)
   

 

Tak w ogóle, to wydaje mi się, że Twoim głównym problemem jest matematyka, a nie jej zapis w Pascalu. Wylicz sobie jak dokładnie zmieniają się kolejne wyrazy ciągu, a wtedy zaimplementowanie tego stanie się proste.

[Bono[UG] 1920]

Proponuję zmienne w programie nazwać tak samo jak we wzorze, będzie łatwiej połapać się co i gdzie jest liczone.

Do liczenia silni i potęgi dobrze mieć zdefiniowane funkcje, raz czytelniej, dwa bez dziwnych kombinacji (na razie niech program działa prawidłowo, potem można kombinować nad wydajnością). Nie pamiętam już ale pewnie są już gotowce w samym pascalu lub dodatkowych bibliotekach (uses math ).

 

Kolejna sprawa, to kiedy przestać liczyć, czyli co porównać z założonym epsilonem.

Nie wiem na jakim poziomie nauki jesteś, liceum czy studia i czy miałeś już coś z metod numerycznych. Ogólna zasada jest taka, że trzeba sprawdzić zbieżność liczonego szeregu, czyli po prostu zbadać różnicę między poprzednią wartością szeregu, a obecnie wyliczoną. Jeżeli będzie <epsilon, to kończy się obliczenia i zwraca wynik.

[Elf_Gektor 0]

i := 1;

 

x - nieznana stała

 

i:=i*(i*1);

i:=1*(1*1);

i = 1; ???

 

Jaki to ma sens?