VPF::Не могу разобраться с разложением функции в ряд

ErioL

Дата 10.10.2007, 16:41 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 9
Регистрация: 10.10.2007

Репутация: нет
Всего: нет

Найти значение функции y=X*arctgX-ln(sqrt(x*x+1)), разложенной в ряд с общей формулой
((-1)^(n+)*(x^(2*n))/((2*n(2*n-1)))

Код


uses crt; 
{----------------x v stepeni 2n----------------------} 
function Power(Base:real;Exponent:integer):real; 
begin 
if Exponent=0 then Power := 1 {n***0+1} 
else if(Base = 0) and (Exponent > 0) then 
Power := 0 {0**n=0,n>0) } 
else Power := Exp(Exponent*Ln(Base)); 
{if (Exponent div 2)and (Base=-1) then power:=1 
else power:=-1;} 
end; 
{-----------------osnovnaya proga--------------------} 
var c,x,y,s,b:real; 
a:real;i,j,n:integer; 
begin clrscr; 
n:=10; 
x:=0.01; 
writeln(' po formule sleduet'); 
writeln; 
for i:=1 to 9 do begin 
x:=x+0.0885; 
y:=X*arctan(X)-ln(sqrt(X*X+1)); 
writeln('x= ',x:5:5,' s:= ',y:5:5); 
end; 
s:=1; 
x:=0.01; 
{ b:=X*arctan(X)-Ln(Sqrt(X));} 
writeln( 'po razlojeniyu v ryad sleduet'); 
writeln; 


for i := 1 to 10 do begin 
if (i=2*(i div 2))and(i<>2) then i:=-i; {почему то здесь у ходит в бесконечный цикл} 
x:=x+0.0885; 

a:=(Power(x,2*i))/(2*i*(2*i-1)); 
s:=s+a; 
writeln(' x= ',x:5:5,' s= ',s:5:5); 
end; 

readkey; 
end.

ПОМОГИТЕ я ужо фпарился с ней.

Теги не забывай...

Это сообщение отредактировал(а) volvo877 - 10.10.2007, 17:38

volvo877

Дата 10.10.2007, 18:06 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Комодератор
Сообщений: 2073
Регистрация: 15.11.2004

Репутация: 2
Всего: 116

И опять же не нужны здесь никакие дополнительные функции... Достаточно знать, как изменяется общий член ряда при переходе от i-го элемента к (i+1)-му...

Код


uses crt;
const n = 10;

var
  X, next, Y, s: real;
  i, k: integer;
begin
  clrscr;
  writeln(' po formule sleduet');
  writeln;

  x := 0.01;
  for i:=1 to 10 do begin
    x := x + 0.0885;
    y := X*arctan(X)-ln(sqrt(X*X+1));
    writeln('x = ', x:7:5,' s = ', y:7:5);
  end;

  writeln( 'po razlojeniyu v ryad sleduet');
  writeln;

  x := 0.01;
  for i := 1 to n do begin
    x:=x+0.0885;
    next := -1; k := 0;
    s := 0;
    repeat
      inc(k);
      next := ((-1) * next * sqr(x));
      s := s + next / (2 * k * (2*k - 1));
    until abs(next / (2 * k * (2*k - 1))) < 0.0001;
    writeln('x = ', x:7:5, ' s = ', s:7:5);
  end;
  readkey;
end.

Что непонятно - спрашивай...

ErioL

Дата 11.10.2007, 12:33 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 9
Регистрация: 10.10.2007

Репутация: нет
Всего: нет

Гигантское спасибо, я как раз хотел енто седня реализовать, потому что кенты надо мной прикалывались))))

Добавлено через 49 секунд
Только один вопрос, как енто в функции выглядеть будет?
Спасибо заранее

ErioL

Дата 14.10.2007, 12:54 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 9
Регистрация: 10.10.2007

Репутация: нет
Всего: нет

Все работает, сделал сам))
Спасибо большое)))

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Object Pascal: кроссплатформенные технологии \| Следующая тема »