Форум программистов, компьютерный форум CyberForum.ru

Изменить код С++ - C++

Восстановить пароль Регистрация
 
Виталий32
1 / 1 / 0
Регистрация: 01.10.2012
Сообщений: 13
11.06.2013, 14:18     Изменить код С++ #1
Вот код
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
#include <cstdlib>
#include <iostream>
#include <math.h>
 
using namespace std;
 
float f1(float x, float y, float z)
{
    return cos(3*x) - 4*x;
}
 
float f2(float x, float y, float z)
{
    return z;
}
 
int main()
{
    // начальные значения
    float xk=0.8, x=0, y=10, z=5, h=0.1;
    float k1, k2, k3, k4, q1, q2, q3, q4;
 
    for (x=0; x<=xk; x+=h)
    {
        // выодим на экран x y z
        cout << "\nx = " << x << "\t\ty = " << y << "\t\tz = " << z;
        
        k1 = h * f1(x, y, z);
        q1 = h * f2(x, y, z);
        
        k2 = h * f1(x + h/2, y + q1/2, z + k1/2);
        q2 = h * f2(x + h/2, y + q1/2, z + k1/2);
        
        k3 = h * f1(x + h/2, y + q2/2, z + k2/2);
        q3 = h * f2(x + h/2, y + q2/2, z + k2/2);
        
        k4 = h * f1(x + h, y + q3, z + k3);
        q4 = h * f2(x + h, y + q3, z + k3);
        
        float new_y = (k1 + 2*k2 + 2*k3 + k4)/6;
        float new_z = (k1 + 2*k2 + 2*k3 + k4)/6;
        
        y += new_y;
        z += new_z;
    }
 
    cin.get();
    return 0;
}
Программа численно решает ДУ методом Рунге-Кутты

он работает для дифференциального уравнения y"+4x=cos(3x) , надо изменить для уравнения y"+4y=cos(3x)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
11.06.2013, 14:18     Изменить код С++
Посмотрите здесь:

Необходимо изменить код программы. C++
C++ нужно изменить код программы
C++ Изменить код программы
C++ Изменить код. ООП
Изменить код C++
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Expocooky
 Аватар для Expocooky
56 / 56 / 2
Регистрация: 06.06.2013
Сообщений: 112
11.06.2013, 14:20     Изменить код С++ #2
Тыкну наугад
Может здесь нужно поменять 4х на 4у
C++
1
    return cos(3*x) - 4*x;
Виталий32
1 / 1 / 0
Регистрация: 01.10.2012
Сообщений: 13
11.06.2013, 14:25  [ТС]     Изменить код С++ #3
Цитата Сообщение от Expocooky Посмотреть сообщение
Тыкну наугад
Может здесь нужно поменять 4х на 4у
C++
1
    return cos(3*x) - 4*x;
менял считает не правильно тогда
Fedorys
 Аватар для Fedorys
487 / 198 / 8
Регистрация: 19.03.2013
Сообщений: 444
11.06.2013, 14:44     Изменить код С++ #4
Ну насколько я понимаю, достаточно поменять функцию f1:
C++
1
2
3
4
float f1(float x, float y, float z)
{
    return cos(3*x) - 4*y;
}
Добавлено через 8 минут
И еще меня смущает то, что у вас в уравенении вторая производная, а метод используется для уравнения вида y' = f(x,y).

Добавлено через 7 минут
Да, пожалуй, что здесь надо менять обе функции. Попробуйте вот так:
C++
1
2
3
4
5
6
7
8
9
float f1(float x, float y, float z)
{
    return z;
}
 
float f2(float x, float y, float z)
{
    return cos(3*x) - 4*y;
}
Виталий32
1 / 1 / 0
Регистрация: 01.10.2012
Сообщений: 13
11.06.2013, 15:05  [ТС]     Изменить код С++ #5
Цитата Сообщение от Fedorys Посмотреть сообщение
Ну насколько я понимаю, достаточно поменять функцию f1:
C++
1
2
3
4
float f1(float x, float y, float z)
{
    return cos(3*x) - 4*y;
}
Добавлено через 8 минут
И еще меня смущает то, что у вас в уравенении вторая производная, а метод используется для уравнения вида y' = f(x,y).

Добавлено через 7 минут
Да, пожалуй, что здесь надо менять обе функции. Попробуйте вот так:
C++
1
2
3
4
5
6
7
8
9
float f1(float x, float y, float z)
{
    return z;
}
 
float f2(float x, float y, float z)
{
    return cos(3*x) - 4*y;
}
спасибо теперь правильно считает
Yandex
Объявления
11.06.2013, 15:05     Изменить код С++
Ответ Создать тему
Опции темы

Текущее время: 02:57. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru