VPF::[Язык любой] Метод численного дифференцирования

eganeev

Дата 3.2.2007, 15:09 (ссылка)

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

Загрузка ...

Новичок

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

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

Кто может напишите плиз прогу по численному дифференцированию! У самого ничего неполучается((( Можете писать на любом языке. Заранее огромное спасибо.

Rodman

Дата 3.2.2007, 15:50 (ссылка)

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

Загрузка ...

CIO

Профиль
Группа: Участник
Сообщений: 6144
Регистрация: 7.5.2006
Где: Ukraine ⇛ Kyiv ci ty

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

Код


program diff;
{$APPTYPE CONSOLE}
const eps=1.0E-6; //заданная точность
function fx(x: real): real; //заданная функция
begin
fx:=sin(x);
end;
function f1(x0: real): real; //вычисление первой производной
var dx,dy,dy2: real;
begin
dx:=1;
repeat
dx:=dx/2;
dy:=fx(x0+dx/2)-fx(x0-dx/2);
dy2:=fx(5*x0/4+dx)-2*fx(5*x0/4);
dy2:=dy2+fx(5*x0/4-dx);
until abs(dy2/(2*dx))<eps;
f1:=dy/dx
end;
Function f2(x0: real): real; // вычисление 2-й производной
var dx,dy,dy3: real;
begin
dx:=1;
repeat
dx:=dx/2;
dy:=fx(x0 + dx)-2*fx(x0)+fx(x0-dx);
dy3:=fx(5*x0/4+2*dx)-2*fx(5*x0/4+dx);
dy3:=dy3-fx(5*x0/4-2*dx)+2*fx(5*x0/4-dx);
until abs(dy3/(6*dx))<eps;
f2:=dy/(dx*dx)
end;
var x: real;
begin
x:=Pi/6;
// пример выполнения
writeln('f1 =',f1(x));
writeln('f` =',cos(x));
writeln('f2 =',f2(x));
writeln('f``=',-sin(x));
readln;
end.

--------------------

CRM с функцией АвтоОбзвона

Лотереи мира

eganeev

Дата 3.2.2007, 18:08 (ссылка)

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

Загрузка ...

Новичок

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

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

Спасибо огромное.

eganeev

Дата 3.2.2007, 19:24 (ссылка)

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

Загрузка ...

Новичок

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

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

А нет такой программы еще на Visual Basic?

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