getwd()
source("plot_im.R") 
source("plot_im.R") 
source("plot_im.R") 
z <- (0+1i)^(1/4) * (1, 1.05, 0.95)
z <- (0+1i)^(1/4) * c(1, 1.05, 0.95)
z
source("plot_im.R") 
a=1; b=-1; c=3
discr = b^2 - 4*a*c
source("eq2_0.R")
discr
source("eq2_0.R")
source("eq2_0.R")
source("eq2_0.R")
source("eq2_0.R")
source("eq2_1.R")
a=1; b=-1; c=3
discr = b^2 - 4*a*c
print(discr)
if(discr <0) discr <- discr * (1+0i)
discr
source("eq2_2.R")
a=1; b=-1; c=-3
plot(function(x) a*x^2+b*x+c, -3, 3)   
abline(h=0)
source("eq2_1.R")
sol
abline(v=sol[1])
abline(v=sol[2])
q()