丸め込み誤差を減少させる方法

Question

今自作の関数を作っているのですが、上記のようにおそらく丸め込み誤差によってグラフがギザギザになってしまっています。for文のrangeを(4,52,4)とすれば滑らかなグラフとなるのですが、
N_5GとN_Wを他の値に変えた時と比較したいのでfor文のrangeを変更したくありません。グラフを滑らかにするためにグラフを近似する以外の何らかの手法やサイトなどを教えて頂けると幸いです。よろしくお願い致します。

コード
import random
import math### 

nCr = {}
def comb(n, r):#コンビネーションを計算する関数
    if r == 0 or r == n: return 1
    if r == 1: return n
    if (n,r) in nCr: return nCr[(n,r)]
    nCr[(n,r)] = comb(n-1,r) + comb(n-1,r-1)
    return nCr[(n,r)]


def Judge(P1,P2):#与えられた確率に応じてユーザーが1,0どちらになるのかを判定する
    j=random.randint(0,100)%10/10
    if(P1>j):
        if(P2>j):
            return(1,1)
        else:
            return(1,0)
    else:
        if(P2>j):
            return(0,1)
        else:
            return(0,0)


def C_C(P_5G,P_W,q_C,q_W):
    Y_5G,Y_W=Judge(P_5G,P_W) 
    if(Y_5G==1):
        if(Y_W==1):
            return (1-(1-q_C)*q_W)
        else:
            return (1-(1-q_C)*(1-q_W))
    else:
        if(Y_W==1):
            return (1-q_C*q_W)
        else:
            return (1-q_C*(1-q_W))

def W_C(P_5G,P_W,q_C,q_W):
    Y_5G,Y_W=Judge(P_5G,P_W)
    if(Y_5G==1):
        if(Y_W==1):
            return q_C*(1-q_W)
        else:
            return q_W*q_C
    else:
        if(Y_W==1):
            return (1-q_C)*(1-q_W)
        else:
            return (1-q_C)*q_W

def Result_cellular(N_5G,N_W,P_5G,P_W,q_C,q_W):
    N=N_5G+N_W
    sum=0
    k=max(0,(N+1)/2-N_W)
    C=C_C(P_5G,P_W,q_C,q_W)
    W=W_C(P_5G,P_W,q_C,q_W)
    for i in range(int(k),N_5G+1):
        l=max(0,(N+1)/2-i)
        for j in range(int(l),N_W+1):
            s1=comb(N_5G,i)*comb(N_W,j)
            s2=pow(C,i)*pow((1-C),N_5G-i)
            s3=pow(W,j)*pow((1-W),N_W-j)
            sum+=s1*s2*s3
    return (sum)

def P_error(P_5G,P_W,c1,c2):
    s1=P_5G*(1-P_W)*(1-c1)  
    s2=(1-P_5G)*P_W*c2
    s3=P_5G*(1-P_W)+(1-P_5G)*P_W
    return ((s1+s2)/s3)


P_5G=0.6
P_W=0.8
P_4G=1
q_C=0.6
q_W=0.8


for i in range(3,53,2):
    N_5G=int(math.ceil(0.25*i))
    N_W=int(math.ceil(0.75*i))
    c1=Result_cellular(N_5G,N_W,1,0,q_C,q_W)
    c2=Result_cellular(N_5G,N_W,0,1,q_C,q_W)
    print(i,P_error(P_5G,P_W,c1,c2))

Answer 1

例えば gnuplot の場合、ベジエ曲線で補間することができます。

plot 'output.dat' smooth bezier