सूर्यसिद्धान्तको गणितीय रहस्य: संस्कृत श्लोक, अहर्गण सूत्र र शून्य-लुकअप (Zero-Lookup) खगोल इन्जिनियरिङ

२०२६ अगस्ट २७ By khumnath

परिचय

सूर्यसिद्धान्त (Surya Siddhanta) मानव इतिहासकै सबैभन्दा प्राचीन र गणितीय दृष्टिले समृद्ध खगोलीय ग्रन्थहरूमध्ये एक हो। वैदिक ज्योतिष तथा खगोल परम्परामा आधारित यस ग्रन्थले लाखौं वर्षसम्मको ग्रहस्थिति, तिथि, नक्षत्र, सूर्य तथा चन्द्रग्रहण, ग्रहयुद्ध र ऋतु चक्रहरूको सटीक भविष्यवाणी गर्न सक्ने एउटा पूर्ण खगोलीय गणितीय प्रणाली प्रस्तुत गरेको छ।

सूर्यसिद्धान्तको प्रथम अध्यायको प्रथम श्लोक, १८४७ संस्कृत पाण्डुलिपि संस्करण

प्रायः आधुनिक पात्रो तथा क्यालेन्डर एपहरूले निश्चित वर्षका तिथि र चाडपर्व देखाउन पहिले नै हिसाब गरेर राखिएका स्थिर तालिकाहरू (Static Lookup Tables) प्रयोग गर्छन्। तर यस्ता हार्डकोडेड तालिकाहरूले खगोल विज्ञानको वास्तविक भौतिक नियमलाई दर्शाउन सक्दैनन् र ती हजारौं वर्षको समयमा काम गर्दैनन्।

सूर्यसिद्धान्तको वास्तविक गणितलाई आधुनिक कम्प्युटर प्रोग्रामिङमा उतार्न मैले Surya Siddhanta Audit Project निर्माण गरेको छु—जसले कुनै पनि पूर्वनिर्धारित तालिका बिना (Zero-Lookup Architecture) शुद्ध सूत्रबाट सम्पूर्ण खगोलीय मानहरू निकाल्दछ।

🌐 प्रत्यक्ष खगोलीय क्यालकुलेटर र अडिट ड्यासबोर्ड: khumnath.com.np/surya-siddhanta-audit
📖 सम्पूर्ण API तथा शास्त्रीय श्लोक प्रमाणहरू: khumnath.com.np/surya-siddhanta-audit/api/index.html

यस लेखमा हामी सूर्यसिद्धान्तका मुख्य गणितीय सिद्धान्तहरू, मूल संस्कृत श्लोकहरू र तिनीहरूको आधुनिक खगोल विज्ञानसँगको तुलनात्मक विश्लेषण गर्नेछौं।

१. काल गणनाका नौ प्रकार (नव मानानि)

सूर्यसिद्धान्तको चौधौं अध्याय (मानाध्याय) मा समय नाप्ने नौवटा आधारभूत प्रणालीहरूको वर्णन गरिएको छ:

