s _ BM( 1)= 100； Hairless percentage

s _ Cao( 1)= 0； % grassland

s _ GM( 1)= 0； % shrub

s _ HW( 1)= 0； Percentage of poverty

s _ GM(2)= 0； s _ HW(2)= 0；

S _ Cao(2)= S _ Cao( 1)+0.3 * S _ BM( 1)；

S _ BM(2)= S _ BM( 1)-S _ Cao(2)-S _ GM(2)-S _ HW(2)；

Because i=3:500

S _ GM(I)= S _ GM(I- 1)+0. 1 * S _ Cao(I- 1)-0.04 * S _ GM(I- 1)；

S _ Cao(I)= S _ Cao(I- 1)+0.3 * S _ BM(I- 1)-0. 1 * S _ Cao(I- 1)-0.05 * S _ Cao(I- 1)+0.04 * S _ GM(I- 1)；

S _ HW(I)= S _ HW(I- 1)+0.05 * S _ Cao(I- 1)；

S _ BM(I)= 100-S _ GM(I)-S _ Cao(I)-S _ HW(I)；

end

t = 1: 1:500；

plot(t，S_bm，' *r '，t，S_cao，' og '，t，S_gm，'。 b '，t，S_hw，' & gtc’)