Решение ОДУ 2го порядка

17.05.2011, 20:15. Показов 2022. Ответов 0

Метки нет (Все метки)

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

Здравствуйте нужна ваша помощь
ОДУ 2го порядка. y''=5y=0, y(0)=0 y(1)=1;
Необходимо решить методом рунге кутта 4го п.т.

Я сделал как систему ду.
z' = -5y;
y'=z;

Не получается реализовать

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
#include "stdafx.h" 
#include "iostream" 
 
using namespace std; 
 
double y(double x, double y,double z) 
{ 
return -5*y;
} 
double z(double x, double y,double z) 
{ 
return z; 
} 
 
 
 
void r_k_4(double xn,double xk,double h,double y0,double z0,int n) 
{ 
int i; 
double x; 
double *yt = new double [n]; 
double *zt = new double [n]; 
double k1,k2,k3,k4; 
double k_1,k_2,k_3,k_4; 
 
for(x=xn,i=0; x<=xk,i<=n; x+=h,i++) 
{ 
k1 = h * y(x, xn, z0); 
k_1 = h * z(x, xn, z0); 
 
k2 =h*y(x+h/2, y0+k1/2, z0+k_1/2); 
k_2=h*z(x+h/2, y0+k1/2, z0+k_1/2); 
 
k3=h*y(x+h/2,y0+k2/2,z0+k_2/2); 
k_3=h*z(x+h/2,y0+k2/2,z0+k_2/2); 
 
k4=h*y(x+h,y0+k3,z0+k_3); 
k_4=h*z(x+h,y0+k3,z0+k_3); 
 
yt[i]=y0+(k1+2*k2+2*k3+k4)/6; 
zt[i]=z0+(k_1+2*k_2+2*k_3+k_4)/6; 
} 
 
for(x=xn,i=0; x<=xk,i<=n; x+=h,i++) 
{ 
 printf("y[%.2lf] = %3.16lf z[%.2lf] = %3.16lf\n",x,yt[i],x,zt[i]); 
} 
} 
 
void main() 
{ 
double xn=0.0,xk=1.0,y0=0,z0=1,n,h; 
cout<<"N = "; cin>>n; 
h=(xk-xn)/n; 
r_k_4(xn,xk,h,y0,z0,n);
 
}

cpp_developer

Эксперт

20123 / 5690 / 1417

Регистрация: 09.04.2010

Сообщений: 22,546

Блог

17.05.2011, 20:15

Ответы с готовыми решениями:

С++ .Задач условий Коши Оду 2-го порядка
Задача Коши Оду 2-го порядка. Заданы условия для Коши y(1)=3, y'(1)=2.Собственно сам глупый вопрос, как их задать, если с y(0) все понятно...

Численно решить ОДУ методом Рунге-Кутта 4 порядка
Доброго времени суток. Мне нужно численно решить ОДУ методом Рунге-Кутта 4 порядка. Уравнение имеет вид: p'=-A-B*p, p(0)=p0. Проблема...

Решение ОДУ в FreePascal
Здравствуйте, буду признателен если не откажите в помощи. Нужно в паскале задать условия программы, чтобы она решала систему Обыкновенных...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.

raxper

Эксперт

30234 / 6612 / 1498

Регистрация: 28.12.2010

Сообщений: 21,154

Блог

17.05.2011, 20:15

Помогаю со студенческими работами здесь

Решение системы ОДУ на С++
В руки попал код на крестах, в которых я не разбираюсь, Эйлер явный и неявный работают исправно, а вот Рунге 4 порядка работать...

Решение системы ОДУ + СЛАУ
Добрый день. Помогите, пожалуйста, с решением системы ОДУ + СЛАУ есть система уравнений 1) dx/dt=f1(x,y), dx/dt(0) = a1 2) y =...

Решение СЛАУ третьего порядка
Разработать алгоритм и программу для решения системы алгебраических уравнений третьего порядка методом обратной матрицы....

Решение системы 3(k+1) о.д.у первого порядка
Здравствуйте, решаю систему 3(k+1) ОДУ первого порядка.Решаю 2 методами - Р.К. 4 порядка и Мерсона. Сама система: ...

Решение системы линейных уравнений третьего порядка по правилу Крамера
Напишите функцию для решения системы линейных уравнений третьего порядка по правилу Крамера

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:

Новые блоги и статьи Все статьи Все блоги /
Подключение Box2D v3 к SDL3 для Android: физика и отрисовка коллайдеров 8Observer8 29.01.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами. Версия v3 была полностью переписана на Си, в. . .	Инструменты COM: Сохранение данный из VARIANT в файл и загрузка из файла в VARIANT bedvit 28.01.2026 Сохранение базовых типов COM и массивов (одномерных или двухмерных) любой вложенности (деревья) в файл, с возможностью выбора алгоритмов сжатия и шифрования. Часть библиотеки BedvitCOM Использованы. . .	Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 28.01.2026 Содержание блога SDL3 имеет собственные средства для загрузки и отображения PNG-файлов с альфа-каналом и базовой работы с ними. В этой инструкции используется функция SDL_LoadPNG(), которая. . .	Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL3_image 8Observer8 27.01.2026 Содержание блога SDL3_image - это библиотека для загрузки и работы с изображениями. Эта пошаговая инструкция покажет, как загрузить и вывести на экран смартфона картинку с альфа-каналом, то есть с. . .
Влияние грибов на сукцессию anaschu 26.01.2026 Бифуркационные изменения массы гриба происходят тогда, когда мы уменьшаем массу компоста в 10 раз, а скорость прироста биомассы уменьшаем в три раза. Скорость прироста биомассы может уменьшаться за. . .	Воспроизведение звукового файла с помощью SDL3_mixer при касании экрана Android 8Observer8 26.01.2026 Содержание блога SDL3_mixer - это библиотека я для воспроизведения аудио. В отличие от инструкции по добавлению текста код по проигрыванию звука уже содержится в шаблоне примера. Нужно только. . .	Установка Android SDK, NDK, JDK, CMake и т.д. 8Observer8 25.01.2026 Содержание блога Перейдите по ссылке: https:/ / developer. android. com/ studio и в самом низу страницы кликните по архиву "commandlinetools-win-xxxxxx_latest. zip" Извлеките архив и вы увидите. . .	Вывод текста со шрифтом TTF на Android с помощью библиотеки SDL3_ttf 8Observer8 25.01.2026 Содержание блога Если у вас не установлены Android SDK, NDK, JDK, и т. д. то сделайте это по следующей инструкции: Установка Android SDK, NDK, JDK, CMake и т. д. Сборка примера Скачайте. . .

Наверх