input: Write a python function to calculate the relative differences between fluxes computed by Wilson and Horvat irradiation methods as a function of phase for given inclinations.
  assert calculate_flux_difference([0, 30, 60, 90]) == 'Expected output based on the provided code logic'
  assert calculate_flux_difference([10, 45, 75]) == 'Expected output based on the provided code logic'
  assert calculate_flux_difference([15, 35, 55, 85]) == 'Expected output based on the provided code logic'
output: def calculate_flux_difference(inclinations):
      def wilson_flux(phase):
          return 1 + 0.1 * phase
  
      def horvat_flux(phase):
          return 1 + 0.15 * phase
  
      def relative_difference(flux1, flux2):
          return (flux2 - flux1) / flux1
  
      phases = [0, 0.25, 0.5, 0.75, 1]
      differences = {}
  
      for inclination in inclinations:
          differences[inclination] = []
          for phase in phases:
              wf = wilson_flux(phase)
              hf = horvat_flux(phase)
              diff = relative_difference(wf, hf)
              differences[inclination].append(diff)
      
      return differences