поправить прогу,метод Эйлера

@lekz · Регистрация: 23.03.2013

Author24 — интернет-сервис помощи студентам

Помогите поправить код,в долгу не останусь)
нужно методом Эйлера получить значения в пределах [0;2] с шагом h=0,1 . y(0)=1

ур-е
y'-2y/(1+x)=((1+x)^2)*e^x

точное решение{ y(x)=((1+x)^2)*e^x }

Код Pascal

Pascal

program rab1;
 
function fxy(x, y: Real): Real;
begin
  fxy :=exp(ln(exp(x)*x))*((1+x)*(1+x))+2*y/(1+x);
end;
 
var
  x, t, a, b, h: Real;
  i, n: Integer;
 
begin
  writeln('vvedite a', ' b', ' h');
  readln(a, b, h);
  n := trunc((b - a) / h) + 1;
  
  Writeln('x':3, ' Euler':14);
  x := a;
  t := 1;
  for i := 1 to n do
  begin
    Writeln
    (x:5:2,  t:18:10);
    t := t + h * fxy(x, t);
    x := x + h;
  end;
end.

Добавлено через 20 часов 45 минут
del

Добавлено через 32 секунды
тему можно удалить!

@ermolay · 31.10.2013, 09:10

lekz, если у Вас есть верное решение - выложите его, наверняка для кого-нибудь оно будет полезным.

@lekz · 10.11.2013, 18:50 **[ТС]**

Pascal

program eiler;
function fxy(x, y: Real): Real;
 begin
  fxy :=exp(ln(2.71)*x)*((1+x)*(1+x))+(2*y)/(1+x);
   end;
   const
   yo=1;
 var
  x, m,a,b,h: Real;
  i,n: integer;
 
 begin
  writeln('vvedite a=   ' , 'vvedite b=   ' , 'vvedite h=  ' );
  readln(a,b,h);
 
  n:= trunc((b-a)/h)+1;
  Writeln('x':4, 'y':4);
  x:= a;
  m:= 1;
 for i := 1 to n do
  begin
Writeln
(x:5:2 , m:5:2);
 
m:= m + h * fxy(x,m);
x := x + h;
end;
end.

@ermolay 3451 / 2389 / 2135 Регистрация: 04.12.2011 Сообщений: 3,966
	31.10.2013, 09:10	2
	lekz, если у Вас есть верное решение - выложите его, наверняка для кого-нибудь оно будет полезным. 1

	10.11.2013, 18:50