Курсовая работа МХ / TAU1

Program TAU;
Uses graph,crt;
label
m1;
Const
hmax=1;
Var
sigm:array [0..3001] of single;
I1,I2,I3,I4:array [0..3001] of single;


fj,Ij:array[1..4] of real;
h:char;
f,n,i,j:integer;
dt,t:real;
k,b1,a1,a2,a3:real;

Begin
Clrscr;
dt:=0.001;
N:=600;
While N<=2500 do Begin
dt:=0.0001;
While dt<=0.01 do Begin

t:=0;
randomize;
For i:=0 to N do Begin
t:=i*dt;
{ sigm[i+1]:=0.1-0.15*exp(-2*t)+0.05*exp(-5*t)+0.008;}
{ sigm[i+1]:=1-2.25*exp(-t)+1.5*exp(-2*t)-0.25*exp(-3*t);}
sigm[i+1]:=48*exp(-1/0.4*t)-50*exp(-1/0.3*t)-46.666666667*t*exp(-1/0.3*t)+2;
end;

{‘Ґ©з б ўлзЁб«пҐ¬ Ї®¤Ё­вҐаЈ «м­лҐ §­ зҐ­Ёп I1,I2,I3,I4}
t:=0;
for i:=1 to N do Begin
t:=(i-1)*dt;
I1[i]:=1-sigm[i];
I2[i]:=(1-sigm[i])*t;
I3[i]:=(1-sigm[i])*t*t;
I4[i]:=(1-sigm[i])*t*t*t;
end;
{‚лзЁб«Ґ­ЁҐ ­ҐЇ®б।б⢥­­® б ¬Ёе I1,I2,I3,I4}
For j:=1 to 4 do Ij[j]:=0;
For i:=1 to N-1 do Begin
Ij[1]:=Ij[1]+(I1[i]+I1[i+1])/2*dt;
Ij[2]:=Ij[2]+(I2[i]+I2[i+1])/2*dt;
Ij[3]:=Ij[3]+(I3[i]+I3[i+1])/2*dt;
Ij[4]:=Ij[4]+(I4[i]+I4[i+1])/2*dt;
end;
{Ќ е®¤Ё¬ Fj}
Fj[1]:=Ij[1];
Fj[2]:=Ij[1]*Ij[1]-Ij[2];
Fj[3]:=Ij[1]*Ij[1]*Ij[1]-2*Ij[1]*Ij[2]+0.5*Ij[3];
Fj[4]:=-Ij[4]/6+0.5*Ij[1]*Ij[3]+Fj[3]*Ij[1]-Fj[2]*Ij[2];

b1:=-Fj[4]/Fj[3];
a1:=b1+Fj[1];
a2:=Fj[2]+b1*Fj[1];
a3:=Fj[3]+b1*Fj[2];
K:=sigm[N+1];
Writeln ('a1=',a1:6:4,', a2=',a2:6:4,', a3=',a3:6:4,', b1=',b1:6:4,', K=',K:6:4,', n=',n, ', dt=',dt:0:4);
If (a3>0.03) and (a3<0.04) and (a1>0.85) and (a1<1.15) then Begin {goto m1;} Readln; end;
dt:=dt+0.0001;
f:=f+1;
If f mod 23=0 then Begin {h:=readkey; If h='r' then exit;} end;
end;
{h:=Readkey;
If h='r' then exit;}
N:=N+10;
end;
end.