Разобраться с функцией Omega

@ZverCi · Регистрация: 07.12.2011

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

Добрый вечер!
У меня такая проблема:есть функция function Omega(w: real; t0: real; t1: real): real; которая должна работать следующим образом. Есть глобальное время в программе,которая считает значение на каждом шаге. Мне нужно это текущее время в программе сравнивать с t1 и t0.

То есть пока t0<=t<=t1 то значение w = 0,когда же значение t будет находиться вне этого отрезка, то должно равняться значению введенному в поле ввода.
Помоги разобраться.

Delphi
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
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Math;
 
type
  TForm1 = class(TForm)
    Button1: TButton;
    Edit1: TEdit;
    Edit2: TEdit;
    Edit3: TEdit;
    Label1: TLabel;
    Label2: TLabel;
    Label3: TLabel;
    Edit4: TEdit;
    Label4: TLabel;
    Edit5: TEdit;
    Label5: TLabel;
    Edit6: TEdit;
    Label6: TLabel;
    Edit7: TEdit;
    Label7: TLabel;
    Button2: TButton;
    Edit8: TEdit;
    Label8: TLabel;
    Edit9: TEdit;
    Label9: TLabel;
    CheckBox1: TCheckBox;
    Edit10: TEdit;
    Label10: TLabel;
    Edit11: TEdit;
    Label11: TLabel;
    Edit13: TEdit;
    Label13: TLabel;
    Button3: TButton;
    procedure Button1Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
    procedure Button3Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
end;
 
type
 
rvector = record
    x : double;
    z : double;
    Vx: double;
    Vz: double;
end;
 
vector = class (TObject)
  private
    x : double;
    z : double;
    Vx: double;
    Vz: double;
 
    protected
    // Âíåøíå äîñòóïíûå ïîëÿ è ìåòîäû
    public
      constructor Create(x_ : double; z_ : double; Vx_ : double; Vz_: double); overload;
      //procedure mult(a: double);
    // Âíåøíå äîñòóïíûå è âèäèìûå ïîëÿ è ìåòîäû
    published
      // Îáðàòèòå âíèìàíèå, ÷òî ñâîéñòâà äîëæíû èñïîëüçîâàòü ðàçíûå èìåíà äëÿ ìåñòíûõ îïðåäåëåíèé
      property x_  : double
        read   x;
      property z_  : double
        read   z;
      property Vx_ : double
        read   Vx;
      property Vz_ : double
        read   Vz;
end;                    // Êîíåö îïðåäåëåíèÿ êëàññà vector
 
 
var
  Form1: TForm1;
  Vx0, Vz0, z0: Double;
  D, ro, Vv, T_, p: Double;
  Ow,Ot0,Ot1: double;
 
 
 
implementation
 
{$R *.dfm}
constructor vector.Create(x_, z_, Vx_, Vz_ : double);
begin
  self.x:=x_;
  self.z:=z_;
  self.Vx:=Vx_;
  self.Vz:=Vz_;
end;
 
function mult(U__: rvector; a: double) : rvector;
begin
  //mult:=vector.Create(0,0,0,0);
  mult.x:=U__.x*a;
  mult.z:=U__.z*a;
  mult.Vx:=U__.Vx*a;
  mult.Vz:=U__.Vz*a;
end;
 
function sum(U1: rvector; U2: rvector) : rvector;
begin
  //sum:=vector.Create(0,0,0,0);
  sum.x:=U1.x+U2.x;
  sum.z:=U1.z+U2.z;
  sum.Vx:=U1.Vx+U2.Vx;
  sum.Vz:=U1.Vz+U2.Vz;
end;
 
function razn(U1: rvector; U2: rvector) : rvector;
begin
  //razn:=vector.Create(0,0,0,0);
  razn.x:=U1.x-U2.x;
  razn.z:=U1.z-U2.z;
  razn.Vx:=U1.Vx-U2.Vx;
  razn.Vz:=U1.Vz-U2.Vz;
end;
 
function WV(z: double): double;
begin
  if Form1.CheckBox1.Checked=true then
    begin
    if z>=0.001 then
      begin
      WV:=(Vv/0.4)*ln(z/0.0003);
      end
    else
      begin
      WV:=0
      end;
    end
  else
    begin
    WV:=Vv;
    end;
 
