Federal Tax in Calca

fedtaxes.calca

brackets = [ 13,150, 50,200, 129,601, 209,851, 411,501, 439,000, 999,999,999,999]
rates =    [    10%,    15%,     25%,     28%,     33%,     35%,           39.6%]

tax (amt, level) =
  let this = brackets[level] in
  let prev = brackets[level-1] in
  let rate = rates[level] in
  if this <= amt     then (this - prev) * rate
  else if prev < amt then (amt - prev) * rate
  else 0

mytax(income) =
  sum(tax(income, l), 0..6)
  
effective tax rate (i) =
  mytax(i) / i * 100

plot(mytax(income), income, 10,000..250,000)
plot(effective tax rate, 10,000..250,000)