|
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: 768 (0x300)
Types: TextFile
Names: »chtesttxt«
└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ. Detaljer om "HC8000" projekt.
└─⟦0364f57e3⟧
└─⟦this⟧ »chtesttxt«
└─⟦667bb35d6⟧ Bits:30007480 RC8000 Dump tape fra HCØ.
└─⟦4334b4c0b⟧
└─⟦this⟧ »chtesttxt«
program for testing chstiff
21 11 79 17 00 00
begin
real t,dt,t1,t2,h,eps;
real array X(1:2);
procedure fjac(bjac,X,t,F,A);
boolean bjac;
real t;
real array X,F,A;
begin F(1):= X(2);
F(2):= -X(1);
if bjac then
begin A(1,1):= 0; A(1,2):= 1;
A(2,1):= -1; A(2,2):= 0;
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(2,X,t,t+dt,fjac,h,eps);
write(out,<:<10>:>,t,h,X(1),X(2),sqrt(X(1)**2+X(2)**2));
if t<t2 then goto l1;
end;
end
▶EOF◀