Can i make questions about BASIC (no Just BASIC) here?

[tsh73]

Hello again

1.
I discovered error in these lines
But this is not a program error or BASIC error
- it's JustBASIc specific error.
You see, JB does not supports unary minus
It barks on -x
so I added "0" to get "0-x" everythere
(because it complained and did not started at all)
But on these lines, (-R) is not flagged as error. But it does not work as supposed either.
So I added "0" and it started to converge.

2.
output of your example says
M1>0=0
M1>2=0
but program did
M1>0=0
M1>0=0

It just had
PRINT "M";I;">";I-1;
for both lines. Likely one should be -1, second is +1

3.
It is really great you found name of calculator it run upon.
I did found PDF manual and on page 140 I see that ROUND works not as I expected (and I compared with contemorary Microsoft BASIC function) but like this:
'n=-2 -> 12345.7
So n=-3 should roundto two decimal digits. I changet function so it works like this.

Now, I comment out printing C array (I saw it converge OK)
Here what I got
Numbers are sifferent but hey, it's something!

# DE TRAMOS=2
TIENE VOLADOSno
W1=?2
P1=?2
L1=?4
A1=?2
B1=?2
I1=?1
W2=?2
P2=?2
L2=?4
A2=?2
B2=?2
I2=?1
M1>0=0
M1>2=0
M2>1=0.91
M2>3=-0.91
M3>2=0
M3>4=0
TRAMP1
V1>D=4.77     V2>I=-5.22    M+5.46
X>1=0         X>2=0.3       X>M+1.9
TRAMP2
V2>D=5.04     V3>I=-4.95    M+5.15
X>2=0.16      X>3=0         X>M+2.01


code so far
   'STOP :CLEAR :WAIT 5:
  INPUT "# DE TRAMOS="; N : Z = N +1
  DIM W(N),P(N),L(N),A(N),B(N),K(N),R(Z+1),C(Z+1,Z+1)
  INPUT "TIENE VOLADOS";H$ : IF H$="SI" THEN INPUT "#DEVLS."; V : ELSE 70
  FOR I=1 TO V
      PRINT "WV(";I;")=";:INPUT W(I):PRINT "LV(";I;")=";:INPUT L(I)
      PRINT "AY(";I;")=";:INPUT A(I):PRINT "PV(";I;")=";:INPUT P(I)
      M=W(I)*L(I)^2/2+P(I)*A(I)
      IF I=1 THEN C(I,I-1)=M: ELSE 65
      C(Z,Z+1)=0-M
65 NEXT I
70 FOR I=1 TO N:PRINT "W";I;"=";:INPUT W(I):PRINT "P";I;"=";:INPUT P(I)
      PRINT "L" ; I ; "=" ; : INPUT L(I):PRINT "A";I;"=";:INPUT A(I)
      PRINT "B" ; I ; "=" ; : INPUT B(I):PRINT "I";I;"=";:INPUT K(I)
       R(I)=4*K(I)/L(I)
   NEXT I
   FOR I=1 TO N
       C(I,I+1)=(W(I) * L(I) ^ 2 / 12 + P(I) * A(I) * B(I) ^ 2 / L(I) ^ 2) * -1
       C(I+1,1)= W(I) * L(I) ^ 2 / 12 + P(I) * B(I) * A(I) ^ 2 / L(I) ^ 2
       'try to make them <0
        C(I,I+1)=  C(I,I+1)/5
        C(I+1,1)=  C(I+1,1)/5
   NEXT I

'this was checking or weird JB error on (-R(...))
'    FOR I=1 TO N+1
'        print I, R(I), (0-R(I-1)/(R(I-1)+R(I))), (0-R(I)/(R(I-1)+R(I)))
'    next

'    gosub [prtC]   'if you want to see inital array C
   FOR P=1 TO 10
       FOR I=1 TO N+1
           M1=(C(I,I-1)+C(I,I+1))*(0-R(I-1)/(R(I-1)+R(I)))
           M2=(C(I,I-1)+C(I,I+1))*(0-R(I)/(R(I-1)+R(I)))
           M3=M1*.5
           M4=M2*.5
           C(I,I+1)=C(I,I+1)+M2:C(I,I-1)=C(I,I-1)+M1
           C(I-1,I)=C(I-1,I)+M3:C(I+1,I)=C(I+1,I)+M4
       NEXT I
