% plot cost as a funciton of constraints 

load('TCcol_min18_min10_RC_03_07_5point')



%price_vector=[D, B, boilup_energy, comp_salesgas, comp_boster];

price_vector=[-4,-2,1,1,1]
cost_variables=cost(2:end,8:12)'
costfunction=price_vector*cost_variables;

[TC_feed,RC]=meshgrid(linspace(-18,-10,5),linspace(0.3,0.7,5));

cost_matrix=zeros(5,5);
k=1;
  for i=1:5
        for j=1:5

            cost_matrix(i,j)=costfunction(k);	   
            k=k+1;
        end
  end
    

figure(2)
surf(TC_feed,RC,cost_matrix)
xlabel('Temperature feed sep')
ylabel('Split liquid cold separator')
zlabel('Cost')  
