function energybanlancemodel()
%主函數
    for i=1:100 %賦初值
        u(i)=-1;
        a(i)=-1;
    end
    u(101)=0;
    a(1)=1;
    a(101)=0;
    r=0.3;
    %pi=3.14159;
    for i=100:-1:1
        u(i)=fu(u,a,i);%調用函數求各節點傳輸的數據量
        a(i)=fa(u,a,i);%調用函數求各節點采用直傳和多跳傳輸的比值
    end
    for i=1:100
        d(i)=1+u(i+1)*(1-a(i+1))/(pi*(2*i-1));%數據密度
        e(i)=r*(d(i)-1)+(a(i)*i*i-a(i)+1)*d(i);%求能量公式
        y(i)=a(i);
    end
    x=1:1:100;  
    plot(x,y);
    ylabel('\alpha');
    xlabel('半徑');
    title('\alpha Curve for Eo=1403.765');
    figure;
    y=e;
    plot(x,y);
    axis([0 120 1403.76 1403.8]);
    title('Eo=1403.765');
    a
    %d
    %u
    format long;
    e
    

%------------------------------function fa(u,a,i)--------------------------
function ra=fa(u,a,i)
%求直傳和多跳傳輸的比值
    if a(i)~=-1
        ra=a(i);return;
    end
    e=1403.765;r=0.3;
    %以下是a(x)和u(x)的關系式
    a(i)=(e*pi*(2*i-1)-r*fu(u,a,i+1)*(1-fa(u,a,i+1))-(pi*(2*i-1)+fu(u,a,i+1)*(1-fa(u,a,i+1))))/((pi*(2*i-1)+fu(u,a,i+1)*(1-fa(u,a,i+1)))*(i*i-1));
    ra=a(i);
    
%----------------------------------function fu(u,a,i)--------------
function ru=fu(u,a,i)
%求傳輸的數據量
    if u(i)~=-1 
        ru=u(i);return;
    end
    %pi=3.14159;
    u(i)=pi*(2*i-1)+(1-fa(u,a,i+1))*fu(u,a,i+1);
    ru=u(i);