CLS : DEFDBL A-Z
DeltaT = 5: IntNum = 24: Sigma = 100
Sum0 = 0: Sum1 = 0: Sum2 = 0: Sum3 = 0: Sum4 = 0
BSum0 = 0: BSum1 = 0: BSum2 = 0
DIM Y(IntNum)
OPEN "Caven.Dat" FOR OUTPUT AS #1
   PRINT #1, "Dati sperimentali letti sull'asta graduata relativi all'esperimento"
   PRINT #1, "della bilancia di Cavendish durante i primi"; DeltaT * IntNum; "secondi:"
   PRINT #1, ""
   PRINT #1, "Tempi (in secondi)", "Valore raggiunto (in centimetri)"
CLOSE
FOR I = 0 TO IntNum
   OPEN "Caven.Dat" FOR APPEND AS #1
      PRINT "Valore raggiunto dopo"; DeltaT * I; "secondi";
      INPUT ""; Y(I)
      PRINT #1, SPACE$(18 - LEN(STR$(DeltaT * I))); DeltaT * I; SPACE$(41 - LEN(STR$(Y(I)))); Y(I)
      Sum0 = Sum0 + Sigma
      Sum1 = Sum1 + Sigma * DeltaT * I
      Sum2 = Sum2 + Sigma * (DeltaT * I) ^ 2
      Sum3 = Sum3 + Sigma * (DeltaT * I) ^ 3
      Sum4 = Sum4 + Sigma * (DeltaT * I) ^ 4
      BSum0 = BSum0 + Sigma * Y(I)
      BSum1 = BSum1 + Sigma * DeltaT * I * Y(I)
      BSum2 = BSum2 + Sigma * (DeltaT * I) ^ 2 * Y(I)
   CLOSE
NEXT
ID = (Sum0 * Sum2 * Sum4 + 2 * Sum1 * Sum2 * Sum3 - Sum1 ^ 2 * Sum4 - Sum2 ^ 3 - Sum0 * Sum3 ^ 2) ^ -1
D1 = Sum2 * Sum4 - Sum3 ^ 2: D2 = Sum2 * Sum3 - Sum1 * Sum4: D3 = Sum1 * Sum3 - Sum2 ^ 2
D4 = Sum0 * Sum4 - Sum2 ^ 2: D5 = Sum1 * Sum2 - Sum0 * Sum3: D6 = Sum0 * Sum2 - Sum1 ^ 2
A = ID * (D1 * BSum0 + D2 * BSum1 + D3 * BSum2)
B = ID * (D2 * BSum0 + D4 * BSum1 + D5 * BSum2)
C = ID * (D3 * BSum0 + D5 * BSum1 + D6 * BSum2)
SigmaA = ID * (Sum2 * Sum4 - Sum3 ^ 2)
SigmaB = ID * (Sum0 * Sum4 - Sum2 ^ 2)
SigmaC = ID * (Sum0 * Sum2 - Sum1 ^ 2)
IR = 1 / D6
AR = IR * (Sum2 * BSum0 - Sum1 * BSum1): BR = IR * (Sum0 * BSum1 - Sum1 * BSum0)
SigmaAR = IR * Sum2: SigmaBR = IR * Sum0
OPEN "Caven.Dat" FOR APPEND AS #1
   PRINT #1, ""
   PRINT #1, "DeltaT="; DeltaT; "Sigma^-2="; Sigma
   PRINT #1, ""
   PRINT #1, "Parametri della retta che meglio approssima l'andamento empirico:"
   PRINT #1, "Sum0="; Sum0; "Sum1="; Sum1; "Sum2="; Sum2
   PRINT #1, "BSum0="; BSum0; "BSum1="; BSum1
   PRINT #1, "InvDet(D)="; IR
   PRINT #1, "D1="; IR * Sum2; "D2="; -IR * Sum1; "D3="; IR * Sum0
   PRINT #1, "A="; AR; "B="; BR
   PRINT #1, "SigmaA="; SigmaAR; "SigmaB="; SigmaBR
   PRINT #1, ""
   PRINT #1, "Parametri della parabola che meglio approssima l'andamento empirico:"
   PRINT #1, "Sum0="; Sum0; "Sum1="; Sum1; "Sum2="; Sum2; "Sum3="; Sum3; "Sum4="; Sum4
   PRINT #1, "BSum0="; BSum0; "BSum1="; BSum1; "BSum2="; BSum2
   PRINT #1, "InvDet(D)="; ID
   PRINT #1, "D1="; D1; "D2="; D2; "D3="; D3
   PRINT #1, "D4="; D4; "D5="; D5; "D6="; D6
   PRINT #1, "A="; A; "B="; B; "C="; C
   PRINT #1, "SigmaA="; SigmaA
   PRINT #1, "SigmaB="; SigmaB; "SigmaC="; SigmaC
CLOSE
PRINT : PRINT "DeltaT="; DeltaT; "Sigma^-2="; Sigma
PRINT : PRINT "Parametri della retta che meglio approssima l'andamento empirico:"
PRINT "Sum0="; Sum0; "Sum1="; Sum1; "Sum2="; Sum2
PRINT "BSum0="; BSum0; "BSum1="; BSum1
PRINT "InvDet(D)="; IR
PRINT "D1="; IR * Sum2; "D2="; -IR * Sum1; "D3="; IR * Sum0
PRINT "A="; AR; "B="; BR
PRINT "SigmaA="; SigmaAR; "SigmaB="; SigmaBR
PRINT : PRINT "Parametri della parabola che meglio approssima l'andamento empirico:"
PRINT "Sum0="; Sum0; "Sum1="; Sum1; "Sum2="; Sum2; "Sum3="; Sum3; "Sum4="; Sum4
PRINT "BSum0="; BSum0; "BSum1="; BSum1; "BSum2="; BSum2
PRINT "InvDet(D)="; ID
PRINT "D1="; D1; "D2="; D2; "D3="; D3
PRINT "D4="; D4; "D5="; D5; "D6="; D6
PRINT "A="; A; "B="; B; "C="; C
PRINT "SigmaA="; SigmaA
PRINT "SigmaB="; SigmaB; "SigmaC="; SigmaC