#!/usr/bin/python import numpy as np import cmath as cm a,b,c=np.loadtxt("sv.dat") print a,b,c def root(a,b,c): dt=cm.sqrt(b**2-4.*a*c) x1=(-b+dt)/(2.*a) x2=(-b-dt)/(2.*a) return x1,x2 print root(a,b,c)