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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
lostman23
0 / 0 / 0
Регистрация: 03.02.2014
Сообщений: 13
#1

Подскажите правильно ли реализована программа решения уравнения? - C++

07.02.2014, 15:31. Просмотров 173. Ответов 0
Метки нет (Все метки)

необходимо решить задачу такую:
http://www.cyberforum.ru/cgi-bin/latex.cgi?\frac{du}{dt}=u(1-u-rv)+\frac{{d}^{2}u}{{dt}^{2}}\frac{dv}{dt}=-buv+\frac{{d}^{2}v}{{dt}^{2}}

вот код. использую явную схему. скажите правильно ли это?
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
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
#include <iostream>
#include <cmath>
#include <fstream>
 
#define a -30.0 //левая граница отрезка
#define b 30.0 //правая граница отрезка
#define h 0.5 // шаг по пространству
#define dt 0.0000625 // шаг по времени
 
 
 
#define T 11
#define X (int)((b - a) / h + 1) 
 
using namespace std;
 
void init(float u[T][X],float v[T][X])
{
        double x=a;
        for(int j = 0; j < X; j++) 
        {
            u[0][j] =-6.0 / (cosh(x) * cosh(x));
            v[0][j] = 1+6.0 / (cosh(x) * cosh(x));
            x += h;
        //  cout<<u[0][j]<<"    "<<v[0][j]<<endl;
        }   
}
 
 
int main()
{
    
    // -- parameters --
    float r = 5,a2,b2;
    float b1 = 1.25f;
 
    float u[T][X], v[T][X];
 
 
    init(u,v);
 
    for(int i = 0; i < T-1; i++) 
    {
        for(int j = 0; j < X; j++) 
        {
            a2=u[i][j]*(1-u[i][j]-r*v[i][j])*dt+dt*(u[i][(j + 1 + X) % X]-2*u[i][j]+u[i][(j - 1 + X) % X])/h/h+u[i][j];
            b2=(-b1*u[i][j]*v[i][j])*dt+dt*(v[i][(j + 1 + X) % X]-2*v[i][j]+v[i][(j - 1 + X) % X])/h/h+v[i][j];
            u[i+1][j] =a2;
            v[i+1][j] =b2;
        }
    }
 
    ofstream out;
    out.open("der.xls");
 
    ofstream out1;
    out1.open("de1.txt");
 
    ofstream out2;
    out2.open("de2.txt");
    //out1<<"<table> <tr> <th>i</th> <th>fabs</th> <th>sqrt</th> </tr>";
    out<<"<table> <tr> <th>i</th> <th>fabs</th> <th>sqrt</th> </tr>";
    
    for(int i = 0; i < T; i++) 
    for(int x = 0; x < X; x++) 
            {
                out << "<tr><td>" << (float)u[i][x]<<"</td><td>" << (float)v[i][x]<<"</td></tr>";
            //  out1<<u[100][x]<<endl;
            //  out2<<v[100][x]<<endl;
 
            }
   
            
system("pause");
        return 0;
}
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
07.02.2014, 15:31
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Подскажите правильно ли реализована программа решения уравнения? (C++):

Программа решения квадратного уравнения - C++
Написать программу которая будет решать квадратное уравнение (любое). Помогите пожалуйста написать.

Правильно реализована функция - C++
//функция удалить с массива и уровнение id //function delete from massive and equalize id AddPluss * DelOnePerson (AddPluss * person...

Программа решения набора квадратных уравнения (ax2+bx+c=0) - C++
Программа решения набора квадратных уравнения (ax2+bx+c=0).

Подскажите как правильно подключать библиотеки и решения. (функции уже определены в *.obj) - C++
Вот тут взял самомопальный aes. Там есть 3 исходных кода. 3 - сама программа. 1 и 2 - библиотека и решение.

Программа для решения квадратного уравнения. Где ошибка? - C++
где ошибка при создании программы для вычисления квадратного уравнения?? #include &lt; iostream&gt; #include &lt;conio.h&gt; #include &lt; math.h&gt; ...

Правильно ли реализована передача указателя на контейнер - C++
Я сомневаюсь, правильно ли я передаю указатель на контейнер void put (multimap &lt;int, int&gt; *map2) { multimap &lt;int, int&gt; map3 =...

0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
07.02.2014, 15:31
Привет! Вот еще темы с ответами:

Программа закрывается после решения первого уравнения. Как продолжить ее выполнение? - C++
Здравствуйте форумчани! Нужна ваша помощь... Задача простая, начал изучать С++ не так давно, но хз как реализовать, вообщем суть: ...

Подскажите правильно ли сделана программа? - C++
Вычислить и вывести на экран в виде таблицы значения функции, заданной с помощью ряда Тейлора, на интервале от xнач. до xкон. с шагом dx...

подскажите что делать при вводимых данных 10 25 и 5 20 программа работает не правильно должна выводить 135 подскажите что делать - C++
Вывод В выходной файл выведите одно число - сколько ударов сделали часы за этот отрезок времени.

"Первый после Hello, world!". Программа решения квадратного уравнения - C++
Добрый день, уважаемые форумчане. Покритикуйте пожалуйста, программу вывода на экран в полиномиальной форме и нахождения корней квадратного...


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

Или воспользуйтесь поиском по форуму:
1
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru