function[R1,R2]=secroot(a,b,c) Det=b^2-4*a*c if (Det<0) R1=(-b+j*sqrt(-Det))/2/a; R2=(-b-j*sqrt(-Det))/2/a; elseif(Det==0) R1=-b/2/a; R2=-b/2/a; else(Det>0) R1=(-b+sqrt(Det))/2/a; R2=(-b-sqrt(Det))/2/a; end