Вылетает программа

@ElKros · Регистрация: 14.02.2018

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

вылетает при нахождении у0 и дальше не идет

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
void func(double y0, double y1, double y2, double y3, double t0, double max_t, double h, int n, double matr[N][N]) {
 
    double t = t0;
    for (int i = 0; i < ((max_t - t0) / h); i++) {
        std::cout << t << "\t" << std::setprecision(5) << y0 << "\t" << y1 << "\t" << y2 << "\t" << y3 << std::endl;
        y0 += matr_inv[N][N]*h*dy0dt(y0,y1,y2,y3,t0);
        y1 += (matr_inv[N][N] * h*dy1dt(y0, y1, y2, y3, t0));
        y2 += (matr_inv[N][N] * h*dy2dt(y0, y1, y2, y3, t0));
        y3 += (matr_inv[N][N] * h*dy3dt(y0, y1, y2, y3, t0));
        t = t + h;
    }
    std::cout << t << "\t" << std::setprecision(5) << y0 << "\t" << y1 << "\t" << y2 << "\t" << y3 << std::endl;
}

@ElKros · 06.02.2019, 00:16 **[ТС]**

zss, решила проверить на системе, которая имеет точное решение. но у меня не сходится значения, которые в программе и мои отдельно..
я брала вот такую систему :
$https://www.cyberforum.ru/cgi-bin/latex.cgi?{{u}_{1}}^{'}=-10{u}_{1}$
$https://www.cyberforum.ru/cgi-bin/latex.cgi?{{u}_{2}}^{'}=-{u}_{2}$
ее точное решение:
$https://www.cyberforum.ru/cgi-bin/latex.cgi?{u}_{1}(t)={u}_{1}(0){e}^{-10t}$
$https://www.cyberforum.ru/cgi-bin/latex.cgi?{u}_{2}(t)={u}_{2}(0){e}^{-t}$
где $https://www.cyberforum.ru/cgi-bin/latex.cgi?{u}_{1}(0)={u}_{2}(0)=1$

в программе получилось, что

C++
1
2
3
4
5
6
7
8
9
10
0       1       1
0.1     0.33333 0.90476
0.2     0.11111 0.81859
0.3     0.037037        0.74063
0.4     0.012346        0.6701
0.5     0.0041152       0.60628
0.6     0.0013717       0.54854
0.7     0.00045725      0.4963
0.8     0.00015242      0.44903
0.9     5.0805e-05      0.40626

но u1 не совпадает совсем, а u2 хоть как-то похоже.
в чем может быть ошибка ?

@ElKros · 07.02.2019, 12:05 **[ТС]**

zss, подскажите, пожалуйста...?

Новые блоги и статьи Все статьи Все блоги /
Рефакторинг программы уравнивания. Massaraksh7 26.05.2026 Пример по предыдущей записи в блоге. Но, надо заметить, что, во-первых, там оптимизация не только математики, но и работы с базой данных, и с графами, а во-вторых, это ещё не всё.	Использование TThread в Lazarus для математических вычислений. Massaraksh7 25.05.2026 Производя рефакторинг своих программ на предмет ускорения их работы, обратил внимание на такой аспект, как сокращение времени матвычислений. Дело в том, что приходится работать с большими матрицами. . .	Модель здравосохранения 18. Чем здоровее работник, тем быстрее выгорает anaschu 24.05.2026 Имитационная модель корпоративного здравоохранения: что показывает математика Сегодня в модели рабочего коллектива на AnyLogic появились три новые механики — выгорание через накопленную усталость,. . .	Модель здравосохранения 17. Планы на выгорание anaschu 23.05.2026 Вот конкретная схема реализации: В классе Работник добавить: накопленнаяУсталость — растёт каждый час работы, снижается в перерывы и болезни коэффициентПрезентеизма — снижает продуктивность. . .
Изменение цветов в палитре gif файла aka фавикона russiannick 23.05.2026 Изменение цветов в палитре gif файла, юзаемого как фавиконка в составе html-файла, помещенная в base64, средствами нативного Java Script, навеянное сном в майский день. Для работы необходим браузер,. . .	Модель здравосохранения 16. Слишком хорошие и здоровые сотрудники уходят, недовольные зарплатой anaschu 23.05.2026 Отладка увольнений и настройка производительности Сегодня во второй половине дня разобрались с механикой увольнений и настроили коэффициент сложности заданий. Вот что было сделано. . . .	Как я стал коммунистом))) Модель сохранения здоровья сотрудников, запись блога номер 15 anaschu 23.05.2026 Внезапно хорошее здоровье сотрудников не нужно капиталистам?))	Модель здравоСохранения 15. Как мы чинили AnyLogic модель рабочего коллектива: сочленение диаграммы состояний болезней и поломок в ресурспул anaschu 23.05.2026 Как мы чинили AnyLogic модель рабочего коллектива Сегодня разобрались с пятью багами, из-за которых модель либо падала с ошибкой, либо давала совершенно бессмысленные результаты. Каждый баг был. . .

@ElKros 5 / 4 / 1 Регистрация: 14.02.2018 Сообщений: 583
	07.02.2019, 12:05 [ТС]
	zss, подскажите, пожалуйста...? 0