потоковый ввод

@Ленка · Регистрация: 01.08.2009

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

всем привет!
ситуация такая: при выполнении считываются только R2 и h, не могу понять, в чем дело. dev-c++ просто выдает "step t = нажмите любую клавишу для выхода", а если запускать через MS Visual C++ выдает ошибку "damage after normal block at ". в чем дело? помогите, пожалуйста, код-то простой вроде.

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
#include <iostream>       
 
 #include <stdlib.h> 
 
     
using namespace std; 
 
 
int main () {
    double R2=0;
    double h=0;
    double t=0;   
    
    cout<<"R2= ";
    cin>>R2;
        
    cout<<"step r = ";
    cin >>h;
        
    cout<<"step t = ";
    cin >>t;   
          system ("PAUSE")  ;
   
   return 0;
}

@outoftime · 07.04.2010, 16:18

wxDev-C++ все хавает

Добавлено через 13 секунд
ошибок в коде нет

@Ленка · 07.04.2010, 16:20 **[ТС]**

я уже два компилятора скачала( может с машиной что-то не то?

@easybudda · 07.04.2010, 17:00

Ленка, а это так задуманно, что программа три числа считывает и ничего не выводит?
Я синтаксических ошибок тоже не увидел, на что VS ругается - не понятно...

@Ленка · 07.04.2010, 17:03 **[ТС]**

да нет, это кусок программы, я уже вырезала и оформила в отдельную, все равно у меня оба компилятора не хотят слушаться. уже на считывании этом косяк, я даже не могу проверить, что там у меня дальше происходит, хоть об стену бейся.

@easybudda · 07.04.2010, 17:18

Ленка, так Вы бы всю программу выложили...

@Ленка · 07.04.2010, 17:22 **[ТС]**

ну вот что-то в этом духе...

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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
#include <iostream>       
 #include <fstream>
 #include <stdlib.h> 
 
     
using namespace std; 
 
double coef(double x, double y);
double coef1 (double x, int i, int j);
double coef2 (double x, int i, int j);
double fi (double x);
 
int main () {
    double R2=0;
    double h=0;
    double t=0;
    const double T=0.5;
    const double y0=1;
    const double y1=1;
    const double s=1;    
    
    ofstream fout("result.txt");
    
    cout<<"R2= ";
    cin>>R2;
        
    cout<<"step r = ";
    cin >>h;
        
    cout<<"step t = ";
    cin >>t;    
    
    
    int n=R2/h; 
    int m=T/t;  
    double* p = new double [n];
    double* q = new double [n];
    double a, b, c, d, x=0; 
 
    int i=0, j=0;
    
//динамический массив
    double** y = new double* [n];                
    for (int k=0; k<n+1; k++) y[k]=new double[m];
   
 
 
    for (i=0; i<n+1; i++) {                                                   
             y[i][0]=y0;
            }
 
for (j=0; j<n; j++){
    
   p[0]=1;
   q[0]=0; 
    
    for(i=1; i<n+1; i++) {                                                   //прямо ход прогонки
             x=i*h;
             
             a=coef2(x,i,j)/h/h/x;
             b=1/t+coef1(x,i,j)/h/h/x+coef2(x,i,j)/h/x/h;
             c=coef1(x,i,j)/h/h/x;
             d=fi(x)*(-1)-y[i][j]/t;
             
             p[i]=c/(b-a*p[i-1]);                                  
             q[i]=(a*q[i-1]-d)/(b-a*p[i-1]);
             }
             
    y[n][j+1]=(y[n][j]/t+y1*2/h/h/x*coef(R2,y[n][j])*s*h-q[n]*2/h/h/x*coef(R2,y[n][j])-fi(x))/(1/t-2/h/h/x*coef(R2,y[n][j])*(h*s+1+p[n]));
             
    for (i=n; i>-1; i--)  {                                                  //обратный ход прогонки
        y[i][j+1]=y[i+1][j+1]*p[i]+q[i];
        fout << y[i][j+1] << "  ";
        }
   fout << endl;     
        }    
        delete [] y;
        delete [] p;        
        delete [] q;        
      system ("PAUSE")  ;
   
   return 0;
}  
 
 
double coef(double x, double y) {
       return x*y;
       }
       
double coef1 (double x, int i, int j) {
       int m,n;    
       double** y = new double* [n];                
       for (int k=0; k<n+1; k++) y[k]=new double[m];
       return (1/2*(coef(x,y[i+1][j])+coef(x,y[i][j])));
       }
       
double coef2 (double x, int i, int j) {
       int m,n;
       double** y = new double* [n];                
       for (int k=0; k<n+1; k++) y[k]=new double[m];
       return (1/2*(coef(x,y[i][j])+coef(x,y[i-1][j])));
       }
       
double fi (double x) {
       return 1;}

@Nike775 · 07.04.2010, 17:26

у меня никаких ошибок таже не пишет
проверь компилятор

@Ленка · 07.04.2010, 17:28 **[ТС]**

да два компилятора стоит, говорю же. ошибок не выдает, только прога не работает.

@Lаmer · 07.04.2010, 18:46

У меня вылетела с ошибкой "Переменная 'n' используется без инициализации".

А что делает эта программа?

@Ленка · 07.04.2010, 19:39 **[ТС]**

решает нелинейное уравнение теплопроводности, если вам это что-то говорит)

