HW 1. JB

#Question 1

x <- 1
for (i in c(2:12)) {
  x <- x*i
}
print(x)

## [1] 479001600

#Question 2


nVector <- seq(20,50,5)
print(nVector)

## [1] 20 25 30 35 40 45 50

#Question 3 

quadratic <- function(a,b,c){
  if(quad(a, b, c) > 0) {
  x1 = (-b+sqrt(quad(a,b,c)))/(2*a)
  x2 = (-b-sqrt(quad(a,b,c)))/(2*a)  
  quadratic = c(x1,x2)
  print(paste(x1,x2))
  }
  else if(quad(a, b, c) == 0){
      x = -b/(2*a)  
  }
  else {"N/A - No roots OR Complex Roots"}
}

quad <- function(a, b, c){
      b^2-4*a*c  
}

quadratic(1, -2, 1)
quadratic(1, 6, 5)

## [1] "-1 -5"

quadratic(1, 1, 1)

## [1] "N/A - No roots OR Complex Roots"