factorial=prod(1:12)
factorial
## [1] 479001600
vect=seq(20,50,by=5)
vect
## [1] 20 25 30 35 40 45 50
quadratic=function(a,b,c){
  x1=-(b+sqrt(b^2-4*a*c))/(2*a)
  x2=-(b-sqrt(b^2-4*a*c))/(2*a)
  return(c(x1,x2))
}

quadratic(1,-2,1)
## [1] 1 1