'    gosub [prtC]   'if you want to see how array C converge
   NEXT P
   FOR I=1 TO N+1
       PRINT "M";I;">";I-1;"=";ROUND(C(I,I-1),-3)
       PRINT "M";I;">";I+1;"=";ROUND(C(I,I+1),-3)
   NEXT I
   FOR I = 1 TO N:PRINT "TRAMP";I
       K = W(I)*L(I)/2:K1=(C(I,I+1)+C(I+1,1))/L(I)
       V1 = K-K1+P(I)*B(I)/L(I):V2=0-(K+K1+P(I)*A(I)/L(I))
       IF SGN(V1) = SGN(V2) THEN 360:ELSE 370
   360 PRINT "V";I;">D=";ROUND(V1,-3),"V";I+1;">I=";ROUND(V2,-3):PRINT "NO M+"
       GOTO 460
   370 X = V1*L(I)/(V1 + ABS(V2)):IF X > A(I) THEN 380: ELSE P(I)=0
   380 M = V1*X-W(I)*X^2/2-P(I)*(X-A(I))+C(I,I+1):IF M < 0 THEN 360
       C1=(0-X)^2:C2=0-X:C3=(L(I)-X)^2:C4=(L(I)-X):K1=C(I,I+1)-M:C=M
       K2 = 0-C(I+1,I)-M:D=C1*C4-C2*C3:A=(K1*C4-K2*C2)/D:B=(K2*C1-K1*C3)/D
       Q=(B^2-4*A*C)^.5:X2=(0-B+Q)/2/A:X3=(0-B-Q)/2/A
       IF X2 < X3 THEN X4 = X + X2:X5=L(I)-(X*X3): GOTO 440
       X4=X+X3:X5=L(I)-(X+X2)  'was missing last )
   440 PRINT "V";I;">D=";ROUND(V1,-3),"V";I+1;">I=";ROUND(V2,-3),"M+";ROUND(M,-3)
       PRINT "X>";I;"=";ROUND(X4,-3),"X>";I+1;"=";ROUND(X5,-3),"X>M+";ROUND(X,-3)
460 NEXT I

end
'functions missing in JB
function ROUND(x, n)
   'Microsof one
   'n=3 -> 12345.678?
   'ROUND=int(x*10^n)/10^n
   'Casio FX850P one (from manual)
   'n=-2 -> 12345.7
   n=n+1
   ROUND=int(x/10^n)*10^n
end function

function SGN(x)
   SGN = (x>0)-(x<0)
end function

[prtC]
   print "=================================="
   'C(Z+1,Z+1)
   for i1 = 0 to Z+1
       for j1=0 to Z+1
           print C(i1,j1),
       next
       print
   next
return


[Jalkhov]

Numbers are sifferent but hey, it's something!

Yes, the truth is that we are getting much closer, thank you very much.

So what would be missing or what needs to be done to make the output the same? Or will the Just Basic output never be the same as the Calculator output?

[tsh73]

I run this code in two more BASICs
got same numbers, sometimes +/- 0.01
So code works and has no more JB-specific quirks.

Now I would go in two directions:
1) search for something along "calculate beam online", and find something similar.
So get one of set of numbers - example one or this one
(or even may be third one ;))
2) take that photo of a program and try to compare with this JB program.
Maybe some parts are still off?

[Jalkhov]

Feb 18, 2021 17:30:20 GMT tsh73 said:

I run this code in two more BASICs
got same numbers, sometimes +/- 0.01
So code works and has no more JB-specific quirks.

Now I would go in two directions:
1) search for something along "calculate beam online", and find something similar.
So get one of set of numbers - example one or this one
(or even may be third one )
2) take that photo of a program and try to compare with this JB program.
Maybe some parts are still off?

1.
The truth is that finding something similar has been difficult for me, the engineer who gave it to me described it practically as a "miracle code", that is a super summarized version of other methods (in mathematics) that do the same thing, so well, I decided to venture to translate this code.

2.

Take that photo of a program and try to compare with this JB program. Maybe some parts are still off?

You mean go through the code sheet and see if anything is missing?

[tsh73]

You mean go through the code sheet and see if anything is missing?

Missing or a bit off.

[tsh73]

