pBin<-function(k,p,n)
{
newton<-factorial(n)/(factorial(k)*factorial(n-k))
return (newton*p^(k)*(1-p)^(n-k));


}
pBin(1,0.5,10);

normal<-function(left,right)
{
x<-seq(from=min(-3,left),to=max(3,right),length=100);
y<-dnorm(x)
plot(main="Rozkład normalny",x,y,type="l",lwd=1,col="red",xlab="Wartość",ylab="gęstość");
lines(x[x>left&x<right],y[x>left&x<right],type="h",col="green",lwd=5);
lines(x,y,lwd=2);
lines(right,dnorm(right),type="h",lwd=1);
lines(left,dnorm(left),type="h",lwd=1);
lines(right,dnorm(right),type="p",lwd=1);
lines(left,dnorm(left),type="p",lwd=1);
legend("topright",1,c(paste("Phi(",left,") = ",round(pnorm(left),3)),paste("Phi(",right,") = ",round(pnorm(right),3))));

}
normal(-3,0)