end;
 
 
 
function myu(T: double): double;
begin
  myu:=1.717*0.00001*exp(0.683*ln(T/273.0));
end;
 
function nyu(T: double; p_: double): double;
begin
  nyu:=myu(T)/(p_/(287.4*T));
end;
 
function Omega(w: real; t0: real; t1: real): real;
begin
   if(t>=t0) and (t<=t1)  then
    w:=0;
     end;
   else
      begin
        w:=Ow;
      end;
 
 
function Re(Vx: double; Vz: double; Ux: double; Uz: double; D: double; nyu_: double): double;
begin
           Re:= D*sqrt((Vx-Ux)*(Vx-Ux)+(Vz-Uz)*(Vz-Uz))/nyu_;
end;
 
 
function modVp(Vx: double; Vz: double; Ux: double; Uz: double): double;
begin
 
      modVp:= sqrt((Vx-Ux)*(Vx-Ux)+(Vz-Uz)*(Vz-Uz));
 
end;
 
 
function PSI(Re: double): double;
var
      Psi1, Psi2: double;
begin
  if    (Re <= 0.2) then
  begin
      Psi1:= 0;
      Psi2:= 0;
      end
      else if (Re > 0.2) and (Re < 10) then
      begin
        Psi1:= 0.15*exp( 0.667*ln(Re) );
        Psi2:= -0.023*ln(Re/5)-0.0017*( ln(Re/5)*ln(Re/5) );
          end
            else if (Re >= 10) then
            begin
              Psi1:= 0.15*exp( 0.667*ln(Re) );
              //Psi2:= -0.023*ln(Re/5)-0.0017*( ln(Re/5)*ln(Re/5) );
             Psi2:= -0.004*exp( 1.315*ln(Re/5) );
            end;
        PSI:= (1+Psi1+Psi2);
end;
 
 
function F(U_: rvector; t: double): rvector;
var
  g: double;
begin
  g:=9.81;
  //F:=vector.Create(0,0,0,0);
  F.x:=U_.Vx;
  F.z:=U_.Vz;
  //F.Vx:=0;
  //F.Vz:=-g;
   F.Vx:= -18*(1.27/(ro*1000))*(0.15*0.0001)*(1/(D*D*(0.000001*0.000001)))*(( (U_.Vx-WV(U_.z)) - 0.85*( (U_.Vz-WV(U_.z))*((D*Omega(Ow,Ot0,Ot1))/modVp(U_.Vx,U_.Vz,WV(U_.z),0)))*PSI(Re(U_.Vx, U_.Vz, WV(U_.z),0,D*0.000001,nyu(T_,p)))));
  F.Vz:= -g -18*(1.27/(ro*1000))*(0.15*0.0001)*(1/(D*D*(0.000001*0.000001)))*(( (U_.Vz-WV(U_.z)) + 0.85* ( (U_.Vx-WV(U_.z))*((D*Omega(Ow,Ot0,Ot1))/modVp(U_.Vx,U_.Vz,WV(U_.z),0)))*PSI(Re(U_.Vx, U_.Vz, WV(U_.z),0,D*0.000001,nyu(T_,p)))));
end;
 
 
function Ff(U_: rvector; t: double): rvector;
var
  g: extended;
  ksi: double;
  Vx, Vz: double;
