complete linear systems setup and solving but does not work now
The total protein mass blows up when solving these equations iteratively