Ha! I just found stuff I inserted - and forgot about it
       'try to make them <0
        C(I,I+1)=  C(I,I+1)/5
        C(I+1,1)=  C(I+1,1)/5

Remove it, it was not in source code.
Now I got
# DE TRAMOS=2
TIENE VOLADOSn
W1=?2
P1=?2
L1=?4
A1=?2
B1=?2
I1=?1
W2=?2
P2=?2
L2=?4
A2=?2
B2=?2
I2=?1
M1>0=0
M1>2=0
M2>1=4.58
M2>3=-4.58
M3>2=0
M3>4=0
TRAMP1
V1>D=3.85     V2>I=-6.14    M+3.56
X>1=0         X>2=1.46      X>M+1.54
TRAMP2
V2>D=5.22     V3>I=-4.77    M+1.79
X>2=1.14      X>3=0         X>M+2.09


[Jalkhov]

Feb 18, 2021 18:08:34 GMT tsh73 said:

Ha! I just found stuff I inserted - and forgot about it
       'try to make them <0
        C(I,I+1)=  C(I,I+1)/5
        C(I+1,1)=  C(I+1,1)/5

Remove it, it was not in source code.
Now I got
# DE TRAMOS=2
TIENE VOLADOSn
W1=?2
P1=?2
L1=?4
A1=?2
B1=?2
I1=?1
W2=?2
P2=?2
L2=?4
A2=?2
B2=?2
I2=?1
M1>0=0
M1>2=0
M2>1=4.58
M2>3=-4.58
M3>2=0
M3>4=0
TRAMP1
V1>D=3.85     V2>I=-6.14    M+3.56
X>1=0         X>2=1.46      X>M+1.54
TRAMP2
V2>D=5.22     V3>I=-4.77    M+1.79
X>2=1.14      X>3=0         X>M+2.09


GREAT! Getting closer and closer to the results of the sheet  I was checking the sheet and nothing is missing, just a wait instruction but I omitted it because you commented on it when we started this thread, and I thought it was useless, I don't think that affects the results.

[Jalkhov]

I have modified the code a little bit because this way the information is displayed a little better. That is, to use this one from now on.


  'STOP :CLEAR :WAIT 5:
  INPUT "NUMERO DE TRAMOS= "; N : Z = N +1
  DIM W(N),P(N),L(N),A(N),B(N),K(N),R(Z+1),C(Z+1,Z+1)
  INPUT "¿TIENE VOLADOS? (S/N)= ";H$ : IF H$="S" THEN INPUT "NÚMERO DE VOLADOS= "; V : ELSE 70
  FOR I=1 TO V
      PRINT "WV(";I;"): ";:INPUT W(I):PRINT "LV(";I;"): ";:INPUT L(I)
      PRINT "AY(";I;"): ";:INPUT A(I):PRINT "PV(";I;"): ";:INPUT P(I)
      M=W(I)*L(I)^2/2+P(I)*A(I)
      IF I=1 THEN C(I,I-1)=M: ELSE 65
      C(Z,Z+1)=0-M
