Программа построения графика рис. 3.

Начните Вашу математическую модель отсюда, не забыв о ссылке и на меня. 'Бейсик, Моя программа уравнений "плюса". Уравнения Абалакина В.К..
'L1=316.20521+30.349253*T, L2=158.29638+12.220685*T,
'L3=99.08674+4.284616*T, L4=194.42558+2.184891*T, отнесённые к эклиптике, 'точке весеннего равно-денствия, стандартной эпохи 1950.0. стр. 65-66, 'Астрономический календарь. Постоянная часть. Москва. Наука. 1981.

CLS : SCREEN 12: LOCATE 10, 1: PRINT "2000": LOCATE 10, 24: PRINT "1900"
LOCATE 10, 49: PRINT "1800": LOCATE 10, 74: PRINT "1700"
S = 1950: Pi = 3.141593 / 180: LINE (20, 200)-(620, 200)
FOR x = 20 TO 620 STEP 20: LINE (x, 200 - 5)-(x, 200 + 5): NEXT x
FOR x = 20 TO 620 STEP 100: LINE (x, 200 - 10)-(x, 200 + 10): NEXT x
FOR x = 20 TO 620 STEP 200: LINE (x, 200 - 15)-(x, 200 + 15): NEXT x
FOR i = -5 TO 305 STEP .01: x1 = x2: y1 = y2: T = 1700 + i
l1 = (316.205 + 30.349253# * (T - S)) * Pi
l2 = (158.296 + 12.220685# * (T - S)) * Pi 
l3 = (99.087 + 4.284616 * (T - S)) * Pi
l4 = (194.426 + 2.184891 * (T - S)) * Pi
h1 = 10 * SIN(l1): g1 = 10 * COS(l1)
h2 = 10 * SIN(l2): g2 = 10 * COS(l2)
h3 = 10 * SIN(l3): g3 = 10 * COS(l3)
h4 = 10 * SIN(l4): g4 = 10 * COS(l4)
l12 = l1 - l2: l13 = l1 - l3: l14 = l1 - l4
l23 = l2 - l3: l24 = l2 - l4: l34 = l3 - l4
f = 20 * SIN(l12 - l34) - 10 * SIN(l12 + l34) 
IF ABS(f) < .2 THEN 1 ELSE 2
1 REM крестовые симметрии диаметров орбит
LINE (620 - i * 2 + g1, 260 - h1)-(620 - i * 2 - g1, 260 + h1)
LINE (620 - i * 2 + g2, 260 - h2)-(620 - i * 2 - g2, 260 + h2)
LINE (620 - i * 2 + g3, 260 - h3)-(620 - i * 2 - g3, 260 + h3)
LINE (620 - i * 2 + g4, 260 - h4)-(620 - i * 2 - g4, 260 + h4)
REM положение планет
CIRCLE (620 - i * 2 + g1, 260 - h1), 1
CIRCLE (620 - i * 2 + g2, 260 - h2), 1
CIRCLE (620 - i * 2 + g3, 260 - h3), 1
CIRCLE (620 - i * 2 + g4, 260 - h4), 1
2 REM не крестовые симметрии
x2 = i * 2: y2 = 200 - f
IF k = 0 THEN 5 ELSE 3
3 LINE (620 - x1, y1)-(620 - x2, y2)
5 k = k + 1: NEXT i