begin
  g:=9.81;
  ksi:=(p/(287.4*T_))/(ro*1000);  //îòíîøåíèå ïëîòíîñòåé âîçäóõà è âåùåñòâà ïåñ÷èíêè
 
  //F:=vector.Create(0,0,0,0);
  Ff.x:=U_.Vx;
  Ff.z:=U_.Vz;
  //F.Vx:=0;
  //F.Vz:=-g;
  Vx:= -18*ksi*nyu(T_, p)*(1/(D*D*(0.000001*0.000001)))*(( (U_.Vx-WV(U_.z)) -0.85*( (U_.Vz-WV(U_.z)) *((D*Omega(Ow,Ot0,Ot1))/modVp(U_.Vx,U_.Vz,WV(U_.z),0)) )*PSI(Re(U_.Vx, U_.Vz, WV(U_.z),0,D*0.000001,nyu(T_,p)))));
  Ff.Vx:= Vx - (27/8)*ksi*(1/(D*0.000001))*(U_.Vx-WV(U_.z))*sqrt((U_.Vx-WV(U_.z))*(U_.Vx-WV(U_.z))+(U_.Vz-0)*(U_.Vz-0));
 
  Vz:= -g -18*ksi*nyu(T_, p)*(1/(D*D*(0.000001*0.000001)))*( ((U_.Vz-WV(U_.z)) +0.85*( (U_.Vx-WV(U_.z))*((D*Omega(Ow,Ot0,Ot1))/modVp(U_.Vx,U_.Vz,WV(U_.z),0)) )*PSI(Re(U_.Vx, U_.Vz, WV(U_.z),0,D*0.000001,nyu(T_,p)))));
  Ff.Vz:= Vz - (27/8)*ksi*(1/(D*0.000001))*(U_.Vz-0)*sqrt((U_.Vx-WV(U_.z))*(U_.Vx-WV(U_.z))+(U_.Vz-0)*(U_.Vz-0));
end;
 
 
procedure TForm1.Button1Click(Sender: TObject);
var
    i : WORD ;
    V0, Vx, Vz : Double;
  x, z: Double;
  Alfa: Double;
  t, dt: Double;
  V: Double;
  U: rvector;
  Uu: rvector;
  k, k1, k2, k3, k4: rvector;
  Outfile : TextFile;
  Re: Double;
begin
t:=0;
 
Vx0:= StrToFloat(Edit1.Text);
Vz0:= StrToFloat(Edit2.Text);
z0:= StrToFloat(Edit5.Text);
D:=StrToFloat(Edit6.Text);
ro:=StrToFloat(Edit7.Text);
dt:=StrToFloat(Edit3.Text);
Vv:=StrToFloat(Edit4.Text);
Ow:= StrToFloat(Edit13.Text);
Ot0:= StrToFloat(Edit10.Text)*0.001;
Ot1:= StrToFloat(Edit11.Text)*0.001;
 
 
 
SetCurrentDir('D:\Temp');
AssignFile(Outfile, 'temp.txt');
Rewrite(Outfile);
 
