Dim ka1, ka2, ka3, ka4, ka1t, ka2t, ka3t, ka4t As Double
  ka1t = 0 : ka1 = 0.0
  ka2t = 15 : ka2 = 0.0
  ka3t = 60 : ka3 = 0.03
  ka4t = 240 : ka4 = 0.03
  Dim cs1, cs2, cs3, cs4, cs1t, cs2t, cs3t, cs4t As Double
  cs1t = 0 : cs1 = 0.01
  cs2t = 15 : cs2 = 0.01
  cs3t = 60 : cs3 = 0.1
  cs4t = 240 : cs4 = 0.1
Start:
  If X >= ka1t Then ka = (ka2 - ka1) / 2 * Cos(PI * (X - (ka2t - ka1t)) / (ka2t - ka1t)) + ka1 + (ka2 - ka1) / 2
  If X >= ka2t Then ka = (ka3 - ka2) / 2 * Cos(PI * (X - (ka3t - ka1t)) / (ka3t - ka2t)) + ka2 + (ka3 - ka2) / 2
  If X >= ka3t Then ka = (ka4 - ka3) / 2 * Cos(PI * (X - (ka4t - ka1t)) / (ka4t - ka3t)) + ka3 + (ka4 - ka3) / 2
  If X >= cs1t Then Cs(2) = (cs2 - cs1) / 2 * Cos(PI * (X - (cs2t - cs1t)) / (cs2t - cs1t)) + cs1 + (cs2 - cs1) / 2
  If X >= cs2t Then Cs(2) = (cs3 - cs2) / 2 * Cos(PI * (X - (cs3t - cs1t)) / (cs3t - cs2t)) + cs2 + (cs3 - cs2) / 2
  If X >= cs3t Then Cs(2) = (cs4 - cs3) / 2 * Cos(PI * (X - (cs4t - cs1t)) / (cs4t - cs3t)) + cs3 + (cs4 - cs3) / 2