MST

n<-1:100
nb<-50
plot((1/sqrt(2*pi*nb))*exp(-((n-nb)^2)/(1*nb)),type="l")

for(i in 1:99){lines((1/sqrt(2*pi*nb))*exp(-((n-nb)^2)/(i*nb)),type="l")}



plot((1/sqrt(2*pi*nb))*exp(-((n-nb)^2)/(50*nb)),type="l")

for(i in 1:100){lines((1/sqrt(2*pi*nb))*exp(-((n-nb)^2)/(i*nb)),type="l")}


n<-1:100
nb<-10
i<-10
plot((1/sqrt(2*pi*nb))*exp(-((n-nb)^2)/(i*nb)),type="l",col="red",xlab="",ylab="")

for(e in 11:90){lines((1/sqrt(2*pi*e))*exp(-((n-e)^2)/(i*e)),type="l")}


n<-1:100
nb<-10
i<-10
plot((1/sqrt(2*pi*nb))*exp(-((n-nb)^2)/(i*nb)),type="l",col="red",xlab="",ylab="")
a<-0
for(e in 1:10){a<-a+10;lines((1/sqrt(2*pi*a))*exp(-((n-a)^2)/(i*a)),type="l")}


n<-1:100
nb<-10
i<-10
plot((1/sqrt(2*pi*nb))*exp(-((n-nb)^2)/(i*nb)),type="l",col="red",xlab="",ylab="")
a<-0
for(e in 1:10){a<-a+10;for(i in 1:10){lines((1/sqrt(2*pi*a))*exp(-((n-a)^2)/(i*a)),type="l")}}