a=[87 63 71 81 67 68 71 78 68 81 67 68 68 79 65 75 72 72 65 72 78 68 73 71 64 92 74 68 76 73 65 66 70 66 63 84 78 75 67 68 84 63 75 74 71 73 83 66 80 68 60 58 54 62 62 82 52 54 53 59 77 58 63 59 57 82 58 59 67 53 56 59 89 62 56 54 59 61 56 66 59 59 66 54 68 63 62 63 64 58 60 72 58 56 67 64 53 60 63 65 70 53 58 55 63 72 52 59 56 66 68 54 75 55 63 64 62 61 75 56 57 64 55 57 75 57 57 60 55 56 57 57 57 59 57 65 52 58 75 60 60 62 52 56 75 65 70 53 58 55 75 70 65 64 68 82 68 68 64 66 68 66 71 67 87 66 91 70 82 66 71 66 66 66 72 65 64 68 66 66 69 66 72 77 62 75 74 89 65 77 74 76 70 65 63 70 74 66 73 78 82 65 67 62 87 70 65 62 68 70 64 68 64 75 64 81 68 66 79 80 71 69 65 72 72 64 71 75 81 88 70 65 84 67 87 76 66 69 66 67 67 76 79 64 65 94 91 73 72 64 63 92 75 65 76 71 80 80 67 80 72 70 67 63 75 72 63 79 69 63 70 65 102 62 76 80 69 90 83 67 79 66 64 66 70 78 76 69 86 64 69 95 85 76 77 83 64 71 76 76 71 76 71 66 82 64 76 75 62 69 63 72 63 77 84 73 70 73 72 77 90 78 75 76 75 72 66 72 66 71 67 62 67 63 75 65 66 66 79 63 71 63 63 71 66 77 68 70 66 74 65 82 70 61 86 81 93 76 68 66 74 70 94 65 75 66 64 73 78 86 70 67 61 69 75 72 69 71 62 84 74 71 73 65 67 61 77 65 66 82 70 66 62 67 65 73 85 79 71 67 69 65 70 83 68 80 66 68 68 68 65 61 72 64 64 64 82 70 87 72 63 68 78 63 75 70 66 69 81 73 67 67 70 73 71 80 65 62 81 77 77 81 82 75 63 70 75 61 73 81 68 83 69 71 62 65 64 87 72 64 64 68 71 67 76 70 78 62 64 65 63 73 68 67 64 67 69 68 86 74 71 73 77 66 84 87 70 66 81 62 78 74 69 66 76 63 72 68 63 61 ];
w=a(1:500,1);
n=500;
pdf = @(x)power(2*pi*8.49*8.49,-n/2)*exp(-0.5*sum(power((w-x)/8.49,2)))*((x-50)/50);
proppdf = @(x,y) unifpdf(y-x,-1000,1000); %proposal pdf q(y|x)
proprnd = @(x)unifrnd(50,100);%the range is 2 times of the one set by the minimun and maximum measured values
nsamples =15000;
[x,accept]= mhsample(78,nsamples,'pdf',pdf,'proprnd',proprnd,'proppdf',proppdf,'thin',21,'burnin',2000);