for now only full substitution works; partial substitution needs more work: e.g. f(a = 5.0, w = 2.5) should generate a new formula