t = Times; n=size(Times,1); % para = cell2mat(Wellpara); nw= 1;%Ñ¡¾®ºÅ % n = length(para); bhp = zeros(n,1); for i = 1 : n % bhp(i) = para(i).bhp(2); % bhp(i) = Wellpara{i}{1,nw}.pwf; bhp(i) = Wellpara{i}{1,nw}.qo; % bhp(i) = Wellpara{i}{1,nw}.pwf; % bhp(i) = Wellpara{i}{1,nw}.qw/(Wellpara{i}{1,nw}.qw+Wellpara{i}{1,nw}.qo); end plot(t,bhp,'b^-') hold off % ylim([0,200]); xlabel('time, day'); % ylabel('water cut'); ylabel('oil production rate, m^3/d');