function [PL1]=main1(D,P,Xl)PL1=zeros(length(D),1);PLL=D(:,5);area1=[1 2 5 7 28 3 4 6 8 9 11];area2=[13 12 16 14 18 17 15 19 20 23];area3=[10 27 22 21 30 24 25 29 26];%例如計算區域1到區域3的ATC則增加區域1發電，增加區域3負荷；while sum(abs(PL1)>PLL)==0&sum(P(:,2)>P(:,4))    P(area1,2)=P(area1,2)+(10/sum(P(area1,2)))* P(area1,2);%增加區域1發電;    P(area3,3)=P(area3,3)+(10/sum(P(area3,3)))* P(area3,3);%增加區域3負荷；    [O,PL]=chaoliu(D,P,Xl);    PL1=PL;end