Как напечатать вывод.bas файла в текст

Я пытаюсь распечатать результаты координат программы в текстовый файл, чтобы использовать ее в другой программе, но я ничего не знаю о GWBASIC и его первом использовании MS-DOS. Мне нужно, чтобы он открыл текстовый файл с именем plot.txt и распечатал там вывод и сохранил его без фактического построения графиков на GWBASIC. Вот программа, которую я нашел в старом журнале.

810 REM   MAKE A GLOBULAR
12 REM
14 R0=20: R2=R0*R0: R3=R2*R0
16 P1=3.14159265#
18 C0=P1*P1*R3/4
20 R1=R0/SQR(2)
22 XM=512: YM=512
24 X2=XM/2: Y2=YM/2: S=5
26 INPUT "HOW MANY STARS ";T
27 RANDOMIZE TIMER
28 CLS: REM  CLEAR SCREEN
30 FOR I=1 TO T
32 C=C0*RND: R=R1
34 REM
36 REM   NOW FIND R
38 FOR K=1 TO 5
40 GOSUB 100
42 R=R+(C-C1)/D
44 NEXT K
46 REM  3-DIMENSIONAL PLACE
48 X=RND-.5
50 Y=RND-.5
52 Z=RND-.5
54 S1=SQR(X*X+Y*Y+Z*Z)
56 IF S1>.5 THEN GOTO 48
58 REM  POINT IS NOW IN SPHERE
60 R=R*S1: X=X*R: Y=Y*R: Z=Z*R
62 GOSUB 200
64 NEXT I
66 END
68 REM
100 REM  NEWTON-RAPHSON ITERATION
105 A=R/R0
110 C1=ATN(A)*.5*R3
115 A=1+A*A
120 C1=C1+R*.5*R2/A
125 C1=P1*(C1-R*R2/(A*A))
130 D=4*P1*R*R/(A*A*A)
135 RETURN
140 REM
200 REM  2-DIMENSIONAL PLOT
203 SCREEN 9
205 X=X*S+X2: Y=Y*S+Y2
210 IF X<0 OR Y<0 THEN 225
215 IF X>=XM OR Y>=YM THEN 225
220 PSET(X,Y)
225 RETURN
230 REM  ------------------------
240 REM  APPEARED IN ASTRONOMICAL
250 REM  COMPUTING, SKY & TELE-
260 REM  SCOPE, APRIL, 1986
270 REM  ------------------------
Источник

2 ответа

Решение

Вот парафраз Python 3:

#globular.py
#Python paraphrase of model.bas from
#http://www.skyandtelescope.com/wp-content/uploads/model.bas

from math import pi, sqrt, atan
from random import uniform, random

#Global variables:

r0 = 20.0
r2 = r0**2
r3 = r0**3
c0 = pi**2*r3/4
r1 = r0/sqrt(2)

def NRI(c,r):
    #Newton-Raphson Iteration
    a = r/r0
    c1 = atan(a)*0.5*r3
    a = 1+a**2
    c1 += r*0.5*r2/a
    c1 = pi*(c1-r*r2/a**2)
    d = 4*pi*r**2/a**3
    return (c1,d)

def makeStars(t):
    stars = []
    for i in range(t):
        c = c0*random()
        r = r1
        for k in range(5):
            c1,d = NRI(c,r)
            r += (c-c1)/d
        while True:
            x = uniform(-0.5,0.5)
            y = uniform(-0.5,0.5)
            z = uniform(-0.5,0.5)
            s1 = sqrt(x**2 + y**2 + z**2)
            if s1 <= 0.5: break
        r *= s1
        x *= r
        y *= r
        z *= r
        stars.append((x,y,z))
    return stars

def starsToFile(t,fname):
    stars = makeStars(t)
    f = open(fname,'w')
    for star in stars:
        print(*star, sep = ', ',file = f)
    f.close()

Я пропустил часть о печати x а также y и вместо этого написал функцию makeStars вернуть список (x,y,z) кортежи, а также связанная функция, которая принимает такой вывод и отправляет его в текстовый файл. Эта последняя функция - единственное, что использовало Python 3 вместо Python 2. Если вы используете Python 2, вы можете импортировать функцию печати Python 3 из будущего.

Typing starsToFile(100,'stars.txt') в оболочке Python дал мне текстовый файл, который начинается:

-0.32838465248713156, -0.3294895266926551, -1.2963580524762535
14.20224408569865, 1.4434961933043464, 6.450969593697097
1.6525937589658193, -0.24447292610082685, 1.0543647986350608
1.5707528567123823, 5.190972598268825, -2.0054790217091134

У меня нет хороших трехмерных графиков рассеяния на кончиках моих пальцев, но вот снимок экрана с 50 точками, сгенерированными функцией и нанесенными с использованием системы компьютерной алгебры под названием Derive:

Последнее замечание: мне интересно, есть ли опечатка в исходном коде. Линия

C0=P1*P1*R3/4

кажется мне подозрительным, так как это довольно редко в математике для pi казаться в квадрате - хотя это случается. Может быть, должен быть только 1 фактор pi там (что будет иметь эффект установки C0 пропорционально объему сферы радиуса R0). С другой стороны, я не знаю точно, что здесь происходит, поэтому я оставил это. Если результаты кажутся проблематичными, вы можете поэкспериментировать с этой строкой.

Источник

Если вы хотите получить копию рассчитанных координат, просто добавьте эти строки:

1 OPEN "PLOT.TXT" FOR OUTPUT AS #1
65 CLOSE #1
221 PRINT #1, X + "," + Y

Программа будет работать как прежде, но в дополнение к этому она выводит координаты в файл с именем plot.txt

Поместите их в изображение размером 640x350 (этот размер требуется SCREEN 9) и вы получите тот же результат.

Источник
Другие вопросы по тегам