Решить задачу Коши для дифференциального уравнения второго порядка методом Эйлера-Коши

@stolos · Регистрация: 01.11.2012

Решить задачу Коши для дифференциального уравнения второго порядка методом Эйлера- Коши. y//=2x-y+y/

на отрезке [0,1] c шагом h=0,2;

Начальные условия: у0=1 z0=1

Написал программу, но сомневаюсь в ее правильности, если не сложно проверьте!!

Pascal

program rabota1;
 
function fxyz(x, y, z: Real): Real;
begin
  fxyz := 2*x-y+y;
end;
 
 
 
var
  x, x0, a, b, h, y, y0, z0, s, z, k: Real;
  i, n: Integer;
 
 
begin
  writeln('vvedite a', ' b', ' n', ' y0 ', ' z0 ' );
  readln(a, b, n, y0, z0);
  Writeln('x':3, 'y':6);
  x0 := a;
  h := (b - a) / n;
  for i := 1 to n do
  begin
    x := x0 + h;
    s := y0 + +h * z0;
    k := z0 + h * fxyz(x0, y0, z0);
    y := y0 + h / 2 * (z0 + k);
    z := z0 + h / 2 * (fxyz(x0, y0, z0) + fxyz(x, s, k));
    x0 := x;
    y0 := y;
    z0 := z;
    writeln(x,'    ', y);
  end;
end.

@Крюгер · 24.09.2013, 16:11

Не по теме:

код в теги следует закидывать...
а вы пробовали этот код запустить? 0_о

@stolos · 24.09.2013, 16:56 **[ТС]**

Да) Но значения которые она выдает не сходятся на 1 десятую с решением. Причем чем больше шагов, тем больше отклонение. В чем проблема?

@Крюгер
	24.09.2013, 16:11 #2
	Не по теме: код в теги следует закидывать... а вы пробовали этот код запустить? 0_о 0

@stolos 0 / 0 / 1 Регистрация: 01.11.2012 Сообщений: 26
	24.09.2013, 16:56 [ТС]	3
	Да) Но значения которые она выдает не сходятся на 1 десятую с решением. Причем чем больше шагов, тем больше отклонение. В чем проблема? 0

Опции темы