program normal_prob
real,parameter::pi=3.141592654,e=2.718281828,interval=0.0625
real::x,y
integer::i
real,dimension(65)::fi,input_x
x=-2-interval
do i=1,65
x=x+interval
input_x(i)=x
y=abs(x)
fi(i)=(1.0/(sqrt(2.0*pi)))*(e**(-((x**2.0)/2.0)))
end do
do i=0,12
print"(a4,f9.5,f9.5,f9.5,f9.5,f9.5)","x  =",input_x(i*5+1:i*5+5)
print"(a4,f9.5,f9.5,f9.5,f9.5,f9.5)","phi=",fi(i*5+1:i*5+5)
print*
enddo
endprogram normal_prob