通过软件Origin进行复杂网络博弈中b-K与r-K曲线的相关绘制（适用于学术论文）

该图表的核心在于如何识别数据点，整体构图较为直观。
在图示中，蓝色曲线上方所覆盖的区域标识为背叛者区域(D)，而红色曲线下方所涵盖的区域则表示为合作者(C)，两者之间的地带被定义为混合区域(C+D)。

以子图(a)为例，横轴表示噪声因子K（即Fermi公式中的参数），纵轴则代表背叛引诱值b（对应PDG模型）。在计算过程中，b的取值范围设定为1至1.3，步长为0.005；而K的取值区间为0至2（不包含0），步长为0.1。这意味着当K等于0.1时，需对b的不同数值进行计算（在程序运行中，当合作率降为0时，将执行break操作以节省运算时间）。总体而言，计算量大约在几百次左右（当时借助超算进行运算，预计仅需一天即可完成全部计算）。

    for(T=1.0;T<=2.01;T=T+0.005)//模型参数b
    {   
        aa=0;
        initial();//初始化函数
        for(steps=0; steps<MC_STEPS; steps++)//MCS=50000