//Vx:= V0*cos(Alfa*3.1415/180.0);
//Vz:= V0*sin(Alfa*3.1415/180.0);
//U:= vector.Create(0, 0, Vx, Vz);
U.x:=0;
U.z:=z0;
U.Vx:=Vx0;
U.Vz:=Vz0;
 
 
Writeln(Outfile, format('%.8f'+#9+'%.8f'+#9+'%.8f'+#9+'%.8f'+ #9+'%.8f'+ #9+'%.8f'+ #9+'%.8f' + #9 +'%.4f', [t*1000, U.x*1000, U.z*1000, U.Vx, U.Vz, Uu.Vx, Uu.Vz , WV(U.z)]));
 
//k:=vector.Create(0,0,0,0);
//k1:=vector.Create(0,0,0,0);
//k2:=vector.Create(0,0,0,0);
//k3:=vector.Create(0,0,0,0);
//k4:=vector.Create(0,0,0,0);
 
while U.z>=0. do
begin
  k1:=mult(F(U,t),dt);
  k2:=mult(F(sum(U, mult(k1, 0.5)), t+0.5*dt), dt);
  k3:=mult(F(sum(U, mult(k2, 0.5)), t+0.5*dt), dt);
  k4:=mult(F(sum(U, k3), t+0.5*dt), dt);
  k:=sum(k1, mult(k2, 2));
  k:=sum(k, mult(k3, 2));
  k:=sum(k, k4);
 
  U:=sum(U, mult(k, 1.0/6.0));
  Uu:=F(U,t);
  t:= t+dt;
  if U.z>=0. then Writeln(Outfile, format('%.8f'+#9+'%.8f'+#9+'%.8f'+#9+'%.8f'+ #9+'%.8f'+ #9+'%.8f'+ #9+'%.8f' + #9 +'%.4f', [t*1000, U.x*1000, U.z*1000, U.Vx, Uu.Vx, Uu.Vz, U.Vz, WV(U.z)]));
 
end;
 
//Append(Outfile);
 
//Writeln(Outfile, '111');
 
//U.Free;
//k.Free;
//k1.Free;
//k2.Free;
//k3.Free;
//k4.Free;
 
 
Flush(Outfile);
CloseFile(Outfile);
end;
 
 
 
procedure TForm1.Button2Click(Sender: TObject);
var
    i : WORD ;
    V0, Vx, Vz : Double;
  x, z: Double;
  Alfa: Double;
  t, dt: Double;
  V: Double;
  U: rvector;
  Uu: rvector;
  k, k1, k2, k3, k4: rvector;
  Outfile : TextFile;
  Re: Double;
begin
 
t:=0;
 
Vx0:= StrToFloat(Edit1.Text);
Vz0:= StrToFloat(Edit2.Text);
z0:= StrToFloat(Edit5.Text);
D:=StrToFloat(Edit6.Text);
ro:=StrToFloat(Edit7.Text);
dt:=StrToFloat(Edit3.Text);
Vv:=StrToFloat(Edit4.Text);
T_:=273+StrToFloat(Edit8.Text);
p:=StrToFloat(Edit9.Text);
Ow:= StrToFloat(Edit13.Text);
Ot0:= StrToFloat(Edit10.Text);
Ot1:= StrToFloat(Edit11.Text);
 
 
SetCurrentDir('D:\Temp');
AssignFile(Outfile, 'temp1.txt');
Rewrite(Outfile);
 
//Vx:= V0*cos(Alfa*3.1415/180.0);
//Vz:= V0*sin(Alfa*3.1415/180.0);
//U:= vector.Create(0, 0, Vx, Vz);
U.x:=0;
U.z:=z0;
U.Vx:=Vx0;
U.Vz:=Vz0;
Uu:=Ff(U,t);
 
Writeln(Outfile, format('%.8f'+#9+'%.8f'+#9+'%.8f'+#9+'%.8f'+#9+'%.8f'+#9+'%.8f'+ #9+'%.8f'+ #9+'%.8f'+ #9+'%.8f' + #9 +'%.4f', [t*1000, U.x*1000, U.z*1000, U.Vx, U.Vz, Uu.Vx, Uu.Vz, WV(U.z), sqrt((U.Vx*U.Vx)+(U.Vz*U.Vz)), sqrt ((Uu.Vx*Uu.Vx)+(Uu.Vz*Uu.Vz)) ]));
 
while U.z>=0. do
begin
  k1:=mult(Ff(U,t),dt);
  k2:=mult(Ff(sum(U, mult(k1, 0.5)), t+0.5*dt), dt);
  k3:=mult(Ff(sum(U, mult(k2, 0.5)), t+0.5*dt), dt);
  k4:=mult(Ff(sum(U, k3), t+0.5*dt), dt);
  k:=sum(k1, mult(k2, 2));
  k:=sum(k, mult(k3, 2));
  k:=sum(k, k4);
 
  U:=sum(U, mult(k, 1.0/6.0));
  Uu:=Ff(U,t);
  t:= t+dt;
  if U.z>=0. then Writeln(Outfile, format('%.8f'+#9+'%.8f'+#9+'%.8f'+#9+'%.8f'+#9+'%.8f'+#9+'%.8f'+ #9+'%.8f'+ #9+'%.8f'+ #9+'%.8f' + #9 +'%.4f', [t*1000, U.x*1000, U.z*1000, U.Vx, U.Vz, Uu.Vx, Uu.Vz, WV(U.z), sqrt((U.Vx*U.Vx)+(U.Vz*U.Vz)), sqrt ((Uu.Vx*Uu.Vx)+(Uu.Vz*Uu.Vz)) ]));
 
end;
 
Flush(Outfile);
CloseFile(Outfile);
end;
 
procedure TForm1.Button3Click(Sender: TObject);
var
      CD,Re_: double;
      Outfile : TextFile;
begin
 
SetCurrentDir('D:\Temp');
AssignFile(Outfile, 'temp2.txt');
Rewrite(Outfile);
 
        Re_:=0.001;
        while Re_ < 500 do
            begin
              CD:= 24/Re_ * PSI(Re_);
              Re_:= Re_ + 0.001;
              Writeln(Outfile, format('%.8f'+#9+'%.8f', [Re_,CD] ));
            end;
 
Flush(Outfile);
CloseFile(Outfile);
 
end;
 
end.

@ZverCi · 03.04.2012, 16:11 **[ТС]**

Да,температуру. Если через дебаг посмотреть значение,которому равна T то это значение 303.

@Mawrat · 03.04.2012, 16:30

Да, точно. Если через отладчик проследить, оказывается, в функции nyu(T: double; p_: double) параметр p_ оказывается равным нулю. А так как в функции p_ тоже оказывается в знаменателе - то происходит деление на ноль. Теперь надо проследить, как нулевое значение попадает в функцию nyu() через второй параметр.

@ZverCi · 03.04.2012, 16:34 **[ТС]**

p это давление и оно задается так же в поле ввода и равно 101325. Не понимаю, как оно может равняться 0

Добавлено через 1 минуту
я сейчас добавил такую строчку

Delphi
1
p:=StrToFloat(Edit9.Text);

@Mawrat · 03.04.2012, 16:36

Ага - вот теперь без ошибок выполняется.

@ZverCi · 03.04.2012, 16:39 **[ТС]**

Но программа не работает так как мне надо!
если вы введете в поля ввода w = 1000, t начальное = 30 и t конечное = 80 то увидете,что она выдает лишь одну строчку и цикл срывается,это прослеживается черед отладчик...а если эти значения равны 0,то тогда опять же через отладчик видно,как она правильно работает...я уже не понимаю в чем дело

@Mawrat · 03.04.2012, 16:40

ZverCi, ещё, когда отлаживаешь проект, отключай оптимизацию. У тебя сейчас оптимизация включена в проекте. Из-за оптимизации в программе некоторые действия могут выполняться неожиданным способом. И на некоторых строках может отсутствовать возможность поставить точку останова.
Как отключить оптимизацию.

@ZverCi · 03.04.2012, 16:44 **[ТС]**

новая версия программы

@Mawrat · 03.04.2012, 16:58

ZverCi, я к сожалению тут вряд-ли что-то посоветую... - Это уже вопросы по конкретной мат. модели. Там, в самом деле, при тех параметрах w = 1000, t = 30..80, уже на первой итерации проекция вектора U на ось Z оказывается отрицательной. Если так не должно быть - это мат. модель надо копать или алгоритм.

@ZverCi · 03.04.2012, 17:00 **[ТС]**

спасибо!за помощь

@Mawrat · 03.04.2012, 17:02

Кстати, на форуме есть научный раздел. В нём есть математический подраздел: Математика. И там есть подразделы: Дифференциальные уравнения и Численные методы. - Там можно попробовать обсудить математическую и расчётную часть задачи.

@ZverCi · 03.04.2012, 17:56 **[ТС]**

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

@ZverCi · 03.04.2012, 17:59 **[ТС]**

вот эти формулы из курсовой,которые использовались при написании старой программы

Новые блоги и статьи Все статьи Все блоги /
Был там один разговор по поводу свободы в материальном мире. kumehtar 19.08.2026 Суть: рассматривается живое существо, оказавшееся внутри довольно странной системы (этого мира) и пытающееся обустроить в ней свой кусок пространства. Жизнь действительно предъявляет каждому. . .	Когда логика программы не спасает от человеческих ошибок Maks 18.08.2026 В последнее время всё чаще и чаще сталкиваюсь с таким явлением, как абсолютная невнимательность (или глупость) пользователей. Проявляется это чаще всего на работе в коллективе. Допустим, человек с. . .	Лето уходит kumehtar 17.08.2026	Мысли в слух kumehtar 17.08.2026 Забавно, насколько сейчас стала доступна информация. Например о магии, духовном развитии, медитациях, и других подобных направлениях, ранее зачастую тайных, передаваемых от учителя к ученику. Хотя. . .
Перемещение строк из ТЧ в другой документ с учетом текущего пробега Maks 17.08.2026 Реализация из решения ниже выполнена на примере нетипового документа "Автозапчасти", с ТЧ "Шины". За основу взят алгоритм отсюда: https:/ / www. cyberforum. ru/ blogs/ 359708/ 10838. html Задача: . . .	Саморегулирующийся социальный контракт для сервера cross-section. Hrethgir 14.08.2026 С кодом конечно таких глубоких размышлений пока не было, впрочем я уже привык к алгоритмизации. Суть предмета записи: снова в диалоге с нейросетью (я взял пока себе ник для учётки админа - Rector). . . .	Часы электронные Uhbif79 12.08.2026 Выкладываю программу часов. Программа позволяет: 1. Использовать системное время и дату, 2. Есть возможность вводить время и дату вручную. 3. Реализованы 2 будильника: начало и конец рабочего дня. . . .	Часы с будильником на основе класса QLCDNumber Uhbif79 12.08.2026 Всем добрый день, выкладываю программу часов с будильником на основе класса QLCDNumber. Здесь я пробовал самостоятельно создавал классы, впервые столкнулся с видимостью переменной одного класса из. . .

@ZverCi 0 / 0 / 0 Регистрация: 07.12.2011 Сообщений: 32
	03.04.2012, 16:11 [ТС]
	Да,температуру. Если через дебаг посмотреть значение,которому равна T то это значение 303. 0

@Mawrat 13117 / 5898 / 1708 Регистрация: 19.09.2009 Сообщений: 8,809
	03.04.2012, 16:30
	Да, точно. Если через отладчик проследить, оказывается, в функции nyu(T: double; p_: double) параметр p_ оказывается равным нулю. А так как в функции p_ тоже оказывается в знаменателе - то происходит деление на ноль. Теперь надо проследить, как нулевое значение попадает в функцию nyu() через второй параметр. 1

@Mawrat 13117 / 5898 / 1708 Регистрация: 19.09.2009 Сообщений: 8,809
	03.04.2012, 16:36
	Ага - вот теперь без ошибок выполняется. 1

@ZverCi 0 / 0 / 0 Регистрация: 07.12.2011 Сообщений: 32
	03.04.2012, 16:39 [ТС]
	Но программа не работает так как мне надо! если вы введете в поля ввода w = 1000, t начальное = 30 и t конечное = 80 то увидете,что она выдает лишь одну строчку и цикл срывается,это прослеживается черед отладчик...а если эти значения равны 0,то тогда опять же через отладчик видно,как она правильно работает...я уже не понимаю в чем дело 0

@Mawrat 13117 / 5898 / 1708 Регистрация: 19.09.2009 Сообщений: 8,809
	03.04.2012, 16:40
	ZverCi, ещё, когда отлаживаешь проект, отключай оптимизацию. У тебя сейчас оптимизация включена в проекте. Из-за оптимизации в программе некоторые действия могут выполняться неожиданным способом. И на некоторых строках может отсутствовать возможность поставить точку останова. Как отключить оптимизацию. 1

@Mawrat 13117 / 5898 / 1708 Регистрация: 19.09.2009 Сообщений: 8,809
	03.04.2012, 16:58
	ZverCi, я к сожалению тут вряд-ли что-то посоветую... - Это уже вопросы по конкретной мат. модели. Там, в самом деле, при тех параметрах w = 1000, t = 30..80, уже на первой итерации проекция вектора U на ось Z оказывается отрицательной. Если так не должно быть - это мат. модель надо копать или алгоритм. 1

@ZverCi 0 / 0 / 0 Регистрация: 07.12.2011 Сообщений: 32
	03.04.2012, 17:00 [ТС]
	спасибо!за помощь 0

@Mawrat 13117 / 5898 / 1708 Регистрация: 19.09.2009 Сообщений: 8,809
	03.04.2012, 17:02
	Кстати, на форуме есть научный раздел. В нём есть математический подраздел: Математика. И там есть подразделы: Дифференциальные уравнения и Численные методы. - Там можно попробовать обсудить математическую и расчётную часть задачи. 1