Private Declare Function SetPixel Lib "gdi32" (ByVal hdc As Long, ByVal x As Long, ByVal y As Long, ByVal crColor As Long) As Long
Const m = 1000
Const mult = 400
Private x(m, 10), y(m, 10), z(m, 10) As Single
Private xr(m), yr(m), f(m), zr(m) As Integer
Private z0, y0, x0, zs, ax, ay, az, n, xmax, ymax, t, t1, st, st0, st1 As Single
Private typ, n_t
Private Sub Form_Activate()
xmax = ScaleWidth
ymax = ScaleHeight
n = 1: t = 1: t1 = 1
x0 = 1: y0 = 0: z0 = -20
zs = 2
ax = 0: ay = 0: az = 0
st0 = 100
st = 0
typ = 0
n_t = 6
''''''''''''''''''''''''''''''''''''''''
n = 1
For i = 0 To 5 Step 0.2
x(n, t) = 0: y(n, t) = i: z(n, t) = 0
n = n + 1
x(n, t) = -i: y(n, t) = -i: z(n, t) = i
n = n + 1
x(n, t) = i: y(n, t) = -i: z(n, t) = i
n = n + 1
x(n, t) = i: y(n, t) = i: z(n, t) = i
n = n + 1
x(n, t) = i: y(n, t) = i: z(n, t) = -i
n = n + 1
x(n, t) = i: y(n, t) = -i: z(n, t) = -i
n = n + 1
x(n, t) = -i: y(n, t) = -i: z(n, t) = -i
n = n + 1
x(n, t) = -i: y(n, t) = i: z(n, t) = -i
n = n + 1
x(n, t) = -i: y(n, t) = -i: z(n, t) = -i
n = n + 1
Next
'''''''''''''''''''''''''''''''''''''''
t = 2
For i = 1 To 300
x(i, t) = Int(Rnd * 20) - 10
y(i, t) = Int(Rnd * 20) - 10
z(i, t) = Int(Rnd * 20) - 10
Next
Print n
''''''''''''''''''''''''''''''''''''
t = 3
n = 1
For i = -5 To 5 Step 0.5
x(n, t) = i: y(n, t) = -5: z(n, t) = 5
n = n + 1
x(n, t) = i: y(n, t) = -5: z(n, t) = -5
n = n + 1
x(n, t) = i: y(n, t) = 5: z(n, t) = 5
n = n + 1
x(n, t) = i: y(n, t) = 5: z(n, t) = -5
n = n + 1
x(n, t) = 5: y(n, t) = i: z(n, t) = 5
n = n + 1
x(n, t) = 5: y(n, t) = i: z(n, t) = -5
n = n + 1
x(n, t) = -5: y(n, t) = i: z(n, t) = 5
n = n + 1
x(n, t) = -5: y(n, t) = i: z(n, t) = -5
n = n + 1
x(n, t) = 5: y(n, t) = 5: z(n, t) = i
n = n + 1
x(n, t) = -5: y(n, t) = 5: z(n, t) = i
n = n + 1
x(n, t) = 5: y(n, t) = -5: z(n, t) = i
n = n + 1
x(n, t) = -5: y(n, t) = -5: z(n, t) = i
n = n + 1
Next
'''''''''''''''''''''''''''''''''''''''''
n = 1
t = 4
For i = 0 To 6 - 1
For j = 0 To 6
 Call pset_line(Cos(i * 6.28 / 6), Sin(i * 6.28 / 6), Cos(j * 6.28 / 6), Sin(j * 6.28 / 6), 6, 5)
Next
Next

'''''''''''''''''''''''''''''''''''''''''
n = 1
t = 5
For i = 0 To 10 Step 0.3
x(n, t) = i - 5: y(n, t) = -5: z(n, t) = -5
n = n + 1
x(n, t) = i - 5: y(n, t) = -5: z(n, t) = 5
n = n + 1
x(n, t) = 5: y(n, t) = -5: z(n, t) = i - 5
n = n + 1
x(n, t) = -5: y(n, t) = -5: z(n, t) = i - 5
n = n + 1
x(n, t) = 5 - i / 2: y(n, t) = i - 5: z(n, t) = i / 2 - 5
n = n + 1
x(n, t) = -5 + i / 2: y(n, t) = i - 5: z(n, t) = i / 2 - 5
n = n + 1
x(n, t) = -5 + i / 2: y(n, t) = i - 5: z(n, t) = -i / 2 + 5
n = n + 1
x(n, t) = 5 - i / 2: y(n, t) = i - 5: z(n, t) = -i / 2 + 5
n = n + 1

