|
DataMuseum.dkPresents historical artifacts from the history of: RC4000/8000/9000 |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about RC4000/8000/9000 Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - download
Length: 1536 (0x600)
Types: TextFile
Names: »cho1txt«
└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ. Detaljer om "HC8000" projekt.
└─⟦0364f57e3⟧
└─⟦913d5b2ff⟧ »bzmodel«
└─⟦667bb35d6⟧ Bits:30007480 RC8000 Dump tape fra HCØ.
└─⟦4334b4c0b⟧
└─⟦913d5b2ff⟧ »bzmodel«
└─⟦this⟧
program for integrating of oregonator with chstiff
22 11 79 17 00 00
begin
real t,dt,t1,t2,h,eps;
real array X(1:3);
procedure fjac(bjac,X,t,F,A);
boolean bjac;
real t;
real array X,F,A;
begin F(1):= 1/0.0002*(X(1)+X(2)-8'-6*X(1)*X(1)-X(1)*X(2));
F(2):= -X(2)+2*0.75*X(3)-X(1)*X(2);
F(3):= 1/300*(X(1)-X(3));
if bjac then
begin A(1,1):= 1/0.0002*(1-2*8'-6*X(1)-X(2));
A(1,2):= 1/0.0002*(1-X(1)); A(1,3):= 0;
A(2,1):= -X(2); A(2,2):= -1-X(1); A(2,3):= 2*0.75;
A(3,1):= 1/300; A(3,2):= 0; A(3,3):= -1/300;
end;
end;
write(out,<:<10>X = :>); setposition(out,0,0); read(in,X);
write(out,<:<10>t1,dt,t2,h,eps = :>); setposition(out,0,0);
read(in,t1,dt,t2,h,eps);
t:= t1;
l1: chstiff(3,X,t,t+dt,fjac,h,eps);
write(out,<:<10>:>,t,h,X(1),X(2),X(3));
setposition(out,0,0);
if t<t2 then goto l1;
end;
end
▶EOF◀