метод Эйлера в С

@marina11 · Регистрация: 15.11.2012

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

Здравствуйте!Нужно написать код решения дифф. уравнения простым методом Эйлера.
Уравнение: y''+(y'/x)-(y/x^2)=8*x
Нач. условия:
h=0.1
x₀=1
y(1)=4
y'(1)=2

Вот код, который попыталась написать, выдаёт не совсем верные значения 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
#include "stdafx.h"
#include <cstdlib>
#include <iostream>
#include <math.h>
 
using namespace std;
 
double myfunc(const double& x, const double& y,const double& z)
{
    return (8*x-(z/x)+(y/x*x));
}
 
int _tmain(int argc, _TCHAR* argv[])
{
    int i, n; 
    double x, y, z, h;
 
    cout<< " h = "; cin>>h;
    cout<< " n = "; cin>>n;
    cout<< " x0 = "; cin>>x;
    cout<< " y0 = "; cin>>y;
    cout<< " z0 = "; cin>>z;
 
    for (i = 0; i <= n; i++)
    {
 
        y += h * z;
        z += h * myfunc(x,y,z); 
        
        cout << "   x = " << x;
        cout << "   y = " << y;
        cout << "   z = " << z<< "\n";
 
        x += h;
    }
 
    system("PAUSE");
    return EXIT_SUCCESS;
}

Помогите пожалуйста найти ошибку. Она скорее всего в цикле в функции main.
Ниже приложу скрин из маткада с правильным решением.

@Том Ардер · 23.04.2013, 01:08

Строка 11 кода:

Сообщение от marina11

y/x*x

Вместо желаемого y/x² вычисляется (по правилам С) y/x*x=(y/x)*x=y
Пренебрежение скобками дорого обходится.

P.S. Это тема для ветки форума, посвящённой программированию на С++. Математика здесь ни при чём.

Новые блоги и статьи Все статьи Все блоги /
Знаешь почему 90% людей редко бывают счастливыми? kumehtar 14.04.2026 Потому что они ждут. Ждут выходных, ждут отпуска, ждут удачного момента. . . а удачный момент так и не приходит.	Фиксация колонок в отчете СКД Maks 14.04.2026 Фиксация колонок в СКД отчета типа Таблица. Задача: зафиксировать три левых колонки в отчете. Процедура ПриКомпоновкеРезультата(ДокументРезультат, ДанныеРасшифровки, СтандартнаяОбработка) / / . . .	Настройки VS Code Loafer 13.04.2026 { "cmake. configureOnOpen": false, "diffEditor. ignoreTrimWhitespace": true, "editor. guides. bracketPairs": "active", "extensions. ignoreRecommendations": true, . . .	Оптимизация кода на разграничение прав доступа к элементам формы Maks 13.04.2026 Алгоритм из решения ниже реализован на нетиповом документе, разработанного в конфигурации КА2. Задачи, как таковой, поставлено не было, проделанное ниже исключительно моя инициатива. Было так:. . .
Контроль заполнения и очистка дат в зависимости от значения перечислений Maks 12.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: реализовать контроль корректности заполнения дат назначения. . .	Архитектура слоя интернета для сервера-слоя. Hrethgir 11.04.2026 В продолжение https:/ / www. cyberforum. ru/ blogs/ 223907/ 10860. html Знаешь что я подумал? Раз мы все источники пишем в голове ветки, то ничего не мешает добавить в голову такой источник, который сам. . .	Подстановка значения реквизита справочника в табличную часть документа Maks 10.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: при выборе сотрудника (справочник Сотрудники) в ТЧ документа. . .	Очистка реквизитов документа при копировании Maks 09.04.2026 Алгоритм из решения ниже применим как для типовых, так и для нетиповых документов на самых различных конфигурациях. Задача: при копировании документа очищать определенные реквизиты и табличную. . .