import numpy as np
import matplotlib.pyplot as plt

def f(x) :
    s = 0 
    for n in range(1000) :
        if int(n**0.5) == n**0.5 :
            s += x**n
    return s

X = np.linspace(0,1,100)
Y = [(1-x)*f(x) for x in X]
plt.plot(X,Y)
plt.show()