Fixed matrix multiplication, stopped at line 219

I've gone through and debugged this code line by line.  The error I'm trying to work out now is 
"Error using griddedInterpolant/subsref
Data points in complex number format are not supported.
Use REAL and IMAG to extract the real and imaginary components.

Error in interp2 (line 161)
    Vq = F(Xq,Yq);

Error in bees (line 224)
    predictedhoney = 0.4 *
    interp2(hsurfX,hsurfY,hsurf,0.8*stage(5),stage(6)-PollenForager);
    
Error in testNOP (line 132)
		     X = bees(X,t);  % call to bees.m function, which
             outputs new state of hive"
Code compiles and runs again, but predictedhoney is a constant, winter still doesn't work

I've simplified the survivorship functions for now- there is only intrinsic survival, no relationship with honey/pollen/nursing for any of the stages.  Making this work will be the next step.

predictedhoney = 0.4 *10000; for now, there seems to be an issue with interp2 somewhere.

The winter clearly still doesn't work.
Pollen dynamics are fine, honeycollection still set to a constant
Merge pull request #2 from aipritsker/trying-to-fix-pollen-storage-issue

Trying to fix pollen storage issue
Merge pull request #3 from aipritsker/Aug27

Aug27: the code runs, but the population still looks very jagged, due to the issues with honey storage.  The amount of honey stored just sort of increases indefinitely as you increase the hive size, and there is a weird jump in the honey stored even though there is no jump in population.  The second issue is figuring out the deaths in the winter. I'm going to work on the winter issue for a while today. because I'm not seeing what's wrong with the honey.