>> x=3+randn(500,1);
>> mean(x)
ans =
2.9648
>> std(x)
ans =
1.0134
>> y=normpdf(x,3,1);
>> plot(x,y,'.')