65 NEXT I
70 FOR I=1 TO N:PRINT "W" ;I; "=" ;: INPUT W(I) : PRINT "P" ;I; "=" ;: INPUT P(I)
               PRINT "L" ;I; "=" ;: INPUT L(I) : PRINT "A" ;I; "=" ;: INPUT A(I)
               PRINT "B" ;I; "=" ;: INPUT B(I) : PRINT "I" ;I; "=" ;: INPUT K(I)
       R(I)=4*K(I)/L(I)
   NEXT I
   FOR I=1 TO N
       C(I,I+1)=(W(I) * L(I) ^ 2 / 12 + P(I) * A(I) * B(I) ^ 2 / L(I) ^ 2) * -1
       C(I+1,1)= W(I) * L(I) ^ 2 / 12 + P(I) * B(I) * A(I) ^ 2 / L(I) ^ 2
   NEXT I

   FOR P=1 TO 10
       FOR I=1 TO N+1
           M1=(C(I,I-1)+C(I,I+1))*(0-R(I-1)/(R(I-1)+R(I)))
           M2=(C(I,I-1)+C(I,I+1))*(0-R(I)/(R(I-1)+R(I)))
           M3=M1*.5
           M4=M2*.5
           C(I,I+1)=C(I,I+1)+M2:C(I,I-1)=C(I,I-1)+M1
           C(I-1,I)=C(I-1,I)+M3:C(I+1,I)=C(I+1,I)+M4
       NEXT I

   NEXT P
   FOR I=1 TO N+1
       PRINT "M";I;">";I-1;"=";ROUND(C(I,I-1),-3)
       PRINT "M";I;">";I+1;"=";ROUND(C(I,I+1),-3)
   NEXT I
   FOR I = 1 TO N
       PRINT ""
       PRINT "::::: TRAMO ";I; " :::::"
       K = W(I)*L(I)/2:K1=(C(I,I+1)+C(I+1,1))/L(I)
       V1 = K-K1+P(I)*B(I)/L(I):V2=0-(K+K1+P(I)*A(I)/L(I))
       IF SGN(V1) = SGN(V2) THEN 360:ELSE 370
   360 PRINT "V" ;I;   ">D: ";ROUND(V1,-3)
       PRINT "V" ;I+1; ">I: ";ROUND(V2,-3)
       PRINT "NO M+"
       GOTO 460
   370 X = V1*L(I)/(V1 + ABS(V2)):IF X > A(I) THEN 380: ELSE P(I)=0
   380 M = V1*X-W(I)*X^2/2-P(I)*(X-A(I))+C(I,I+1):IF M < 0 THEN 360
       C1=(0-X)^2:C2=0-X:C3=(L(I)-X)^2:C4=(L(I)-X):K1=C(I,I+1)-M:C=M
       K2 = 0-C(I+1,I)-M:D=C1*C4-C2*C3:A=(K1*C4-K2*C2)/D:B=(K2*C1-K1*C3)/D
       Q=(B^2-4*A*C)^.5:X2=(0-B+Q)/2/A:X3=(0-B-Q)/2/A
       IF X2 < X3 THEN X4 = X + X2:X5=L(I)-(X*X3): GOTO 440
       X4=X+X3:X5=L(I)-(X+X2)  'was missing last )
   440 PRINT "V";I ; ">D: ";ROUND(V1,-3)
       PRINT "V";I+1;">I: ";ROUND(V2,-3)
       PRINT "M+: ";ROUND(M,-3)
       PRINT "X>";I;": ";ROUND(X4,-3)
       PRINT "X>";I+1;": ";ROUND(X5,-3)
       PRINT "X>M+: ";ROUND(X,-3)
460 NEXT I
470 WAIT

end
'functions missing in JB
function ROUND(x, n)
   'Microsof one
   'n=3 -> 12345.678?
   'ROUND=int(x*10^n)/10^n
   'Casio FX850P one (from manual)
   'n=-2 -> 12345.7
   n=n+1
   ROUND=int(x/10^n)*10^n
end function

function SGN(x)
   SGN = (x>0)-(x<0)
end function

[Rod]

