DataMuseum.dk

top - metrics - download

    Length: 1536 (0x600)
    Types: TextFile
    Names: »øv5f1«

└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
    └─⟦b2ec5d50f⟧ 
        └─⟦this⟧ »øv5f1« 

TextFile

program for fitting of solubility against ionic strength
begin
  integer n,i,ip,in,cp,cn,d0;
  real q,fa,fb;
READC:
  cp:=readi(<:charge of cation:>);
  cn:=readi(<:charge of anion:>);
  if cn<0 then cn:=-cn;
  if cp<=0 or cn-0 then
  begin
    write(out,<:illegal charge:>);
    goto READC
  end;
  ip:=d0//cp;
  d0:=0;
  repeat
     d0:=d0+cp
  until (do mod cn)=0;
  in:=d0//cn;
  n:=read(<:number of datapoints:>);
  begin
    integer r;
    real array X,Y(1:n),C,S(0:r);
    r:=readi(<:degree of approx. pol.:>);
    write(out,<:<'nl'>added ionic strength, conc. of cation:>);
    outendcur('nl'); 
    for i:= 1 step 1 until n do read(in,X(i),Y(i));
    for i:=1 step 1 until n do
    begin
      X(i):=sqrt(X(i)+.5*(cp*cp*ip+cn*cn*in)*Y(i)/ip);
      Y(i):=ln(Y(i)/ip)
    end;
    q:=polfit1(n,X,Y,C,S,r);
    if q>0 then q:=sqrt(q);
    fb:=(ip+in)/ln10;
    fa:=fb(ip*cp*cp+in*cn*cn);
    write(out,<:<'nl'> log L =  :>,
         <<-ddd.dd>,C(0*fb+ln(ip**ip+in**in)/ln10,
         <:  +-:>,<<ddd.dd>,S(0)*fb));
    for i:=1 step1 until m do
      write(out,
           <:<'nl'> c(:>,<<d>,i,<:) =  :>,<<-dd.ddd>,C(i)*fa,
           <:  +-:>,<<dd.ddd>,S(i)*fa);
    write(out,<:<'nl'> s.t.d on log s :>.<<dd.ddd>,q/ln10);
  end
end
1
1
5
2
.00490, .001930
.01257, .002025
.02565, .002158
.04081, .002266
.05422, .002359
▶EOF◀