code: func calcObliqCorr(meanObliqEcliptic []float64, julianCentury []float64) (obliqCorr []float64) {
  	if len(meanObliqEcliptic) != len(julianCentury) {
  		return
  	}
  
  	for index := 0; index < len(julianCentury); index++ {
  		temp := meanObliqEcliptic[index] + 0.00256*math.Cos(deg2rad(125.04-1934.136*julianCentury[index]))
  		obliqCorr = append(obliqCorr, temp)
  	}
  	return
  }
docs: // Calculate the Obliq Corr in degrees based on the formula: meanObliqEcliptic+0.00256*cos(deg2rad(125.04-1934.136*julianCentury))
  // meanObliqEcliptic - Mean Obliq Ecliptic calculated by the calcMeanObliqEcliptic function
  // julianCentury - Julian century calculated by the calcJulianCentury function
  // Return the Obliq Corr slice