quantity.Symbolic is buggy for some constants
The following code does not work:
z = quantity.Domain("z", linspace(0, 1, 4)); K1 = quantity.Symbolic.ones(1, z); K1.diff("z", 1).on()
Error using quantity.Symbolic>@()double(f)+quantity.Function.zero() Too many input arguments.