x<-c(4.775591,4.129769,10.820459,2.926689,8.806668,4.783637,8.418380,7.833450,4.537405,6.730402,3.698800,6.864923,3.225586,5.634067,5.504755,6.331493,3.726112,6.496017,3.682233,6.843869,4.193416,6.293272,7.568595,5.758989, 8.720617,3.732230,4.492452,3.948485,7.427152,4.360651, 5.411284,1.328836,6.282210,7.916329,6.750422,5.529266, 3.653097,7.401630,3.742837,8.358330)
xprom=mean(x)
xsd<-sd(x)
n<-length(x)
print(xprom)
## [1] 5.71601
print(xsd)
## [1] 1.999193
mu=40.0
tcalc<-(xprom-mu)/(xsd/sqrt(n))
print(tcalc)
## [1] -108.4593
var(x)
## [1] 3.996771
median(x)
## [1] 5.581666
hist(x)

pt(-108.4593,n-1)
## [1] 2.664416e-50
print
## function (x, ...)
## UseMethod("print")
## <bytecode: 0x0000000013896c18>
## <environment: namespace:base>