findequationshigh.m

电网机制设计

clc

global Tlimit

% guess the initialization
t=cputime;
firstvalue = 75;
ntimes = 150
t0=linspace(1, firstvalue, ntimes);

arr = zeros(1,ntimes);
for i=1:ntimes    
    arr(i)=equationshigh(t0(i),Icost);
    disp('---Be patient, it is computing!---')
end
% (1-lambda)+lambda.*sqrt(pi)./2-((1-lambda).*alpha+1)./(alpha.*capa)
index=find(arr==min(arr))
Tlimit=t0(index)

ComputingTime=cputime-t