To take into acount the exponential nature of your data, you can use the loglog scale as follows:
sage: points(zip(x,expected), color='blue', marker='o', scale='loglog') + points(zip(x,simulated), color='red', marker='x', scale='loglog')

sage: error = [abs(a-b)/a*100 for a,b in zip(expected,simulated)]