@easybudda · 07.04.2010, 22:42

Сообщение от Ленка

C++
1
2
double* p = new double [n]; 
double* q = new double [n];

А действительно - переменной n даже объявления нигде не увидел, тем более - инициализации...

@Ленка · 07.04.2010, 22:52 **[ТС]**

там выше буквально строкой. dev c++ съел. на ввод пользователем я уже подзабила, вбила руками в константы. теперь еще с памятью что-то((

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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
#include <iostream>       
 #include <fstream>
 #include <stdlib.h> 
 
     
using namespace std; 
 
double coef(double x, double y);
double fi (double x);
 
int main () {
    double R2=1;
    double h=0.05;
    double t=0.01;
    const double T=0.5;
    const double y0=1;
    const double y1=1;
    const double s=1;    
    
    
    
    int n; 
    int m;  
    m=T/t+1;
    n=R2/h+1;
    float* p = new float [n];
    float* q = new float [n];
    float a, b, c, d, x=0; 
    int i=0, j=0;
    
 
    float** y = new float* [n];                
    for (int k=0; k<n+1; k++) y[k]=new float[m];
   
 
   for (i=0; i<n; i++) {                                                   
        for (j=0; j<m; j++)
             y[i][j]=0;
            }
 
    for (i=0; i<n+1; i++) {                                                   
             y[i][0]=y0;
            }           
      x=h;            
                  
 
    ofstream fout("result.txt");     
for (j=0; j<n; j++){
    
   p[0]=1;
   q[0]=0; 
 
    for(i=1; i<n+1; i++) {                                                   
             x=h*i;
             
             a=1/2*(x*y[i][j]+x*y[i-1][j])/h/h/x;
             b=1/t+(1/2*(coef(x,y[i][j])+coef(x,y[i+1][j])))/h/h/x+(1/2*(coef(x,y[i][j])+coef(x,y[i-1][j])))/h/x/h;
             c=(1/2*(coef(x,y[i][j])+coef(x,y[i+1][j])))/h/h/x;
             d=fi(x)*(-1)-y[i][j]/t;
             
             p[i]=c/(b-a*p[i-1]);                                  
             q[i]=(a*q[i-1]-d)/(b-a*p[i-1]);
             }
             
    y[n][j+1]=(y[n][j]/t+y1*2/h/h/x*coef(R2,y[n][j])*s*h-q[n]*2/h/h/x*coef(R2,y[n][j])-fi(x))/(1/t-2/h/h/x*coef(R2,y[n][j])*(h*s+1+p[n]));
             
    for (i=n-1; i>-1; i--)  {                                                  
        y[i][j+1]=y[i+1][j+1]*p[i]+q[i];
        fout << y[i][j+1] << "  ";
        }
   fout << endl;     
        }    
   for (i = 0; i < n; i++) {
    delete [] y[i];
    delete [] y;
     }
 
        delete [] p;        
        delete [] q;        
      system ("PAUSE")  ;
   
   return 0;
}  
 
 
double coef(double x, double y) {
       return y *x;
       }
 
double fi (double x) {
       return 1;}

при отладке переменные a,b,c,d не меняются, постоянно остаются 0, 100, 0, -101.
если посчитать вручную, то там и близко такого нет

@easybudda · 07.04.2010, 23:04

Сообщение от Ленка

for (int k=0; k<n+1; k++) y[k]=new float[m];

условие неверное
for ( int k = 0; k < n; ++k )
И как-то очень много переменных на мой взгляд, да ещё и с такими содержательными именами...

@Ленка · 07.04.2010, 23:43 **[ТС]**

не помогло((

а переменные все нужны, тут уж никуда не деться. все они - либо коэффициенты перед функциями, либо счетчики, как их называть по-другому не знаю, у меня фантазии не хватает.

@easybudda · 08.04.2010, 00:10

Сообщение от Ленка

не помогло

там в трёх местах такое. точно всё поправили?

Добавлено через 7 минут
хотя там много всего... например:

Сообщение от Ленка

a=1/2*(x*y[i][j]+x*y[i-1][j])/h/h/x;

a, а так же b, c, d... нигде не объявлены. выражение 1/2 даст ноль (целочисленное деление), да и с индексами всё мутно...

@Ленка · 08.04.2010, 01:21 **[ТС]**

тут я уже заменила 1 на 1., с этим теперь все норм, ошибка с память возникает в 17 строке,

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
for (j=0; j<m-1; j++){
    
   p[0]=1;
   q[0]=0; 
    for(i=1; i<n-1; i++) {                                                   
             x=h*i;
             
             a=1./2*(x*y[i][j]+x*y[i-1][j])/h/h/x;
             b=1./t+(1./2*(coef(x,y[i][j])+coef(x,y[i+1][j])))/h/h/x+(1./2*(coef(x,y[i][j])+coef(x,y[i-1][j])))/h/x/h;
             c=(1./2*(coef(x,y[i][j])+coef(x,y[i+1][j])))/h/h/x;
             d=fi(x)*(-1.)-y[i][j]/t;
             
             p[i]=c/(b-a*p[i-1]);                                  
             q[i]=(a*q[i-1]-d)/(b-a*p[i-1]);
             }
             
    y[n][j+1]=(y[n][j]/t+y1*2./h/h/x*coef(R2,y[n][j])*s*h-q[n]*2./h/h/x*coef(R2,y[n][j])-fi(x))/(1./t-2./h/h/x*coef(R2,y[n][j])*(h*s+1+p[n]));
             
    for (i=n-1; i>-1; i--)  {                                                 
        y[i][j+1]=y[i+1][j+1]*p[i]+q[i];
        fout << y[i][j+1] << "  ";
        }
   fout << endl;     
        }

Добавлено через 1 час 8 минут
всем спасибо, все заработало) с алгоритмом косяк вышел небольшой)

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .
И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»	SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .

@outoftime ║XLR8║ 1212 / 909 / 270 Регистрация: 25.07.2009 Сообщений: 4,360 Записей в блоге: 5
	07.04.2010, 16:18
	wxDev-C++ все хавает Добавлено через 13 секунд ошибок в коде нет 1

@Ленка 0 / 0 / 0 Регистрация: 01.08.2009 Сообщений: 24
	07.04.2010, 16:20 [ТС]
	я уже два компилятора скачала( может с машиной что-то не то? 0

@easybudda Модератор 12843 / 7592 / 1766 Регистрация: 25.07.2009 Сообщений: 13,973
	07.04.2010, 17:00
	Ленка, а это так задуманно, что программа три числа считывает и ничего не выводит? Я синтаксических ошибок тоже не увидел, на что VS ругается - не понятно... 0

@Ленка 0 / 0 / 0 Регистрация: 01.08.2009 Сообщений: 24
	07.04.2010, 17:03 [ТС]
	да нет, это кусок программы, я уже вырезала и оформила в отдельную, все равно у меня оба компилятора не хотят слушаться. уже на считывании этом косяк, я даже не могу проверить, что там у меня дальше происходит, хоть об стену бейся. 0

@easybudda Модератор 12843 / 7592 / 1766 Регистрация: 25.07.2009 Сообщений: 13,973
	07.04.2010, 17:18
	Ленка, так Вы бы всю программу выложили... 0

@Nike775 3 / 2 / 0 Регистрация: 01.04.2010 Сообщений: 30
	07.04.2010, 17:26
	у меня никаких ошибок таже не пишет проверь компилятор 0

@Ленка 0 / 0 / 0 Регистрация: 01.08.2009 Сообщений: 24
	07.04.2010, 17:28 [ТС]
	да два компилятора стоит, говорю же. ошибок не выдает, только прога не работает. 0

@Lаmer 7 / 7 / 7 Регистрация: 24.03.2010 Сообщений: 25
	07.04.2010, 18:46
	У меня вылетела с ошибкой "Переменная 'n' используется без инициализации". А что делает эта программа? 0

@Ленка 0 / 0 / 0 Регистрация: 01.08.2009 Сообщений: 24
	07.04.2010, 19:39 [ТС]
	решает нелинейное уравнение теплопроводности, если вам это что-то говорит) 0

@Ленка 0 / 0 / 0 Регистрация: 01.08.2009 Сообщений: 24
	07.04.2010, 23:43 [ТС]
	не помогло(( а переменные все нужны, тут уж никуда не деться. все они - либо коэффициенты перед функциями, либо счетчики, как их называть по-другому не знаю, у меня фантазии не хватает. 0