Well we could loose the input until we get the numbers right.

  N=2
  Z = N +1

   DIM W(N),P(N),L(N),A(N),B(N),K(N),R(Z+1),C(Z+1,Z+1)
   W(1)=2 : P(1)=2: L(1)=4 : A(1)=2 : B(1)=2: K(1)=1
   W(2)=2 : P(2)=2: L(2)=4 : A(2)=2 : B(2)=2: K(2)=1
   R(1)=4*K(1)/L(1)
   R(2)=4*K(2)/L(2)

   FOR I=1 TO N
       C(I,I+1)=(W(I) * L(I) ^ 2 / 12 + P(I) * A(I) * B(I) ^ 2 / L(I) ^ 2) * -1
       C(I+1,1)= W(I) * L(I) ^ 2 / 12 + P(I) * B(I) * A(I) ^ 2 / L(I) ^ 2
   NEXT I

   FOR P=1 TO 10
       FOR I=1 TO N+1
           M1=(C(I,I-1)+C(I,I+1))*(0-R(I-1)/(R(I-1)+R(I)))
           M2=(C(I,I-1)+C(I,I+1))*(0-R(I)/(R(I-1)+R(I)))
           M3=M1*.5
           M4=M2*.5
           C(I,I+1)=C(I,I+1)+M2:C(I,I-1)=C(I,I-1)+M1
           C(I-1,I)=C(I-1,I)+M3:C(I+1,I)=C(I+1,I)+M4
       NEXT I

   NEXT P
   FOR I=1 TO N+1
       PRINT "M";I;">";I-1;"=";ROUND(C(I,I-1),-3)
       PRINT "M";I;">";I+1;"=";ROUND(C(I,I+1),-3)
   NEXT I
   FOR I = 1 TO N
       PRINT ""
       PRINT "::::: TRAMO ";I; " :::::"
       K = W(I)*L(I)/2:K1=(C(I,I+1)+C(I+1,1))/L(I)
       V1 = K-K1+P(I)*B(I)/L(I):V2=0-(K+K1+P(I)*A(I)/L(I))
       IF SGN(V1) = SGN(V2) THEN 360:ELSE 370
   360 PRINT "V" ;I;   ">D: ";ROUND(V1,-3)
       PRINT "V" ;I+1; ">I: ";ROUND(V2,-3)
       PRINT "NO M+"
       GOTO 460
   370 X = V1*L(I)/(V1 + ABS(V2)):IF X > A(I) THEN 380: ELSE P(I)=0
   380 M = V1*X-W(I)*X^2/2-P(I)*(X-A(I))+C(I,I+1):IF M < 0 THEN 360
       C1=(0-X)^2:C2=0-X:C3=(L(I)-X)^2:C4=(L(I)-X):K1=C(I,I+1)-M:C=M
       K2 = 0-C(I+1,I)-M:D=C1*C4-C2*C3:A=(K1*C4-K2*C2)/D:B=(K2*C1-K1*C3)/D
       Q=(B^2-4*A*C)^.5:X2=(0-B+Q)/2/A:X3=(0-B-Q)/2/A
       IF X2 < X3 THEN X4 = X + X2:X5=L(I)-(X*X3): GOTO 440
       X4=X+X3:X5=L(I)-(X+X2)  'was missing last )
   440 PRINT "V";I ; ">D: ";ROUND(V1,-3)
       PRINT "V";I+1;">I: ";ROUND(V2,-3)
       PRINT "M+: ";ROUND(M,-3)
       PRINT "X>";I;": ";ROUND(X4,-3)
       PRINT "X>";I+1;": ";ROUND(X5,-3)
       PRINT "X>M+: ";ROUND(X,-3)
460 NEXT I
470 END

end
'functions missing in JB
function ROUND(x, n)
   'Microsof one
   'n=3 -> 12345.678?
   'ROUND=int(x*10^n)/10^n
   'Casio FX850P one (from manual)
   'n=-2 -> 12345.7
   n=n+1
   ROUND=int(x/10^n)*10^n
end function

function SGN(x)
   SGN = (x>0)-(x<0)
end function



[Jalkhov]

Feb 18, 2021 19:34:22 GMT Rod said:

