#' # Ellipse with semi-axes a=8, b=2, centered in (0,0), with orientation angle = -pi/3
C <- ellipseToConicMatrix(c(8,2),c(0,0),-pi/3)
# line
l <- c(0.25,0.85,-3)
# intersection conic C with line l:
p_Cl <- intersectConicLine(C,l)
# plot
plot(ellipse(c(8,2),c(0,0),-pi/3),type="l",asp=1)
addLine(l,col="red")
points(t(p_Cl), pch=20,col="red")
Run the code above in your browser using DataLab