Next
''''''''''''''''''''''''''''''''''''''''''''
n = 2
t = 6
mm = 2
Call pset_line(-5, -2, -5, 2, 30, mm) 'n
Call pset_line(-5, 2, -1, -2, 30, mm)
Call pset_line(-1, -2, -1, 2, 30, mm)
Call pset_line(0, -2, 0, 2, 30, mm) 'i
Call pset_line(1, -2, 1, 2, 30, mm) 'c
Call pset_line(1, -2, 4, -2, 30, mm) 'c
Call pset_line(1, 2, 4, 2, 30, mm) 'c

n = 300
Call calc
1
Call erased
Call calc
Call pointed
DoEvents
GoTo 1
End Sub
Sub pset_line(x1, y1, x2, y2, stt, mm)
For k = 0 To stt
 x(n, t) = (x1 + (x2 - x1) * k / stt) * mm
 y(n, t) = (y1 + (y2 - y1) * k / stt) * mm
 n = n + 1
Next
End Sub
Sub calc()
For i = 1 To n
If typ = 0 Then
XX = x(i, t) + (x(i, t1) - x(i, t)) * st / st0
YY = y(i, t) + (y(i, t1) - y(i, t)) * st / st0
zz = z(i, t) + (z(i, t1) - z(i, t)) * st / st0
Else
XX = x(i, typ)
YY = y(i, typ)
zz = z(i, typ)
End If
 'GoTo 1
 S = Sqr((XX) ^ 2 + (YY) ^ 2)
 If XX = 0 Then
 XX = S * Cos(ax - 3.14 * (YY < 0) + 1.57)
 YY = S * Sin(ay - 3.14 * (YY < 0) + 1.57)
 Else
 Alfa = Atn(YY / XX)
 YY = S * Sin(ay - 3.14 * (XX < 0) + Alfa)
 XX = S * Cos(ax - 3.14 * (XX < 0) + Alfa)
 End If
1
 'GoTo 1
 S = Sqr((XX) ^ 2 + (zz) ^ 2)
 If XX = 0 Then
 XX = S * Cos(ax - 3.14 * (zz < 0) + 1.57)
 zz = S * Sin(az - 3.14 * (zz < 0) + 1.57)
 Else
 Alfa = Atn(zz / XX)
 zz = S * Sin(az + Alfa - 3.14 * (XX < 0))
 XX = S * Cos(ax + Alfa - 3.14 * (XX < 0))
End If

 xr(i) = Int((zs) * (XX - x0) / (zz - z0) * mult + xmax / 2)
 yr(i) = Int(ymax / 2 - (zs) * (YY - y0) / (zz - z0) * mult)
 zr(i) = zz - z0 - zs
Next
 ax = ax + 0.01
 ay = ay + 0.01
 az = az + 0.01
 If st < st0 Then
 st = st + 1
 Else
 If st1 < st0 Then
 st1 = st1 + 1
 Else
 t = t1: t1 = Int(Rnd * n_t) + 1: st = 0: st1 = 0
 End If
 End If
End Sub
Sub pointed()
For i = 1 To n
PSet (xr(i), yr(i)), &HFFFFFF - (zr(i) - 10) * (256 + 65536 + 1) * 10
'PSet (xr(i) + 1, yr(i)), &HFFFFFF - (zr(i) - 10) * (256 + 65536 + 1) * 10
'PSet (xr(i), yr(i) + 1), &HFFFFFF - (zr(i) - 10) * (256 + 65536 + 1) * 10

'retval = SetPixel(Form1.hdc, xr(i), yr(i), &HFFFFFF)
retval = SetPixel(Form1.hdc, xr(i) + 1, yr(i), &HFFFFFF)
retval = SetPixel(Form1.hdc, xr(i), yr(i) + 1, &HFFFFFF)

Next
DoEvents
End Sub
Sub erased()
For i = 1 To n
retval = SetPixel(Form1.hdc, xr(i), yr(i), BackColor)
retval = SetPixel(Form1.hdc, xr(i) + 1, yr(i), BackColor)
retval = SetPixel(Form1.hdc, xr(i), yr(i) + 1, BackColor)
'PSet (xr(i), yr(i)), BackColor
'PSet (xr(i) + 1, yr(i)), BackColor
'PSet (xr(i), yr(i) + 1), BackColor
Next
End Sub

Private Sub Form_KeyPress(KeyAscii As Integer)
If KeyAscii = 27 Then End
End Sub

после прорисовки random-ом звёзд надо чтобы они стали сходиться или же расходиться из определённой точки экрана

а потом по таймеру начать каждую двигать по прямой к центру