可這樣：

x=[220.8,799.1,559.9,510.4,786.9,404.4,477.2,461.2,333.1,345.6,137.9,367,227.6,160.1,144.3,260.1 ];

y=[4798,14975,9134,6426,24426,6829,7320,8321,8311,7809,3883,9365,5164,3700,2926,5047 ];

x1=cumsum(x)./sum(x)*100;

y1=cumsum(y)./sum(y)*100;

plot(x1,y1)

hold on

plot([0,100],[0,100])

hold off

xlabel("人口數序列 % ")

ylabel("護士數 % ")

title("衛生資源分配洛倫茲曲線")

legend("實際分配曲線","理論公平曲線","Location","NW")

grid on

B=trapz(x1,y1);

A=trapz(x1,x1-y1);

G=A/(A+B)%吉尼係數