+-------------------------------------------------------------------------+
|                  शास्त्रीय प्रमाण: मानाध्याय (अध्याय १४, श्लोक १-२)       |
+-------------------------------------------------------------------------+
| संस्कृत (देवनागरी):                                                     |
| ब्राह्मं दिव्यं तथा पित्र्यं प्राजापत्यं च गौरवम्‌।                     |
| सौरञ्च सावनं चान्द्रमाक्षं मानानि वै नव।।                                |
|                                                                         |
| नेपाली अनुवाद:                                                          |
| "ब्राह्म, दिव्य, पित्र्य, प्राजापत्य, गौरव (बृहस्पति), सौर, सावन,       |
| चान्द्र र नाक्षत्र—यी नौवटा काल गणनाका प्रमाण (मान) हुन्।"             |
+-------------------------------------------------------------------------+

                        ┌───────────────────────────────┐
                        |      कालका ९ मान (Systems)     |
                        └───────────────┬───────────────┘
                                        │
        ┌───────────────┬───────────────┼───────────────┬───────────────┐
        ▼               ▼               ▼               ▼               ▼
 ┌─────────────┐ ┌─────────────┐ ┌─────────────┐ ┌─────────────┐ ┌─────────────┐
 │   सौर मान   │ │  चान्द्र मान │ │  सावन मान   │ │  नाक्षत्र मान│ │   गौरव मान  │
 │ (सौर वर्ष)  │ │ (चान्द्र मास)│ │ (सावन दिन)  │ │ (नाक्षत्र दिन│ │(बृहस्पतिवर्ष│
 └─────────────┘ └─────────────┘ └─────────────┘ └─────────────┘ └─────────────┘

यी नौ प्रकारका मानहरू विभिन्न व्यावहारिक र खगोलीय प्रयोजनमा प्रयोग हुन्छन्:

सौर मान (Solar Time): सूर्यले १२ राशि पार गर्ने समय। यसबाट १२ सौर महिना, ६ ऋतु र १ सौर वर्ष बन्दछ।
चान्द्र मान (Lunar Time): चन्द्रमाको कला (तिथि तथा शुक्ल/कृष्ण पक्ष)। धार्मिक चाडपर्व, व्रत र अधिकमास (मलमास) गणना गर्न यो मान चाहिन्छ।
सावन मान (Civil Terrestrial Time): एक सूर्योदयदेखि अर्को सूर्योदयसम्मको समय (अहोरात्र)। दैनिक व्यवहार, बार (Vara), र अहर्गण गणना गर्न सावन दिन प्रयोग हुन्छ।
नाक्षत्र मान (Sidereal Time): स्थिर ताराहरूको सापेक्षमा पृथ्वीको एक परिक्रमा (६० घडी = करिब २३ घण्टा ५६ मिनेट ४ सेकेन्ड)। ग्रहहरूको मध्यम गति निकाल्न यो प्रयोग हुन्छ।
गौरव / बार्हस्पत्य मान (Jovian Time): बृहस्पति ग्रहले एक राशि पार गर्ने समय (करिब १ वर्ष), जसबाट ६० संवत्सरको चक्र बन्दछ।
दिव्य, पित्र्य, प्राजापत्य, र ब्राह्म मान: सृष्टिको उत्पत्ति र भूगर्भीय कालखण्ड मापन गर्ने महा-एकाइहरू (जस्तै: १ दिव्य वर्ष = ३६० सौर वर्ष; १ कल्प = ४ अर्ब ३२ करोड वर्ष = ब्रह्माको १ दिन)।

२. महायुगको संरचना र भगण (Revolutions)

सूर्यसिद्धान्त अनुसार ग्रहहरूको गतिलाई एक महायुग (४३,२०,००० सौर वर्ष) मा पूरा हुने परिक्रमा सङ्ख्या (भगण) का रूपमा परिभाषित गरिएको छ:

+-------------------------------------------------------------------------+
|             शास्त्रीय प्रमाण: मध्यमाधिकार (अध्याय १, श्लोक २९)            |
+-------------------------------------------------------------------------+
| संस्कृत (देवनागरी):                                                     |
| युगे सूर्यबुधशुक्राणां खचतुष्करदार्णवाः।                                |
| कुजार्किगुरुशीघ्राणां भगणाः पूर्वयायिनाम्॥                              |
|                                                                         |
| नेपाली अनुवाद:                                                          |
| "एक महायुगमा सूर्य, बुध र शुक्रको भगण सङ्ख्या ४३,२०,००० हुन्छ।          |
| मङ्गल, शनि र बृहस्पतिको शीघ्र विन्दुको भगण पनि यति नै हुन्छ।"          |
+-------------------------------------------------------------------------+

२.१ शब्द-सङ्ख्या सङ्केत (भूतसङ्ख्या पद्धति)

श्लोकमा ४३,२०,००० सङ्ख्यालाई ख-चतुष्-करद-अर्णवाः भनिएको छ:

ख = शून्य (००)
चतुष् = चार (४)
रद = दाँत = ३२
अर्णव = समुद्र = ४
अङ्कको वाम गति नियम (Ankasya Vama Gati) अनुसार उल्टो क्रममा पढ्दा: ४३,२०,०००।

२.२ महायुगमा विभिन्न पिण्डहरूको भगण तालिका

पिण्ड / विन्दु	संस्कृत नाम	महायुग भगण (R_i)
सूर्य (Sun)	सूर्य	४३,२०,०००
चन्द्रमा (Moon)	चन्द्र	५,७७,५३,३३६
चन्द्र मन्दोच्च (Apogee)	मन्दोच्च	४,८८,२०३
राहु (Lunar Node)	पातः / राहु	-२,३२,२३८ (वक्री)
बुध शीघ्र (Mercury Sighra)	बुध-शीघ्र	१,७९,३७,०६०
शुक्र शीघ्र (Venus Sighra)	शुक्र-शीघ्र	७०,२२,३७६
मङ्गल (Mars)	भौम / कुज	२२,९६,८३२
बृहस्पति (Jupiter)	बृहस्पति / गुरु	३,६४,२२०
शनि (Saturn)	शनि / मन्द	१,४६,५६८
सावन दिन (Civil Days in Yuga)	सावन दिन	१,५८,२२,३७,८२८

यसबाट १ सौर वर्षको लम्बाइ शुद्ध गणितबाट निस्कन्छ:

१ नाक्षत्र वर्ष = १,५८,२२,३७,८२८ / ४३,२०,००० = ३६५.२५८७५६४८ सावन दिन
               = ३६५ दिन, ६ घण्टा, १२ मिनेट, ३६.५६ सेकेन्ड

३. अहर्गण (Ahargana) गणना विधि

शून्य-लुकअप खगोल विज्ञानको मेरुदण्ड अहर्गण (दिनहरूको थुप्रो) हो। कलियुगको सुरुवाती क्षण (१८ फेब्रुअरी ३१०२ ई.पू. मध्यरात) देखि आजसम्म बितेका कुल सावन दिनहरूको सङ्ख्यालाई अहर्गण भनिन्छ।

+-------------------------------------------------------------------------+
|             शास्त्रीय प्रमाण: मध्यमाधिकार (अध्याय १, श्लोक ४५-५१)         |
+-------------------------------------------------------------------------+
| संस्कृत (देवनागरी):                                                     |
| सौरमासाः कृतगुणा युगेऽधिमासकैर्युताः।                                  |
| भवन्ति चान्द्रमासास्ते तिथिभिर्गुणिताः पृथक्॥                            |
| अधिमासावमोनास्ताः सावनं द्युगणं विदुः॥                                   |
+-------------------------------------------------------------------------+

३.१ अहर्गण निकाल्ने पाँच चरणहरू:

बितेका सौर महिना (M_solar):

M_solar = (कलियुग वर्ष × १२) + चालू सौर महिना

अधिकमास (मलमास - M_adhika): एक महायुगमा सौर महिनाभन्दा बढी हुने चान्द्र महिना:
```
अधिकमास = floor( (M_solar × १,५९३,३३६) / ५१,८४०,००० )
```

बितेका तिथिहरू (T_elapsed):

चान्द्र महिना = M_solar + M_adhika
T_elapsed = (चान्द्र महिना × ३०) + चालू तिथि

अवम दिन (क्षय तिथि - D_avama):

D_avama = floor( (T_elapsed × १५०,३६२,२५२) / १,७३२,६००,०८० )

अन्तिम अहर्गण (Ahargana: A):
```
अहर्गण (A) = T_elapsed - D_avama
```

अहर्गण $A$ प्राप्त भएपछि कुनै पनि ग्रहको मध्यम भोगांश (λ̄) सिधै निस्कन्छ:

मध्यम भोगांश (λ̄_i) = ( (A × R_i) / १,५८,२२,३७,८२८ × ३६०° ) mod ३६०°

४. मन्द र शीघ्र फल (Equation of Center & Parallax)

मध्यम ग्रह केवल गोलाकार औसत गति हो। आकाशमा देखिने वास्तविक स्पष्ट ग्रह निकाल्न सूर्यसिद्धान्तले दुईवटा वृत्त (Epicycles) को प्रयोग गर्दछ:

                      मध्यम ग्रह (Mean Planet: M)
                                 │
                 ┌───────────────┴───────────────┐
                 ▼                               ▼
        मन्द संस्कार (कक्षीय अण्डाकारता)   शीघ्र संस्कार (दृष्टि भ्रम/Parallax)
                 │                               │
            Equation of Center              Heliocentric to Geocentric
                 ▼                               ▼
        मन्द स्पष्ट ग्रह (M')              अन्तिम स्पष्ट ग्रह (λ_true)

+-------------------------------------------------------------------------+
|             शास्त्रीय प्रमाण: स्पष्टाधिकार (अध्याय २, श्लोक १-४, ४३)      |
+-------------------------------------------------------------------------+
| संस्कृत (देवनागरी):                                                     |
| अदृश्यरूपाः कालस्य मूर्तयो भगणाश्रिताः।                                 |
| शीघ्रमन्दोच्चपाताख्या ग्रहाणां गतिहेतवः॥                                |
| तद्वातवशगाद् बद्धाः प्रागपक्कृष्टिमूर्तयः।                            |
| खेचरा विचरन्त्येवं स्वमार्गेण यथायथम्॥                                  |
+-------------------------------------------------------------------------+

मन्द फल (Manda Phala - कक्षीय अण्डाकारता):

मन्द केन्द्र (θ_manda) = मध्यम ग्रह - मन्दोच्च
sin(मन्द फल μ) = (मन्द परिधि / ३६०°) × sin(θ_manda)
मन्द स्पष्ट ग्रह = मध्यम ग्रह - μ

शीघ्र फल (Sighra Phala - दृष्टि भ्रम र पृथ्वीको गति):

शीघ्र केन्द्र (θ_sighra) = शीघ्र विन्दु - मन्द स्पष्ट ग्रह
tan(शीघ्र फल σ) = (r × sin(θ_sighra)) / (R + r × cos(θ_sighra))
अन्तिम स्पष्ट ग्रह = मन्द स्पष्ट ग्रह + σ

५. पञ्चाङ्गका पाँच अङ्गको प्रत्यक्ष कोड गणना (Zero-Lookup TypeScript)

अहर्गणबाट सूर्य र चन्द्रमाको स्पष्ट भोगांश निस्किएपछि पञ्चाङ्गका पाँच अङ्गहरू कुनै तालिका बिना सिधै हिसाब हुन्छन्:

// सूर्यसिद्धान्त अडिट इन्जिनको वास्तविक TypeScript कोड उदाहरण

export function calculateTithi(sunLong: number, moonLong: number): TithiInfo {
  // चन्द्रमा र सूर्यको कोणीय दूरी (Elongation)
  const elongation = normalizeAngle(moonLong - sunLong);
  
  // प्रत्येक तिथि बराबर १२ डिग्री (७२० कला)
  const tithiFraction = elongation / 12.0;
  const tithiIndex = Math.floor(tithiFraction) + 1;
  const remainingFraction = 1.0 - (tithiFraction - Math.floor(tithiFraction));

  return {
    index: tithiIndex,
    name: TITHI_NAMES[tithiIndex - 1],
    paksha: tithiIndex <= 15 ? 'Shukla' : 'Krishna',
    pakshaTithiIndex: tithiIndex <= 15 ? tithiIndex : tithiIndex - 15,
    remainingFraction
  };
}

export function calculateNakshatra(moonLong: number): NakshatraInfo {
  // प्रत्येक नक्षत्र बराबर १३ डिग्री २० मिनेट (८०० कला)
  const nakshatraSpan = 13.333333333333334;
  const nakshatraFraction = normalizeAngle(moonLong) / nakshatraSpan;
  const nakshatraIndex = Math.floor(nakshatraFraction) + 1;
  const pada = Math.floor((nakshatraFraction - Math.floor(nakshatraFraction)) * 4) + 1;

  return {
    index: nakshatraIndex,
    name: NAKSHATRA_NAMES[nakshatraIndex - 1],
    pada,
    ruler: NAKSHATRA_RULERS[nakshatraIndex - 1]
  };
}

६. आधुनिक खगोल विज्ञानका आवश्यक प्रतिस्थापन सूत्रहरू (Modern Theoretical Replacements)

आधुनिक दृक (Drigganita) खगोलशास्त्रले कसरी अत्यन्त सूक्ष्म स्तरको सटिकता प्राप्त गर्छ र प्राचीन सूर्यसिद्धान्तका कुन-कुन सूत्रलाई आधुनिक भौतिकशास्त्रले विस्थापन गरेको छ भन्ने कुरा तलका मुख्य सूत्रहरूबाट स्पष्ट हुन्छ:

+-------------------------------------------------------------------------+
|        शास्त्रीय अनुमान र आधुनिक खगोलीय भौतिकशास्त्र बीचको तुलना         |
+------------------------------+------------------------------------------+
| शास्त्रीय सिद्धान्त विधि      | आधुनिक वैज्ञानिक प्रतिस्थापन             |
+------------------------------+------------------------------------------+
| • स्थिर मन्द परिधि (Epicycle)| • केपलरको समीकरण र वास्तविक कोण (ν)      |
| • चन्द्रमाको एउटै मन्द वृत्त | • बहु-पिण्ड चन्द्र सिद्धान्त (ELP-2000)  |
|   (अन्य आकर्षण समावेश नभएको) |   (Evection, Variation, Annual Equation) |
| • सरल अहर्गण अनुपात          | • VSOP87 गुरुत्वाकर्षण विचलन शृङ्खला     |
| • स्थिर भूकेन्द्रीय ढाँचा     | • न्युटेसन (Δψ) र स्टेलर एबरेसन (κ)      |
| • स्थिर शून्य विन्दु         | • IAU / लाहिरी उच्च सटिक अयनचलन (Ayanamsa)|
+------------------------------+------------------------------------------+

६.१ केपलरको समीकरण र अण्डाकार कक्ष (मन्द फलको आधुनिक प्रतिस्थापन)

सूर्यसिद्धान्तको मन्द वृत्त कक्षीय अण्डाकारताको सरल पहिलो श्रेणीको अनुमान मात्र हो। आधुनिक खगोलशास्त्रमा ग्रहको वास्तविक स्थिति निकाल्न केपलरको ट्रान्सिडेन्टल समीकरण (Kepler’s Equation) प्रयोग गरिन्छ:

मध्यम कोणबाट उत्केन्द्र कोण (Eccentric Anomaly: E):
```
M = E - e × sin(E)
```
यसलाई Newton-Raphson विधि द्वारा कम्प्युटरमा निरन्तर दोहोर्याएर निकालिन्छ:
```
E_{k+1} = E_k - ( (E_k - e × sin(E_k) - M) / (1 - e × cos(E_k)) )
```
उत्केन्द्र कोणबाट वास्तविक स्पष्ट कोण (True Anomaly: ν):
```
tan(ν / 2) = sqrt( (1 + e) / (1 - e) ) × tan(E / 2)
```
अथवा विश्लेषणात्मक समीकरण विस्तार (Equation of Center Expansion):
```
केन्द्र समीकरण (ν - M) = (2e - (1/4)e³) × sin(M) + (5/4)e² × sin(2M) + (13/12)e³ × sin(3M) + ...
```

६.२ चन्द्रमाका प्रमुख गुरुत्वाकर्षण विचलनहरू (Lunar Perturbations)

चन्द्रमाले पृथ्वीको परिक्रमा गर्दा सूर्यको तीव्र गुरुत्वाकर्षणले गर्दा चन्द्रमाको कक्ष स्थिर रहँदैन। आधुनिक चन्द्र सिद्धान्त (ELP-2000/82) ले चन्द्रमाको स्पष्ट भोगांश निम्न प्रमुख कम्पनात्मक पदहरू (Harmonic Terms) को जोडबाट निकाल्दछ:

λ_चन्द्रमा = L_मध्यम 
           + ६.२८८८° × sin(M')         [१. मुख्य केन्द्र समीकरण (e = ०.०५४९)]
           + १.२७४०° × sin(२D - M')    [२. एवेक्सन (Evection): सौर्य आकर्षणले चन्द्र कक्षा तन्किने]
           + ०.६५८३° × sin(२D)         [३. रिएसन (Variation): चन्द्र कला अनुसारको गति परिवर्तन]
           - ०.१८५६° × sin(M)          [४. वार्षिक समीकरण (Annual Equation): पृथ्वी-सूर्य दूरी अन्तर]
           - ०.११४३° × sin(२F)         [५. क्रान्तिवृत्त रूपान्तरण: ५.१४° को कक्षीय झुकाव]
           + (सयौं साना सूक्ष्म गुरुत्वाकर्षण प्रभावहरू...)

जहाँ कोणहरू:

M’: चन्द्रमाको मध्यम केन्द्र (Moon’s Mean Anomaly)
M: सूर्यको मध्यम केन्द्र (Sun’s Mean Anomaly)
D: सूर्य र चन्द्रमाको कोणीय दूरी (Elongation: λ_{चन्द्रमा} - λ_{सूर्य})
F: चन्द्रमाको अयन विन्दुबाट कोणीय दूरी (Argument of Latitude)

६.३ सौर्यमण्डलका ग्रहहरूको आपसी खिचावट (VSOP87 Series)

ग्रहहरूको कक्ष स्थिर नभई एक-अर्काको गुरुत्वाकर्षणले गर्दा दोहोरो कम्पनमा हुन्छन् (जस्तै: बृहस्पति र शनिको ५:२ को कक्षीय अनुनाद)। आधुनिक खगोलशास्त्रले यसलाई पोइसन शृङ्खला (Poisson Series) द्वारा हिसाब गर्छ:

भोगांश L(t) = ∑_{k=0}^{n} t^k × [ ∑_j A_{j,k} × cos(B_{j,k} + C_{j,k} × t) ]

६.४ पृथ्वीको न्युटेसन (Nutation) र स्टेलर एबरेसन (Aberration)

पृथ्वीको घुमिरहेको अक्ष र प्रकाशको निश्चित गतिको प्रभावलाई हटाउन आधुनिक गणनामा निम्न संशोधन गरिन्छ:

न्युटेसन (Δψ - अक्षको सूक्ष्म कम्पन): चन्द्रमाको १८.६ वर्षे राहु चक्रले गर्दा पृथ्वीको अक्ष हल्लिने प्रभाव:
```
Δψ = -१७.२०" × sin(Ω) - १.३२" × sin(२L_सूर्य) - ०.२३" × sin(२L_चन्द्र) + ०.२१" × sin(२Ω)
```
स्टेलर एबरेसन (κ - प्रकाशको गति र पृथ्वीको वेग):
```
Δλ_aberration = -२०.४९५५" × cos(λ_स्पष्ट - λ_सूर्य) / cos(β)
```

६.५ आधुनिक निरयण रूपान्तरण (Chitrapaksha / Lahiri Ayanamsa)

आधुनिक सायन (Tropical) मानलाई वैदिक निरयण (Sidereal) मानमा बदल्न उच्च सटिक अयनचलन घटाइन्छ:

λ_निरयण = ( λ_सायन (Apparent Tropical) - अयनचलन (Ayanamsa) ) mod ३६०°

जहाँ J2000.0 मानक समयबाट जुलियन शताब्दी ($T$) अनुसार लाहिरी अयनचलन:

Ayanamsa_Lahiri(T) = २३° ५१' २५.५३" + ५०२९.०९६६" × T + १.११२०" × T² + ...

निष्कर्ष

सूर्यसिद्धान्त केवल पुरानो ग्रन्थ मात्र नभएर गणित र खगोलशास्त्रको एउटा अद्वितीय वैज्ञानिक धरोहर हो। यसलाई शून्य-लुकअप TypeScript इन्जिनका रूपमा पुनर्जीवित गरेर हामीले प्राचीन वैदिक खगोलशास्त्रको सम्मान गर्दै आधुनिक अनुसन्धानकर्ताहरूका लागि एउटा खुला र पारदर्शी सफ्टवेयर तयार गरेका छौं।

थप अध्ययन र प्रत्यक्ष गणनाका लागि:

🌐 सूर्यसिद्धान्त अडिट ड्यासबोर्ड: khumnath.com.np/surya-siddhanta-audit
📖 API र शास्त्रीय श्लोक प्रमाणहरू: khumnath.com.np/surya-siddhanta-audit/api/index.html
💻 नेपडेट C++ कोर इन्जिन: github.com/khumnath/nepdate

सन्दर्भ सामग्रीहरू

Surya-Siddhanta: A Text-Book of Hindu Astronomy, रेभ. एबनेजर बर्गेस (Ebenezer Burgess) अनुवाद, फणीन्द्रलाल गाङ्गुली सम्पादन (१८६० / १९३५)।
सूर्यसिद्धान्त, पण्डित बापूदेव शास्त्री (Bapudeva Sastri) विरचित मूल संस्कृत तथा व्याख्या (१८६१)।
केतकर, व्यङ्कटेश बापूजी (बापू केतकर - V. B. Ketkar). केतकी ग्रहगणितम् (Ketaki Grahaganitam) एवं ज्योतिर्गणित (१८९० / १९३०)।
Cross, Donald. Astronomy Engine (NASA JPL भेक्टर गणना), २०१८–२०२४।
Rao, S. Balachandra. Indian Astronomy: An Introduction, Universities Press (२०००)।