epley(w,r) = w*(1+r/30)

brzycki(w,r)=w*36/(37-r)

mcglothin(w,r) = 100*w/(101.3 - 2.67123*r)

lombardi(w,r) = w * r^0.10

mayhew(w,r) = 100*w / (52.2 + 41.9 * exp(-0.055*r))

oconner(w,r) = w*(1+r/40)

wathan(w,r) = 100*w / (48.8 + 53.8 * exp(-0.075*r))

wendler(w,r) = w*r*0.0333 + w

maxFunctions = (epley, brzycki, mcglothin, lombardi, mayhew, oconner, wathan, wendler)

max(w, r) = collect(map(x -> [x; x(w,r)], maxFunctions))