Create cubicsqrt.py

Deterministic algorithm that calculates the square root of the discriminant mod p of a given cubic polynomial that is irreducible modulo p where p is a prime > 3.