DataMuseum.dk

top - metrics - download

    Length: 2304 (0x900)
    Types: TextFile
    Names: »gofuncfit«

└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ.  Detaljer om "HC8000" projekt.
    └─⟦0364f57e3⟧ 
        └─⟦1248b0c55⟧ »gobib« 
            └─⟦this⟧ 

TextFile

;gosav
lookup svdecomp
if ok.no
(i extsvdecomp)
ftxt=set 20
ftxt=edit algfuncfit
l./integer N/,r/nt, nf/n1, nt/,r/i,/i, ii,/,
l./in,N,p/,g/nf/n1/,
l./lo(1:nt)/,r/lo/lo,sgm/,
l./S*A(j,j)/,r/A(j,j)/dy(j)/,i!
    dy(j):= sum((sgm(ii)*A(ii,j)/A(j,n+1))**2,ii,n1,n);
!,l./C(1:n/,l1,d1,i!
   write(res,nl,3,<:
Uncorrelated parameters, X = V*w*pm, with uncertainties::>,nl,2);
   ii:= 0;
   for i:=1 step 1 until n do begin
     y:= sqrt(Sold*sgm(i));
     h:= s:= sum(A(i,j)*A(j,n+1)*pm(no(j)),j,1,n);
     if i<n1 then fx:= real <<      -d.dd'-d> else begin
       j:= entier(ln(abs h)/ln10);
       k:= entier(ln(y/n)/ln10);
       j:= j-k+1; l:= 14;
       if j<1 then begin j:= 2; s:= y end;
       if j>ii then ii:= j;
       fx:= format(l,j,abs s,fy);
     end;   exp(i):= sign(h); write(res,string fx,abs h);
     if i<n1 then write(res,<:  undetermined:>,nl,1)
             else write(res,string fx,y,nl,1)
   end;
   write(res,nl,2,<:w, V =:>,nl,1);
   l:= ii+1;
   for i:=1 step 1 until n do begin
     j:= cif(no(i)); if j>l then l:=j end;
   j:= l+4; i:= 72//j; l:= 72//i;
   fx:= format(l,ii,-0.9,fy); ii:= i-1; k:= 1;
   repwr2: m:= if n>k+ii then k+ii else n; write(res,nl,1);
   for j:=k step 1 until m do begin
     h:= A(j,n+1); i:= cif(no(j));
     write(res,string format(l,i,h,fy),h) end;
   write(res,nl,2);
   for i:=1 step 1 until n do begin
     for j:=k step 1 until m do
       write(res,string fx,A(i,j)*exp(i));
     write(res,nl,1)
   end;
   if m<n then begin k:= m+1; goto repwr2 end;
   for i:=n1 step 1 until n do
   for j:=1  step 1 until n do C(i,j):= A(i,j)*sgm(i)/A(j,n+1);
   for i:=1  step 1 until n do
   for j:=1  step 1 until i do
     A(j,i):= A(i,j):= sum(C(k,i)*C(k,j),k,n1,n)*S;
!,l2,r/C/A/,
l./if B/,i/
   j:= 1; if B then
   write(res,<:<12>:>,nl,2,string head(increase(j)),nl,2);
/,l./1-t/,r/*(1-t)/-t/,l./t*Sold/,r/*Sold//,
l./ortho/,d,i/
  svdecomp(a,A,sgm,dy,m,n,n1);
/,l./pm/,g1/A(j,n+1)/dy(j)/,f
i ftxt
▶EOF◀