LET I = I + di

25 LET I = I + di

LET Uc = Uc + dUc

LET X = X + dx

LET V = V + dv

LET integral2 = integral2 + ABS(I) * h

LET integral = integral + h * (ABS(I)) ^ 2

LET a = (Bo * Pi * D * w * I - nu * V - ksi * X) / m1

LET Ua = La * di / h + Bo * Pi * D * w * V + I * Ra

LET tall = tall + h

26 IF q = 1 THEN

LOCATE 1, 1: PRINT "t="; tall * 1000, " "

LOCATE 2, 1: PRINT "I="; I, " "

LOCATE 3, 1: PRINT "Uб="; Uc, " "

LOCATE 4, 1: PRINT "X="; X * 1000, " "

LOCATE 5, 1: PRINT "V="; V * 1000, " "

LOCATE 6, 1: PRINT "a="; a, " "

LOCATE 7, 1: PRINT "Ua="; Ua, " "

END IF

PSET (t * 20000 * .01 / t2, 250 - I * .5)

PSET (t * 20000 * .01 / t2, 250 + Uc * .1), 3

PSET (t * 20000 * .01 / t2, 250)

PSET (t * 20000 * .01 / t2, 250 - V * 100), 4

PSET (t * 20000 * .01 / t2, 250 - X * 100000), 5

PSET (t * 20000 * .01 / t2, 250 - a / 3), 2

PSET (t * 20000 * .01 / t2, 250 - Ua * .1), 6

LET t = t + h

IF t > ti + t1 + t56 THEN GOTO 27

GOTO 20

27i5 = I

V5 = V

x5 = X

28LET k1i = -E / La - (R / La) * i5 - (Bo * Pi * D * w * V / La)

LET i5 = I + h * .5 * k1i

LET k2i = -E / La - (R / La) * i5 - (Bo * Pi * D * w * V / La)

LET i5 = I + k2i * h * .5

LET k3i = -E / La - (R / La) * i5 - (Bo * Pi * D * w * V / La)

LET i5 = I + k3i * h

LET k4i = -E / La - (R / La) * i5 - (Bo * Pi * D * w * V / La)

LET di = h * (2 * k2i + k1i + 2 * k3i + k4i) / 6

29 LET k1x = V

LET x5 = X + k1x * h * .5

LET k2x = V

LET x5 = X + k2x * h * .5

LET k3x = V

LET x5 = X + k3x * h

LET k4x = V

LET dx = (k1x + 2 * k2x + 2 * k3x + k4x) * h / 6

30 LET k1V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V5 - (ksi / m1) * X

LET V5 = V + k1V * h * .5

LET k2V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V5 - (ksi / m1) * X

LET V5 = V + k2V * h * .5

LET k3V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V5 - (ksi / m1) * X

LET V5 = V + k3V * h

LET k4V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V5 - (ksi / m1) * X

LET dv = h * (2 * k2V + k1V + 2 * k3V + k4V) / 6

31 LET I = I + di

LET X = X + dx

LET V = V + dv

LET integral2 = integral2 + ABS(I) * h

LET integral = integral + h * (ABS(I)) ^ 2

LET a = (Bo * Pi * D * w * I - nu * V - ksi * X) / m1

LET Ua = La * di / h + Bo * Pi * D * w * V + I * Ra

LET tall = tall + h

32 IF q = 1 THEN

LOCATE 1, 1: PRINT "t="; tall * 1000, " "

LOCATE 2, 1: PRINT "I="; I, " "

LOCATE 3, 1: PRINT "Uc="; Uc, " "

LOCATE 4, 1: PRINT "X="; X * 1000, " "

LOCATE 5, 1: PRINT "V="; V * 1000, " "

LOCATE 6, 1: PRINT "a="; a, " "

LOCATE 7, 1: PRINT "Ua="; Ua, " "

END IF

PSET (t * 20000 * .01 / t2, 250 - I * .5)

PSET (t * 20000 * .01 / t2, 250 + Uc * .1), 3

PSET (t * 20000 * .01 / t2, 250)

PSET (t * 20000 * .01 / t2, 250 - V * 100), 4

PSET (t * 20000 * .01 / t2, 250 - X * 100000), 5

PSET (t * 20000 * .01 / t2, 250 - a / 3), 2

PSET (t * 20000 * .01 / t2, 250 - Ua * .1), 6

LET t = t + h

IF t > (ti + t1 + t2 + t56) THEN GOTO 33

GOTO 27

33 LET Uc = -Uc

331 LET Uc6 = Uc

LET i6 = I

LET x6 = X

LET V6 = V

34 LET k1i = E / La - (R / La) * i6 - (Uc / La) - (Bo * Pi * D * w / La) * V

LET i6 = I + h * .5 * k1i

LET k2i = E / La - (R / La) * i6 - (Uc / La) - (Bo * Pi * D * w / La) * V

LET i6 = I + k2i * h * .5

LET k3i = E / La - (R / La) * i6 - (Uc / La) - (Bo * Pi * D * w / La) * V

LET i6 = I + k3i * h

LET k4i = E / La - (R / La) * i6 - (Uc / La) - (Bo * Pi * D * w / La) * V

LET di = h * (2 * k2i + k1i + 2 * k3i + k4i) / 6

35 LET k1Uc = I / C

LET Uc6 = Uc + h * .5 * k1Uc

LET k2Uc = I / C

LET Uc6 = Uc + h * .5 * k2Uc

LET k3Uc = I / C

LET Uc6 = Uc + h * k3Uc

LET k4Uc = I / C

LET dUc = h * (2 * k2Uc + k1Uc + 2 * k3Uc + k4Uc) / 6

36 LET k1x = V

LET x6 = X + k1x * h * .5

LET k2x = V

LET x6 = X + k2x * h * .5

LET k3x = V

LET x6 = X + k3x * h

LET k4x = V

LET dx = (k1x + 2 * k2x + 2 * k3x + k4x) * h / 6

37 LET k1V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V6 - (ksi / m1) * X

LET V6 = V + k1V * h * .5

LET k2V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V6 - (ksi / m1) * X

LET V6 = V + k2V * h * .5

LET k3V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V6 - (ksi / m1) * X

LET V6 = V + k3V * h

LET k4V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V6 - (ksi / m1) * X

LET dv = h * (2 * k2V + k1V + 2 * k3V + k4V) / 6