Add comments; copy in blueprint from mean cache functions
Call solve instead of mean; rename variables and functions