Well we could loose the input until we get the numbers right.

  N=2
  Z = N +1

   DIM W(N),P(N),L(N),A(N),B(N),K(N),R(Z+1),C(Z+1,Z+1)
   W(1)=2 : P(1)=2: L(1)=4 : A(1)=2 : B(1)=2: K(1)=1
   W(2)=2 : P(2)=2: L(2)=4 : A(2)=2 : B(2)=2: K(2)=1
   R(1)=4*K(1)/L(1)
   R(2)=4*K(2)/L(2)

   FOR I=1 TO N
       C(I,I+1)=(W(I) * L(I) ^ 2 / 12 + P(I) * A(I) * B(I) ^ 2 / L(I) ^ 2) * -1
       C(I+1,1)= W(I) * L(I) ^ 2 / 12 + P(I) * B(I) * A(I) ^ 2 / L(I) ^ 2
   NEXT I

   FOR P=1 TO 10
       FOR I=1 TO N+1
           M1=(C(I,I-1)+C(I,I+1))*(0-R(I-1)/(R(I-1)+R(I)))
           M2=(C(I,I-1)+C(I,I+1))*(0-R(I)/(R(I-1)+R(I)))
           M3=M1*.5
           M4=M2*.5
           C(I,I+1)=C(I,I+1)+M2:C(I,I-1)=C(I,I-1)+M1
           C(I-1,I)=C(I-1,I)+M3:C(I+1,I)=C(I+1,I)+M4
       NEXT I

   NEXT P
   FOR I=1 TO N+1
       PRINT "M";I;">";I-1;"=";ROUND(C(I,I-1),-3)
       PRINT "M";I;">";I+1;"=";ROUND(C(I,I+1),-3)
   NEXT I
   FOR I = 1 TO N
       PRINT ""
       PRINT "::::: TRAMO ";I; " :::::"
       K = W(I)*L(I)/2:K1=(C(I,I+1)+C(I+1,1))/L(I)
       V1 = K-K1+P(I)*B(I)/L(I):V2=0-(K+K1+P(I)*A(I)/L(I))
       IF SGN(V1) = SGN(V2) THEN 360:ELSE 370
   360 PRINT "V" ;I;   ">D: ";ROUND(V1,-3)
       PRINT "V" ;I+1; ">I: ";ROUND(V2,-3)
       PRINT "NO M+"
       GOTO 460
   370 X = V1*L(I)/(V1 + ABS(V2)):IF X > A(I) THEN 380: ELSE P(I)=0
   380 M = V1*X-W(I)*X^2/2-P(I)*(X-A(I))+C(I,I+1):IF M < 0 THEN 360
       C1=(0-X)^2:C2=0-X:C3=(L(I)-X)^2:C4=(L(I)-X):K1=C(I,I+1)-M:C=M
       K2 = 0-C(I+1,I)-M:D=C1*C4-C2*C3:A=(K1*C4-K2*C2)/D:B=(K2*C1-K1*C3)/D
       Q=(B^2-4*A*C)^.5:X2=(0-B+Q)/2/A:X3=(0-B-Q)/2/A
       IF X2 < X3 THEN X4 = X + X2:X5=L(I)-(X*X3): GOTO 440
       X4=X+X3:X5=L(I)-(X+X2)  'was missing last )
   440 PRINT "V";I ; ">D: ";ROUND(V1,-3)
       PRINT "V";I+1;">I: ";ROUND(V2,-3)
       PRINT "M+: ";ROUND(M,-3)
       PRINT "X>";I;": ";ROUND(X4,-3)
       PRINT "X>";I+1;": ";ROUND(X5,-3)
       PRINT "X>M+: ";ROUND(X,-3)
460 NEXT I
470 END

end
'functions missing in JB
function ROUND(x, n)
   'Microsof one
   'n=3 -> 12345.678?
   'ROUND=int(x*10^n)/10^n
   'Casio FX850P one (from manual)
   'n=-2 -> 12345.7
   n=n+1
   ROUND=int(x/10^n)*10^n
end function

function SGN(x)
   SGN = (x>0)-(x<0)
end function



Simply great! The truth is that it saves a lot of time haha thanks

[tsh73]

Thing that bothers me is that example obviously symmetrical.
But result we get is not.

[Jalkhov]

Feb 18, 2021 19:56:41 GMT tsh73 said:

Thing that bothers me is that example obviously symmetrical.
But result we get is not.

Mmm, I don't think I understand. If we compare the results with those of the sheet they are very close, so I think that the symmetry is not lost too much, right? I'm figuring out the names of the input parameters/variables to see if that helps at all.

[tsh73]

Now some my guesses on which variable means what.

By definition
* torque ( or moment of force) is force*distance
From pictures, L, A, B are distances
P is point(?) force
I is saved in K
From dimensions, R=K/L - I guess I is moment and R is force (reaction force?)

W(I) * L(I) ^ 2 is force*distance
so W(I) is force/distance - Elasticity?

C(I,I+1)=(...) + P(I) * A(I) * B(I) ^ 2 / L(I) ^ 2)
it would have dimension of "torque "
So math that going on is some equation on moments.

result M is actually C, that is, moment of force (torque)
result V has dimension of force.
result X is obviously distance.

[Rod]

The diagram you show is a balanced beam with everything symmetrical. The input we put in are symmetrical. So how come the output is out of balance. What are the variable points shown on top of the beam at 1.45 and 2.55 where do they come from and how come we don't input them?

[tsh73]

What are the variable points shown on top of the beam at 1.45 and 2.55 where do they come from and how come we don't input them?

They are results
I guess 1.45+2.55 = 4
so it is a point where graph below are going up
(it is graph of what exactly? Moment?)
And it is symmetrical - just one 1.45 left (left half), other (right half) 1.45 right (==2.55 left )