ОДУ методом Эйлера (проверить алгоритм Делфи)

[xenon1000]

Что не так в данном алгоритме? Почему при данном условии решение выводит все 0, и почему решений 11 всесто 10?
Условие задачи:
Используя метод Эйлера, составить таблицу приближенных значений интеграла дифференциального уравнения y’=3*y^(2/3) на отрезке [0;1] с шагом интегрирования h=0,1 и начальным условием y(0)=0.
Алгоритм в ДЕлфи:

procedure TForm1.Button1Click(Sender: TObject);
var
x,y,a,b,h:real;
function f(x,y:real):real;
 begin
   f:=power(3*y,2/3);
 end;
begin
 ListBox1.Clear;
h:=StrToFloat(Edit1.Text);
y:=StrToFloat(Edit2.Text);
a:=StrToFloat(Edit3.Text);
b:=StrToInt(Edit4.Text);
x:=a;
repeat
y:=y+h*f(x,y);
x:=x+h;
ListBox1.Items.Add(FormatFloat('y(0.###) = ', x) + FormatFloat('0.###', y));
until  (x>b);
   end;
end.

[alexander13]

А как Вы сами думаете, что там может быть кроме нуля с таким y(0)?
При цикле с таким постусловием он и должен выводить 11 строк, посчитайте сами.

Да, сам алгоритм правильный.

Upd. Я имею ввиду, сам метод Эйлера реализован правильно. А вот этот порядок операции неверный. Угадайте, почему.

y:=y+h*f(x,y);
x:=x+h;
ListBox1.Items.Add(FormatFloat('y(0.###) = ', x) + FormatFloat('0.###', y));

Upd. 2 А, да.. 3*y^(2/3) это не совсем то же, что power(3*y,2/3).

[xenon1000]

Просто я не превык решая 1 задачу получать нулевой ответ) при этом этот ответ надо скриншотить.
to Upd. Как не переставлял местами все равно 11 ответов, если поставить вывод в листбокс перед Х то значения у(0) выдает вместо у(1,1). В чем косяк?
to Upd2. f:=3*Power(y,2/3); так да? Иначе 3 тоже возводит в степень.

И последний вопрос. Я скачал готовую прогу написано решает методом эйлера тоже, но там еще навороты - график, оформлено и т.д. Вставляю свою формулу - результат отличается от моего, при у(0)=1 как такое возможно?

[alexander13]

Сообщение от xenon1000

Просто я не превык решая 1 задачу получать нулевой ответ)

Вы подумайте.. и увидите, что там никакого другого ответа и не могло быть при таких начальных данных.

Сообщение от xenon1000

Как не переставлял местами все равно 11 ответов, если поставить вывод в листбокс перед Х то значения у(0) выдает вместо у(1,1). В чем косяк?

Ответов и ДОЛЖНО быть 11. Я же говорю, посчитайте. Просто если оставить Ваш порядок операций, то последним будет что-то типа
y(1,1) = ... При том, что у Вас функция вычисляется на отрезке [0,1] - как такое может быть?

Сообщение от xenon1000

to Upd2. f:=3*Power(y,2/3); так да? Иначе 3 тоже возводит в степень.

Да.

Сообщение от xenon1000

И последний вопрос. Я скачал готовую прогу написано решает методом эйлера тоже, но там еще навороты - график, оформлено и т.д. Вставляю свою формулу - результат отличается от моего, при у(0)=1 как такое возможно?

Я не телепат, не знаю, что Вы там скачали. Ради проверки я тоже нашел исходники другой проги, она мне выдала те же самые результаты. Может, в той, что Вы скачали, ошибка. Может еще что-то, не знаю.

[xenon1000]

Спасибо!

[ViktorR]

Обратите внимание на порядок вывода сообщения и увеличения x:
У Вас:

repeat
y:=y+h*f(x,y);
x:=x+h;
ListBox1.Items.Add(FormatFloat('y(0.###) = ', x) + FormatFloat('0.###', y));
until  (x>b);

Но должно быть так:

repeat
y:=y+h*f(x,y);
ListBox1.Items.Add(FormatFloat('y(0.###) = ', x) + FormatFloat('0.###', y));
x:=x+h;
until  (x>b);

Выводиться должно одиннадцать сообщений: 0.0, 0.1, 0.2, ... 0.9, 1.0 - счет с нуля ...

Как-то так ...

[xenon1000]

Возник вопрос. Я хочу вывести ответ не в листбоксе перечнем, как выше написано в алгоритме, а раскидать значения в таблицу, используя объект StringGrid. Я прочитал кучу всего, но не нашел как циклом заполнить таблицу значениями.

 
For i:=1 to 10 do
begin
StringGrid1.Cells[i, 1]:=(FormatFloat('0.###', y));
end;

Этот цикл просто заполняет всеми значениями каждую ячейку, а надо раскидать по порядку.
Вкратце как правильно составить цикл вывода значений по ячейкам? Или придется переписывать полностью цикл repeat?