n<-12
res<-1
for(i in 0:11)
{
res<-res * (n-i)
}
print (res)
## [1] 479001600
nvec<-c(20:50)
nvc<-nvec[seq(0,length(nvec),5)]
nvc
## [1] 24 29 34 39 44 49
calc_factorial = function(a,b,c){
sol<-b^2- (4*a*c)
q<-0
if(sol<0)
{
return ("no solution")
}
else{
q<-q + sqrt(sol)
}
x1<-(b-q)/(2*a)
print(x1)
x2<-(b+q)/(2*a)
print (x2)
}
calc_factorial(1,2,1)
## [1] 1
## [1] 1