Решение систем дифференциальных уравнений методом Рунге-Куты 4 порядка
| Категория реферата: Рефераты по математике
| Теги реферата: конспект, реферат методы
| Добавил(а) на сайт: Juferev.
Предыдущая страница реферата | 1 2 3 4 5 6 7 8 9 10 11 | Следующая страница реферата
Difur;
FOR i:=1 TO n DO BEGIN r2[i]:=dC[i];
C[i]:=cPr[i]+r2[i]*(dX/2);
END;
Difur;
FOR i:=1 TO n DO BEGIN r3[i]:=dC[i];
C[i]:=cPR[I]+r3[i]*dX;
END;
Difur;
FOR i:=1 TO n DO r4[i]:=dC[i];
FOR i:=1 TO n DO rSR[i]:=((r1[i]+r2[i])*(r2[i]+r3[i])*(r3[i]+r4[i]))/6;
END;
PROCEDURE STROKA;
BEGIN
WRITE(f2,'|',x:4:1,'|',c[1]:7:3,'|',c[2]:7:3,'|',c[3]:7:3,'|');
WRITE(f2,sum:3:0,'|',dc[1]:7:3,'|',dc[2]:7:3,'|',dc[3]:7:3,'|');
WRITELN(f2);
END;
PROCEDURE RUN;
BEGIN
WRITE('Step 3: Calculating data and writting results to file : out.rez');
X:=Xn; dX:=0.05;
REPEAT
IF (ABS(x-p)Xk);
WRITELN(' - done.');
END;
PROCEDURE INIT;
BEGIN
ClrScr;
WRITELN('Smith-04: v1.0 (c) 1998 by Mike Smith smith01@home.bar.ru ');
WRITELN;
WRITELN;
WRITE('Step 1: Read data from file : in.dat');
ASSIGN(f1,'in.dat');
RESET(f1);
READLN(f1,C[1],C[2],C[3]);
READLN(f1,k1,k2,k3,k4);
READLN(f1,Xn,Xk,dp,n,eps,p);
WRITELN(' - done.');
ASSIGN(f2,'out.rez');
REWRITE(f2);
WRITE('Step 2: Write header to file : out.rez');
WRITELN(f2,'==========================================================');
WRITELN(f2,'| t,c| Ca,% | Cb,%| Cc,% | SUM | dCa | dCb | dCc |');
WRITELN(f2,'==========================================================');
WRITELN(' - done.');
END;
PROCEDURE DONE;
BEGIN
WRITELN('Step 4: Close all files and exiting...');
CLOSE(f1);
WRITELN(f2,'============================================================');
CLOSE(f2);
WRITELN;
END;
BEGIN
INIT;
RUN;
DONE;
END.
6. Идентификация переменных
Рекомендуем скачать другие рефераты по теме: сочинение описание, отечественная война реферат.
Категории:
Предыдущая страница реферата | 1 2 3 4 5 6 7 8 9 10 11 | Следующая страница реферата