Программирование -компьютерное моделирование

Для решения дифференциального уравнения выберем наиболее распространенный язык программирования Паскаль (Турбо Паскаль).

Алгоритм программы:

program biocheml; uses

crt, GRAPH; Label metl,met2; VAR

I: integer;

ch : char;

s : string;

C1,C2,CATP,CATP_MAX,CCRP,CCRP_MAX:REAL;

V1,V2:REAL;

VC Al, VCA2,TEND: REAL;

ml,m2,m3,mash,t,dt:real;

procedure InpGraph; var

Driver, Mode, Code : Integer; begin

Driver: = Detect;

InitGraph(Driver, Mode, 'c:\tp\bgi'); end;

function DigStr(S: Longint; A:Byte): string;

var SI: string[10];

begin

Str(S:A, SI);

DigStr:= SI; end;{function DigStr}

function DigStrR(S: Real; A,B:Byte): string;

var SI: stringflO];

begin

Str(S:A:B, SI);

DigStrR :« SI; end; {function DigStrR}

BEGIN ClrScr; InpGraph;

C1:=3;C2: = 15;C3:=1; {THE SPEED EXPENDITER ATP ,

CRP}

CATP:=5;CATP_MAX:=5;CCRP:=15;CCRP_MAX:=15;DT:=0.001

T:=0;

gotoXY(6,3);writeln('Input ');

gotoXY(3,5);writeln('ATP - ");

gotoXY(3,7);writeln('CrP - ");

gotoXY(3,9);writeln('Scale - ');

gotoXY(3,ll);writeln('time activity - ');

gotoXY(3,13);writeln( Intensity - ');

gotoXY(3,15);writeln('time investignation - ');

gotoXY(49,5); readln(ml); gotoXY(49,7); readln(m2); gotoXY(49,9); readln(mash); gotoXY(49,ll); readln(ta); gotoXY(49,13); readln(i); gotoXY(49,15); readln(tend);

cl:=cl*ml;c2:=c2*m2;

clrscr;

InpGraph;

line(100,300,600,300);

line(100,50,100,300); textcolor(2); setcolor(l); setbkcolor(O); MET1:

Г I****************************'1

T:=0;

WHILE T<TEND DO BEGIN

V1: =C 1 *i*(C ATP/C ATP_MAX);

V2:=C2*(1-CATP/CATP_MAX)*CCRP/CCRP_MAX; CATP:=CATP+(V2-V1)*DT; CCRP: =CCRP+(V3-V2)*DT;

OUTTEXTXY(3,360,'Time =' +DigStrR(T,9,3));
OUTTEXTXY(3,370,'Power =*+DigStrR(vl,19,3));
OUTTEXTXY(3,380,'ATP ='+DigStrR(CATP,19,3));
OUTTEXTXY(3,390,'CrP ='+DigStrR(CCRP,19,3));

putpixel(round(100+t*500*mash), Round(300-vl*100),7); putpixel(round(100+t*500*mash), Round(300-catp*100),4); putpixel(round(100+t*500*mash), Round(300-ccrp*20),3); T:=T+DT;

if t>=tend then I:=0;

ch:=readkey;

BAR(3,350,200,410); END; {WHILE}

GOTO MET1; CloseGraph; END.