Here is the fixed script
# Generator point coordinates
x=55066263022277343669578718895168534326250603453777594175500187360389116729240
y=32670510020758816978083085130507043184471273380659243275938904335757337482424
# Order
n=115792089237316195423570985008687907853269984665640564039457584007908834671663
# from the curve equation y^2=x^3+7
should_be_zero= (y**2)-(x**3+7)
print("by applying the equation, this should be equal to 0 :",should_be_zero)
# The remainder from the division by n should be 0
print("should_be_zero divided by n leaves:",should_be_zero%n)
the output
by applying the equation, this should be equal to 0 : -166977061698153803977729810299616665720111080589888563362701662779994291659332409807309461070447932090244771419528434792678509158779752908144538176572381887934774683088169260414743338484604182122883788458741320363571878334796108231
should_be_zero divided by n leaves: 0