Представить алгоритм решения дифференциального уравнения методом Эйлера и модифицированным методом Эйлера в виде блок-схемы.
Представить алгоритм решения дифференциального уравнения методом Эйлера и модифицированным методом Эйлера в виде блок-схемы. Числовой результат можно получить с использованием решателей МАТЛАБа. Вариант 1 y'=x+cos(y/e), y1,4=2.5, [1.4;2.4]
Блок-схема метода Эйлера:
Блок-схема модифицированного метода Эйлера:
Программа в MATLAB:
clear all
close all
a=1.4; b=2.4; y0=2.5;
h=0.1;%шаг
n=(b-a)/h;% число шагов
%Решение методом Эйлера
[Xe,Ye]= Euler(a,b,h,y0);
%выводим таблицу с решением
fprintf('X\t\t\t Y\t\n');
for i=1:n+1
fprintf('%f\t %f\t\n',Xe(i),Ye(i));
end
%строим график
plot(Xe, Ye);
xlabel('x')
ylabel('y')
grid on
hold on
%Решение модифицированным методом Эйлера
[XM,YM]= EulerM(a,b,h,y0);
%выводим таблицу с решением
fprintf('X\t\t\t Y\t\n');
for i=1:n+1
fprintf('%f\t %f\t\n',XM(i),YM(i));
end
%строим график
plot(XM, YM, '*-');
legend('Решение уравнения методом Эйлера','Решение модифицированным методом Эйлера','Location','northwest');
%Решение дифференциального уравнения Методом Эйлера на отрезке [a, b]
%начальное условие y0 шагом h
function [X,Y] = Euler(a,b,h,y0)
X=a:h:b;
n=(b-a)/h;
Y=zeros([1,n+1]);
Y(1)=y0;
y=y0;
for i = 1:n
Y(i+1) = Y(i)+h*f(X(i),Y(i));
end
end
%Решение дифференциального уравнения модифицированным методом Эйлера
%начальное условие y0 шагом h
function [X,Y] = EulerM(a,b,h,y0)
X=a:h:b;
n=(b-a)/h;
Y=zeros([1,n+1]);
Y(1)=y0;
for i = 1:n
y1=Y(i)+h*f(X(i),Y(i));
Y(i+1) = Y(i)+(h/2)*(f(X(i),Y(i))+f(X(i)+h,y1));
end
end
%Правая часть дифференциального уравнения y'=f(x,y)
function Z = f(x,y)
Z=x+cos(y/exp(1));
end
Результаты:
Расчётная таблица по методу Эйлера:
X Y
1.400000 2.500000
1.500000 2.700606
1.600000 2.905182
1.700000 3.113304
1.800000 3.324579
1.900000 3.538658
2.000000 3.755234
2.100000 3.974053
2.200000 4.194914
2.300000 4.417671
2.400000 4.642237
Расчётная таблица по модифицированному методу Эйлера:
X Y
1.400000 2.500000
1.500000 2.702591
1.600000 2.908878
1.700000 3.118457
1.800000 3.330962
1.900000 3.546069
2.000000 3.763507
2.100000 3.983052
2.200000 4.204537
2.300000 4.427851
2.400000 4.652940

- Представить данную выборку в виде вариационного ряда. Построить полигон частот, гистограмму, кумуляту и график
- Представить данную функцию в виде интеграла Фурье: fx=-1, если 1≪х≪22, если 0≪х<11, если-1≪х<0 , в
- Представить заданную ФАЛ: 1.1. таблицей истинности; 1.2. совершенной дизъюнктивной нормальной формой; 1.3. совершенной конъюнктивной нормальной формой. 2. Минимизировать
- Представить заданную ФАЛ: 4.1.1. таблицей истинности; 4.1.2. совершенной дизъюнктивной нормальной формой; 4.1.3. совершенной конъюнктивной нормальной формой. 4.2. Минимизировать
- Представить заданные десятичные числа B1, B2, B3 в ЭВМ в форме с фиксированной запятой
- Представить запись результата измерений напряжения на ванне, при условии получения результата измерений 35 В
- Представить запись результата измерений напряжения на ванне, при условии получения результата измерений 35 В. 2
- Представитель ООО «Новая эра» — резидент инновационного центра «Сколково» — обратился в юридическую консультацию
- Представитель ответчика в жалобе на решение суда, по делу о разделе имущества между супругами
- Представитель Президента РФ в Поволжском федеральном округе потребовал от руководителя одного из деревообрабатывающих предприятий
- Представитель Рязанской областной психиатрической больницы обратился в суд с административным исковым заявлением о недобровольной
- Представитель таможенного перевозчика АО «Line&Road» Лошадкин обратился к начальнику пограничного таможенного поста в целях
- Представителю коммерческой организации специалистом экспедиции арбитражного суда было разъяснено, что подаваемые документы должны быть
- Представить 5 определений понятия «процесс» из различных областей знаний, сформированных деятельностью людей. Выбрать наиболее