Delphi. Считаем вторую производную+график.

[mishka-medved]

Проблема вот в чем. Задали мне, на мой взгляд, умопомрачительную курсовую на этот семестр. Нужно найти вторую производную от очень сложной функции+нарисовать график этой функции НЕ теми возможностями, которые предусмотрены в делфи. Сначала я делал всё по муторному, но препод сказал чтобы я эту функцию оформил в отдельном модуле в виде функции, а не считал сначала y потом y2 и потом y3. На данный момент это меня завело в ступор. КАК это сделать(оформить в виде функции)? Вот дана функция
y:=power(x,1/3)*(6*sin(12*x)+12*cos(6*x));
Вторую производную выводим по формуле
z:=(y2-2*y+y3)/(dx*dx)

где
dx-приращение функции
y2, y3 это Y от приращения x+dx и x-dx

Изначально я успешно считал это всё вот так:

PHP код:

x:= strtofloat (edit1.Text);
dx:=strtofloat (edit2.Text);
y:=power(x,1/3)*(6*sin(12*x)+12*cos(6*x));
x2:=x+dx;
y2:=power(x2,1/3)*(6*sin(12*x2)+12*cos(6*x2));
x3:=x-dx;
y3:=power(x3,1/3)*(6*sin(12*x3)+12*cos(6*x3));
z:=(y2-2*y+y3)/(dx*dx);
label4.Caption:=floattostr(z); 


Но что и как делать теперь?

[s-andriano]

Оформите Вашу функцию (математическую) отдельной функцией (программной) и при вычислении производных обращайтесь к ней столько раз, сколько необходимо.

[mishka-medved]

Это понятно, но сколько бы я ни читал информацию насчет этого, не могу сообразить как это делать.

[s-andriano]

Примерно так:
описываем функцию:

function MyFunc(x : extended) : extended;
begin
  MyFunc := power(x,1/3)*(6*sin(12*x)+12*cos(6*x));
end;

а потом пользуемся ею:

x:= strtofloat (edit1.Text);
dx:=strtofloat (edit2.Text);
y := MyFunc(x);
y2 := MyFunc(x+dx);
y3 := MyFunc(x-dx);

[mishka-medved]

Спасибо, сейчас попробую

[mishka-medved]

Запихнул в отдельный модуль эту функцию

unit Unit2;

{$mode objfpc}{$H+}

interface
function proizv(x:real):real;

implementation
uses
  Classes, SysUtils, Math, Unit1;

 function proizv(x:real):real;
 begin
y:=x*ln(x);
 end;

end.

Считаю вторую производную в программе для простейшей функции. Выводит полнейшую ахинею. Вторая производная логарифма это 1/x. А он мне сумасшедшие числа начинает выводить. Вот код основной программы

procedure TForm1.Button1Click(Sender: TObject);
begin
x:=strtofloat (edit1.Text);
dx:=strtofloat (edit2.Text);
y:=proizv(x);
y2:=proizv(x+dx);
y3:=proizv(x-dx);
z:=(y2-2*y+y3)/(dx*dx);
label1.Caption:=floattostr(z);
end;

Что опять не так сделал?

[s-andriano]

Ваша функция proizv вместо результата возвращает мусор.
В коде должно быть:

proizv := ...

[mishka-medved]

Вместо чего? Я извиняюсь, мне пока не показать как это должно быть на деле я не пойму.

Добавил: Кажется получается. Сейчас попробую для курсовой править...

[s-andriano]

Возьмите любой учебник по программированию и прочитайте, что такое функция.
Любое общение подразумевает некоторый минимальный уровень познаний. Иначе мы так и будем разговаривать на разных языках.