P=[0,1,0,0;0.5,0,0.5,0;0,0.5,0,0.5;0,0,1,0]
P =
0 1.0000 0 0
0.5000 0 0.5000 0
0 0.5000 0 0.5000
0 0 1.0000 0
[V,D]=eig(P')
V =
-0.3162 -0.5000 0.3162 -0.5000
0.6325 0.5000 0.6325 -0.5000
-0.6325 0.5000 0.6325 0.5000
0.3162 -0.5000 0.3162 0.5000
D =
-1.0000 0 0 0
0 -0.5000 0 0
0 0 1.0000 0
0 0 0 0.5000
D(3,3)
ans =
1.0000
format long e
D(3,3)
ans =
1.000000000000000e+00
for i = 1:4
if (abs(D(i,i) - 1) < eps)
disp('trovato'),i
end
end
D(3,3)-1
ans =
2.220446049250313e-16
for i = 1:4
if (abs(D(i,i) - 1) < 10*eps)
disp('trovato'),i
end
end
trovato
i =
3
d=diag(D)
d =
-1.000000000000001e+00
-5.000000000000003e-01
1.000000000000000e+00
5.000000000000000e-01
find(abs(d-1)<10*eps)
ans =
3
find(abs(d-1)<10*eps,1)
ans =
3
uiopen('/home/accounts/personale/clrmrc90/aa1011/sistemi_stocastici/randstocsym.m', true);
cd aa1011/sistemi_stocastici/
randstocsym(4)
Q =
Columns 1 through 2
2.950958351198172e-01 3.280810814516226e-01
3.280810814516226e-01 4.213556508396189e-01
4.599507934921470e-02 6.499343700121384e-02
3.308280040793454e-01 1.855698307075447e-01
Columns 3 through 4
4.599507934921470e-02 3.308280040793454e-01
6.499343700121384e-02 1.855698307075447e-01
3.231771565298578e-01 5.658343271197137e-01
5.658343271197137e-01 -8.223216190660376e-02
ans =
Columns 1 through 2
2.772012767618457e-01 2.970975889632979e-01
2.970975889632979e-01 3.533597109535815e-01
1.269466339262489e-01 1.384062195651334e-01
2.987545003486074e-01 2.111364805179873e-01
Columns 3 through 4
1.269466339262489e-01 2.987545003486074e-01
1.384062195651334e-01 2.111364805179873e-01
2.941395991919183e-01 4.405075473166993e-01
4.405075473166993e-01 4.960147181670604e-02
Q=randstocsym(4)
Q =
Columns 1 through 2
6.131019199473993e-02 3.775526064530262e-01
3.775526064530262e-01 3.652190284565637e-01
3.723300676770774e-01 6.475178735355465e-02
1.888071338751565e-01 1.924765777368555e-01
Columns 3 through 4
3.723300676770774e-01 1.888071338751565e-01
6.475178735355465e-02 1.924765777368555e-01
3.018158023781478e-01 2.611023425912201e-01
2.611023425912201e-01 3.576139457967679e-01
Q =
Columns 1 through 2
6.131019199473993e-02 3.775526064530262e-01
3.775526064530262e-01 3.652190284565637e-01
3.723300676770774e-01 6.475178735355465e-02
1.888071338751565e-01 1.924765777368555e-01
Columns 3 through 4
3.723300676770774e-01 1.888071338751565e-01
6.475178735355465e-02 1.924765777368555e-01
3.018158023781478e-01 2.611023425912201e-01
2.611023425912201e-01 3.576139457967679e-01
Q=randstocsym(4)
Q =
Columns 1 through 2
2.649912479299894e-01 1.443143956631153e-02
1.443143956631153e-02 3.069084464376695e-01
3.431414806271793e-01 3.426326817092664e-01
3.774358318765199e-01 3.360274322867525e-01
Columns 3 through 4
3.431414806271793e-01 3.774358318765199e-01
3.426326817092664e-01 3.360274322867525e-01
1.176360467770594e-01 1.965897908864949e-01
1.965897908864949e-01 8.994694495023281e-02
Q =
Columns 1 through 2
2.649912479299894e-01 1.443143956631153e-02
1.443143956631153e-02 3.069084464376695e-01
3.431414806271793e-01 3.426326817092664e-01
3.774358318765199e-01 3.360274322867525e-01
Columns 3 through 4
3.431414806271793e-01 3.774358318765199e-01
3.426326817092664e-01 3.360274322867525e-01
1.176360467770594e-01 1.965897908864949e-01
1.965897908864949e-01 8.994694495023281e-02
Q=randstocsym(4)
Q =
Columns 1 through 2
6.654701691099847e-01 3.000343748303010e-02
3.000343748303010e-02 5.832390152105719e-02
2.610084366436531e-01 4.944547775716034e-01
4.351795676333220e-02 4.172178834243091e-01
Columns 3 through 4
2.610084366436531e-01 4.351795676333220e-02
4.944547775716034e-01 4.172178834243091e-01
6.118612038055479e-02 1.833506654041886e-01
1.833506654041886e-01 3.559134944081701e-01
Q =
Columns 1 through 2
6.654701691099846e-01 3.000343748303010e-02
3.000343748303010e-02 5.832390152105718e-02
2.610084366436531e-01 4.944547775716033e-01
4.351795676333219e-02 4.172178834243090e-01
Columns 3 through 4
2.610084366436531e-01 4.351795676333219e-02
4.944547775716033e-01 4.172178834243090e-01
6.118612038055478e-02 1.833506654041885e-01
1.833506654041885e-01 3.559134944081701e-01
Q=randstocsym(4)
Q =
Columns 1 through 2
1.842674360895647e-01 1.602500341243840e-01
1.602500341243840e-01 1.398351621587482e-01
3.215079868974403e-01 3.664864834946933e-01
3.339745428886109e-01 3.334283202221746e-01
Columns 3 through 4
3.215079868974403e-01 3.339745428886109e-01
3.664864834946933e-01 3.334283202221746e-01
1.487471631667545e-01 1.632583664411120e-01
1.632583664411120e-01 1.693387704481024e-01
Q =
Columns 1 through 2
1.842674360895647e-01 1.602500341243840e-01
1.602500341243840e-01 1.398351621587482e-01
3.215079868974403e-01 3.664864834946933e-01
3.339745428886110e-01 3.334283202221746e-01
Columns 3 through 4
3.215079868974403e-01 3.339745428886110e-01
3.664864834946933e-01 3.334283202221746e-01
1.487471631667546e-01 1.632583664411120e-01
1.632583664411120e-01 1.693387704481024e-01
Q=randstocsym(4)
Q =
Columns 1 through 2
1.556440968936663e-01 3.832684712760450e-01
3.832684712760450e-01 4.653437143226678e-02
3.693944706548420e-01 1.948866364429943e-01
9.169296117544656e-02 3.753105208486939e-01
Columns 3 through 4
3.693944706548420e-01 9.169296117544656e-02
1.948866364429943e-01 3.753105208486939e-01
1.602246365020409e-01 2.754942564001228e-01
2.754942564001228e-01 2.575022615757367e-01
Q =
Columns 1 through 2
1.556440968936664e-01 3.832684712760451e-01
3.832684712760451e-01 4.653437143226679e-02
3.693944706548421e-01 1.948866364429943e-01
9.169296117544658e-02 3.753105208486940e-01
Columns 3 through 4
3.693944706548421e-01 9.169296117544658e-02
1.948866364429943e-01 3.753105208486940e-01
1.602246365020409e-01 2.754942564001229e-01
2.754942564001229e-01 2.575022615757368e-01
Q=randstocsym(4)
Q =
Columns 1 through 2
3.397252579828171e-01 1.153547900658029e-01
1.153547900658029e-01 3.287436135616699e-01
2.287953235484987e-01 3.539788733323658e-01
3.161246284028814e-01 2.019227230401613e-01
Columns 3 through 4
2.287953235484987e-01 3.161246284028814e-01
3.539788733323658e-01 2.019227230401613e-01
2.008822109698540e-01 2.163435921492815e-01
2.163435921492815e-01 2.656090564076760e-01
Q =
Columns 1 through 2
3.397252579828171e-01 1.153547900658029e-01
1.153547900658029e-01 3.287436135616699e-01
2.287953235484987e-01 3.539788733323658e-01
3.161246284028814e-01 2.019227230401613e-01
Columns 3 through 4
2.287953235484987e-01 3.161246284028814e-01
3.539788733323658e-01 2.019227230401613e-01
2.008822109698540e-01 2.163435921492815e-01
2.163435921492815e-01 2.656090564076760e-01
Q=randstocsym(4)
Q =
Columns 1 through 2
3.905210198683901e-01 1.181164477825926e-01
1.181164477825926e-01 5.552775634556000e-01
3.782429092259864e-01 2.091313250361725e-01
1.131196231230309e-01 1.174746637256350e-01
Columns 3 through 4
3.782429092259864e-01 1.131196231230309e-01
2.091313250361725e-01 1.174746637256350e-01
1.194790219498355e-01 2.931467437880057e-01
2.931467437880057e-01 4.762589693633283e-01
Q =
Columns 1 through 2
3.905210198683901e-01 1.181164477825926e-01
1.181164477825926e-01 5.552775634556000e-01
3.782429092259864e-01 2.091313250361725e-01
1.131196231230309e-01 1.174746637256350e-01
Columns 3 through 4
3.782429092259864e-01 1.131196231230309e-01
2.091313250361725e-01 1.174746637256350e-01
1.194790219498355e-01 2.931467437880057e-01
2.931467437880057e-01 4.762589693633283e-01
Q=randstocsym(4)
Q =
Columns 1 through 2
1.517424699073784e-01 3.585058431085063e-01
3.585058431085063e-01 3.001553399877258e-01
2.525437730415233e-01 9.354197620330169e-02
2.372079139425921e-01 2.477968407004663e-01
Columns 3 through 4
2.525437730415233e-01 2.372079139425921e-01
9.354197620330169e-02 2.477968407004663e-01
4.345662895819380e-01 2.193479611732371e-01
2.193479611732371e-01 2.956472841837046e-01
Q =
Columns 1 through 2
1.517424699073784e-01 3.585058431085063e-01
3.585058431085063e-01 3.001553399877259e-01
2.525437730415233e-01 9.354197620330171e-02
2.372079139425921e-01 2.477968407004663e-01
Columns 3 through 4
2.525437730415233e-01 2.372079139425921e-01
9.354197620330171e-02 2.477968407004663e-01
4.345662895819381e-01 2.193479611732371e-01
2.193479611732371e-01 2.956472841837046e-01
Q=randstocsym(4)
Q =
Columns 1 through 2
5.268503453142457e-02 3.747136623778635e-02
3.747136623778635e-02 5.506153600717013e-01
3.686685033972116e-01 7.658194382791572e-02
5.411750958335776e-01 3.353313298625966e-01
Columns 3 through 4
3.686685033972116e-01 5.411750958335776e-01
7.658194382791572e-02 3.353313298625966e-01
5.410309412354743e-01 1.371861153939830e-02
1.371861153939830e-02 1.097749627644274e-01
Q =
Columns 1 through 2
5.268503453142457e-02 3.747136623778635e-02
3.747136623778635e-02 5.506153600717013e-01
3.686685033972116e-01 7.658194382791572e-02
5.411750958335776e-01 3.353313298625966e-01
Columns 3 through 4
3.686685033972116e-01 5.411750958335776e-01
7.658194382791572e-02 3.353313298625966e-01
5.410309412354743e-01 1.371861153939830e-02
1.371861153939830e-02 1.097749627644274e-01
Q=randstocsym(4)
Q =
Columns 1 through 2
9.029114321337552e-02 4.421990290680108e-01
4.421990290680108e-01 8.965539184237223e-02
1.732615736233241e-01 3.258155193907905e-01
2.942482540952896e-01 1.423300596988266e-01
Columns 3 through 4
1.732615736233241e-01 2.942482540952896e-01
3.258155193907905e-01 1.423300596988266e-01
2.439103432141307e-01 2.570125637717548e-01
2.570125637717548e-01 3.064091224341290e-01
Q =
Columns 1 through 2
9.029114321337552e-02 4.421990290680108e-01
4.421990290680108e-01 8.965539184237223e-02
1.732615736233241e-01 3.258155193907905e-01
2.942482540952896e-01 1.423300596988266e-01
Columns 3 through 4
1.732615736233241e-01 2.942482540952896e-01
3.258155193907905e-01 1.423300596988266e-01
2.439103432141307e-01 2.570125637717548e-01
2.570125637717548e-01 3.064091224341290e-01
Q=randstocsym(4)
Q =
Columns 1 through 2
2.687109777756991e-01 4.999255233018701e-02
4.999255233018701e-02 9.545445424176846e-02
1.365659139056634e-01 5.173181107294291e-01
5.447305559884504e-01 3.372348826986155e-01
Columns 3 through 4
1.365659139056634e-01 5.447305559884504e-01
5.173181107294291e-01 3.372348826986155e-01
3.209297691452105e-01 2.518620621969684e-02
2.518620621969684e-02 9.284835509323730e-02
Q =
Columns 1 through 2
2.687109777756991e-01 4.999255233018701e-02
4.999255233018701e-02 9.545445424176846e-02
1.365659139056634e-01 5.173181107294291e-01
5.447305559884504e-01 3.372348826986155e-01
Columns 3 through 4
1.365659139056634e-01 5.447305559884504e-01
5.173181107294291e-01 3.372348826986155e-01
3.209297691452105e-01 2.518620621969684e-02
2.518620621969684e-02 9.284835509323730e-02
Q=randstocsym(4)
Q =
Columns 1 through 2
5.770955444621647e-02 5.204807111416343e-01
5.204807111416343e-01 2.213653317351558e-01
2.507567403341779e-03 2.352846370477091e-01
4.193021670088074e-01 2.286932007550082e-02
Columns 3 through 4
2.507567403341779e-03 4.193021670088074e-01
2.352846370477091e-01 2.286932007550082e-02
4.619359385480312e-01 3.002718570009177e-01
3.002718570009177e-01 2.575566559147741e-01
Q =
Columns 1 through 2
5.770955444621647e-02 5.204807111416343e-01
5.204807111416343e-01 2.213653317351558e-01
2.507567403341779e-03 2.352846370477091e-01
4.193021670088074e-01 2.286932007550082e-02
Columns 3 through 4
2.507567403341779e-03 4.193021670088074e-01
2.352846370477091e-01 2.286932007550082e-02
4.619359385480312e-01 3.002718570009177e-01
3.002718570009177e-01 2.575566559147741e-01
Q=randstocsym(4)
Q =
Columns 1 through 2
2.413218573731981e-01 5.093929654323317e-01
5.093929654323317e-01 6.204054504591330e-02
1.017205750818202e-01 5.800348142870333e-02
1.475646021126500e-01 3.705630080930516e-01
Columns 3 through 4
1.017205750818202e-01 1.475646021126500e-01
5.800348142870333e-02 3.705630080930516e-01
4.312384944629450e-01 4.090374490265314e-01
4.090374490265314e-01 7.283494076776686e-02
Q =
Columns 1 through 2
2.413218573731981e-01 5.093929654323317e-01
5.093929654323317e-01 6.204054504591330e-02
1.017205750818202e-01 5.800348142870333e-02
1.475646021126500e-01 3.705630080930516e-01
Columns 3 through 4
1.017205750818202e-01 1.475646021126500e-01
5.800348142870333e-02 3.705630080930516e-01
4.312384944629450e-01 4.090374490265314e-01
4.090374490265314e-01 7.283494076776686e-02
Q=randstocsym(4)
Q =
Columns 1 through 2
3.646541408082222e-01 2.659164182812208e-01
2.659164182812208e-01 4.109856892600525e-01
1.500252877768920e-01 7.770183899469055e-02
2.194041531336649e-01 2.453960534640361e-01
Columns 3 through 4
1.500252877768920e-01 2.194041531336649e-01
7.770183899469055e-02 2.453960534640361e-01
3.099856261272313e-01 4.622872471011861e-01
4.622872471011861e-01 7.291254630111288e-02
Q =
Columns 1 through 2
3.646541408082222e-01 2.659164182812208e-01
2.659164182812208e-01 4.109856892600525e-01
1.500252877768920e-01 7.770183899469055e-02
2.194041531336649e-01 2.453960534640361e-01
Columns 3 through 4
1.500252877768920e-01 2.194041531336649e-01
7.770183899469055e-02 2.453960534640361e-01
3.099856261272313e-01 4.622872471011861e-01
4.622872471011861e-01 7.291254630111288e-02
Q=randstocsym(4)
Q =
Columns 1 through 2
1.802897502896524e-01 2.145432867215529e-02
2.145432867215529e-02 3.644860829538384e-01
3.900390759237341e-01 3.632895143106232e-01
4.082168451144581e-01 2.507700740633833e-01
Columns 3 through 4
3.900390759237341e-01 4.082168451144581e-01
3.632895143106232e-01 2.507700740633833e-01
1.749132918901111e-01 7.175811787553155e-02
7.175811787553155e-02 2.692549629466271e-01
Q =
Columns 1 through 2
1.802897502896524e-01 2.145432867215529e-02
2.145432867215529e-02 3.644860829538384e-01
3.900390759237342e-01 3.632895143106232e-01
4.082168451144582e-01 2.507700740633833e-01
Columns 3 through 4
3.900390759237342e-01 4.082168451144582e-01
3.632895143106232e-01 2.507700740633833e-01
1.749132918901112e-01 7.175811787553156e-02
7.175811787553156e-02 2.692549629466272e-01
Q=randstocsym(4)
Q =
Columns 1 through 2
4.297502214636057e-01 2.146618487759404e-01
2.146618487759404e-01 6.471663032790362e-02
1.331196517950703e-01 8.854811918725378e-02
2.224682779653837e-01 6.320734017089020e-01
Columns 3 through 4
1.331196517950703e-01 2.224682779653837e-01
8.854811918725378e-02 6.320734017089020e-01
4.859722884956435e-01 2.923599405220325e-01
2.923599405220325e-01 -1.469016201963182e-01
Q =
Columns 1 through 2
3.326356888448396e-01 2.337541648342585e-01
2.337541648342585e-01 1.648205840322328e-01
1.962671700335309e-01 1.757765176720100e-01
2.373429762873710e-01 4.256487334614987e-01
Columns 3 through 4
1.962671700335309e-01 2.373429762873710e-01
1.757765176720100e-01 4.256487334614987e-01
3.584823843295173e-01 2.694739279649418e-01
2.694739279649418e-01 6.753436228618849e-02
sum(Q,2)
ans =
1.000000000000000e+00
9.999999999999999e-01
1.000000000000000e+00
9.999999999999999e-01
format short e
sum(Q,2)
ans =
1.0000e+00
1.0000e+00
1.0000e+00
1.0000e+00
help repmat
REPMAT Replicate and tile an array.
B = repmat(A,M,N) creates a large matrix B consisting of an M-by-N
tiling of copies of A. The size of B is [size(A,1)*M, size(A,2)*N].
The statement repmat(A,N) creates an N-by-N tiling.
B = REPMAT(A,[M N]) accomplishes the same result as repmat(A,M,N).
B = REPMAT(A,[M N P ...]) tiles the array A to produce a
multidimensional array B composed of copies of A. The size of B is
[size(A,1)*M, size(A,2)*N, size(A,3)*P, ...].
REPMAT(A,M,N) when A is a scalar is commonly used to produce an M-by-N
matrix filled with A's value and having A's CLASS. For certain values,
you may achieve the same results using other functions. Namely,
REPMAT(NAN,M,N) is the same as NAN(M,N)
REPMAT(SINGLE(INF),M,N) is the same as INF(M,N,'single')
REPMAT(INT8(0),M,N) is the same as ZEROS(M,N,'int8')
REPMAT(UINT32(1),M,N) is the same as ONES(M,N,'uint32')
REPMAT(EPS,M,N) is the same as EPS(ONES(M,N))
Example:
repmat(magic(2), 2, 3)
repmat(uint8(5), 2, 3)
Class support for input A:
float: double, single
See also bsxfun, meshgrid, ones, zeros, nan, inf.
Overloaded methods:
categorical/repmat
Reference page in Help browser
doc repmat
v=rand(1,3);
v=v/sum(v)
v =
1.6661e-01 2.5062e-01 5.8277e-01
repmat(v,3,1)
ans =
1.6661e-01 2.5062e-01 5.8277e-01
1.6661e-01 2.5062e-01 5.8277e-01
1.6661e-01 2.5062e-01 5.8277e-01
repmat(v,3,2)
ans =
Columns 1 through 4
1.6661e-01 2.5062e-01 5.8277e-01 1.6661e-01
1.6661e-01 2.5062e-01 5.8277e-01 1.6661e-01
1.6661e-01 2.5062e-01 5.8277e-01 1.6661e-01
Columns 5 through 6
2.5062e-01 5.8277e-01
2.5062e-01 5.8277e-01
2.5062e-01 5.8277e-01
Vrepmat(v,3,1)
??? Undefined function or method 'Vrepmat' for input
arguments of type 'double'.
V=repmat(v,3,1)
V =
1.6661e-01 2.5062e-01 5.8277e-01
1.6661e-01 2.5062e-01 5.8277e-01
1.6661e-01 2.5062e-01 5.8277e-01
W=repmat(v',1,3)
W =
1.6661e-01 1.6661e-01 1.6661e-01
2.5062e-01 2.5062e-01 2.5062e-01
5.8277e-01 5.8277e-01 5.8277e-01
V'
ans =
1.6661e-01 1.6661e-01 1.6661e-01
2.5062e-01 2.5062e-01 2.5062e-01
5.8277e-01 5.8277e-01 5.8277e-01
V
V =
1.6661e-01 2.5062e-01 5.8277e-01
1.6661e-01 2.5062e-01 5.8277e-01
1.6661e-01 2.5062e-01 5.8277e-01
W
W =
1.6661e-01 1.6661e-01 1.6661e-01
2.5062e-01 2.5062e-01 2.5062e-01
5.8277e-01 5.8277e-01 5.8277e-01
index = find(V>=W)
index =
1
4
5
7
8
9
P=zeros(3)
P =
0 0 0
0 0 0
0 0 0
Q=randstocsym(3)
Q =
6.7732e-02 1.8918e-01 7.4308e-01
1.8918e-01 3.8115e-01 4.2966e-01
7.4308e-01 4.2966e-01 -1.7274e-01
Q =
2.0291e-01 2.6255e-01 5.3454e-01
2.6255e-01 3.5682e-01 3.8063e-01
5.3454e-01 3.8063e-01 8.4826e-02
P(index)=Q(index)
P =
2.0291e-01 2.6255e-01 5.3454e-01
0 3.5682e-01 3.8063e-01
0 0 8.4826e-02
v
v =
1.6661e-01 2.5062e-01 5.8277e-01
P=metropolis(v)
P =
3.4840e-01 4.2264e-01 2.2895e-01
2.8097e-01 6.0218e-01 1.1685e-01
6.5457e-02 5.0251e-02 8.8429e-01
invariantPotenze(P,1e-6,100)
ans =
1.6661e-01 2.5062e-01 5.8276e-01
w=rand(1,3);
w=w/sum(w)
w =
4.3570e-01 4.5303e-02 5.1900e-01
w=w*P
w =
1.9850e-01 2.3751e-01 5.6399e-01
w=w*P
w =
1.7281e-01 2.5526e-01 5.7193e-01
w=w*P
w =
1.6937e-01 2.5549e-01 5.7515e-01
w=w*P
w =
1.6844e-01 2.5433e-01 5.7723e-01
w=w*P
w =
1.6793e-01 2.5335e-01 5.7872e-01
w=w*P
w =
1.6757e-01 2.5262e-01 5.7981e-01
w=w*P
w =
1.6732e-01 2.5208e-01 5.8060e-01
w=w*P
w =
1.6713e-01 2.5169e-01 5.8119e-01
w=w*P
w =
1.6699e-01 2.5140e-01 5.8161e-01
w=w*P
w =
1.6689e-01 2.5119e-01 5.8192e-01
w=w*P
w =
1.6681e-01 2.5104e-01 5.8215e-01
w=w*P
w =
1.6676e-01 2.5093e-01 5.8231e-01
w=w*P
w =
1.6672e-01 2.5084e-01 5.8244e-01
w=w*P
w =
1.6669e-01 2.5078e-01 5.8252e-01
v
v =
1.6661e-01 2.5062e-01 5.8277e-01
P^100
ans =
1.6661e-01 2.5062e-01 5.8277e-01
1.6661e-01 2.5062e-01 5.8277e-01
1.6661e-01 2.5062e-01 5.8277e-01
exp(-(-20/(1/10)))
ans =
7.2260e+86
exp(-(-20/(1/100)))
ans =
Inf
v=[1e400,2e400]
v =
Inf Inf
sum(v)
ans =
Inf
v/sum(v)
ans =
NaN NaN
exp(-(-20/(1/100)))
ans =
Inf
exp(-(-21/(1/100)))
ans =
Inf
exp(-((-21+20)/(1/100)))
ans =
2.6881e+43
exp(-((-20+21)/(1/100)))
ans =
3.7201e-44
sqrt((3e200)^2+(4e200)^2)
ans =
Inf
5e200
ans =
5.0000e+200
norm([3e200,4e200])
ans =
5.0000e+200
epsilonmetropolis
epsilon =
10
ans =
Columns 1 through 4
3.6340e-01 7.6340e-02 3.7931e-02 7.5504e-02
9.7703e-02 2.4094e-01 5.2816e-02 9.7906e-02
9.6342e-02 1.0482e-01 1.9103e-01 9.8502e-02
9.6634e-02 9.7906e-02 4.9635e-02 2.4414e-01
3.1022e-02 2.4372e-02 1.3976e-02 2.5898e-02
9.8283e-02 1.0798e-01 4.8001e-02 9.5081e-02
9.8344e-02 1.0641e-01 1.2339e-01 1.1552e-01
1.0654e-01 9.8335e-02 4.7176e-02 9.9056e-02
4.9149e-03 3.6298e-03 1.7363e-03 3.2940e-03
1.0191e-01 9.9250e-02 4.6620e-02 1.1490e-01
Columns 5 through 8
9.2882e-02 7.6793e-02 8.3400e-03 8.3247e-02
9.3391e-02 1.0798e-01 1.1550e-02 9.8335e-02
1.0628e-01 9.5259e-02 2.6578e-02 9.3622e-02
9.9240e-02 9.5081e-02 1.2538e-02 9.9056e-02
7.2449e-01 2.6207e-02 2.5118e-03 2.4516e-02
1.0042e-01 2.4065e-01 1.1908e-02 1.1421e-01
8.8678e-02 1.0972e-01 8.8368e-02 8.8642e-02
9.3944e-02 1.1421e-01 9.6210e-03 2.1400e-01
1.3800e-02 3.4040e-03 3.5390e-04 3.4973e-03
1.0601e-01 8.9569e-02 9.8750e-03 1.2064e-01
Columns 9 through 10
1.0593e-01 7.9628e-02
1.0013e-01 9.9250e-02
9.5050e-02 9.2520e-02
9.0865e-02 1.1490e-01
9.9343e-02 2.7666e-02
9.3899e-02 8.9569e-02
8.9947e-02 9.0983e-02
9.6473e-02 1.2064e-01
9.6023e-01 5.1428e-03
1.4187e-01 1.6936e-01
guarda P(5,5) e P(9,9)
epsilon =
1.0000e-01
ans =
Columns 1 through 4
8.0118e-01 1.8798e-12 3.1765e-42 1.8593e-12
9.7703e-02 3.0531e-01 1.7962e-31 9.7906e-02
9.6342e-02 1.0482e-01 2.1761e-01 9.8502e-02
9.6634e-02 9.7906e-02 1.6880e-31 3.0632e-01
2.1989e-49 4.2540e-60 8.2960e-90 4.5204e-60
9.8283e-02 1.0798e-01 1.6324e-31 9.5081e-02
9.8344e-02 1.0641e-01 1.2339e-01 1.1552e-01
1.0654e-01 9.8335e-02 1.6044e-31 9.9056e-02
4.7100e-135 8.5656e-146 1.3934e-175 7.7732e-146
1.0191e-01 9.9250e-02 1.5855e-31 1.1490e-01
Columns 5 through 8
9.2882e-02 1.8910e-12 6.8372e-109 2.0499e-12
9.3391e-02 1.0798e-01 3.8451e-98 9.8335e-02
1.0628e-01 9.5259e-02 2.6018e-68 9.3622e-02
9.9240e-02 9.5081e-02 4.1742e-98 9.9056e-02
9.0066e-01 4.5743e-60 1.4596e-156 4.2792e-60
1.0042e-01 3.0055e-01 3.9645e-98 1.1421e-01
8.8678e-02 1.0972e-01 8.8368e-02 8.8642e-02
9.3944e-02 1.1421e-01 3.2030e-98 2.7080e-01
1.8657e-87 8.0327e-146 2.7804e-242 8.2529e-146
1.0601e-01 8.9569e-02 3.2876e-98 1.2064e-01
Columns 9 through 10
1.0593e-01 1.9608e-12
1.0013e-01 9.9250e-02
9.5050e-02 9.2520e-02
9.0865e-02 1.1490e-01
9.9343e-02 4.8290e-60
9.3899e-02 8.9569e-02
8.9947e-02 9.0983e-02
9.6473e-02 1.2064e-01
1.0000e+00 1.2136e-145
1.4187e-01 2.2585e-01
guarda P(5,5) e P(9,9)
epsilon =
1.0000e-02
ans =
Columns 1 through 4
8.0118e-01 6.7927e-109 0 6.7183e-109
9.7703e-02 3.0531e-01 2.2892e-299 9.7906e-02
9.6342e-02 1.0482e-01 2.1761e-01 9.8502e-02
9.6634e-02 9.7906e-02 2.1513e-299 3.0632e-01
0 0 0 0
9.8283e-02 1.0798e-01 2.0805e-299 9.5081e-02
9.8344e-02 1.0641e-01 1.2339e-01 1.1552e-01
1.0654e-01 9.8335e-02 2.0447e-299 9.9056e-02
0 0 0 0
1.0191e-01 9.9250e-02 2.0206e-299 1.1490e-01
Columns 5 through 8
9.2882e-02 6.8330e-109 0 7.4072e-109
9.3391e-02 1.0798e-01 0 9.8335e-02
1.0628e-01 9.5259e-02 0 9.3622e-02
9.9240e-02 9.5081e-02 0 9.9056e-02
9.0066e-01 0 0 0
1.0042e-01 3.0055e-01 0 1.1421e-01
8.8678e-02 1.0972e-01 8.8368e-02 8.8642e-02
9.3944e-02 1.1421e-01 0 2.7080e-01
0 0 0 0
1.0601e-01 8.9569e-02 0 1.2064e-01
Columns 9 through 10
1.0593e-01 7.0853e-109
1.0013e-01 9.9250e-02
9.5050e-02 9.2520e-02
9.0865e-02 1.1490e-01
9.9343e-02 0
9.3899e-02 8.9569e-02
8.9947e-02 9.0983e-02
9.6473e-02 1.2064e-01
1.0000e+00 0
1.4187e-01 2.2585e-01
guarda P(5,5) e P(9,9)
epsilonmetropolis
epsilon =
10
ans =
Columns 1 through 4
3.1970e-01 8.7963e-02 1.8180e-02 6.9766e-02
1.1258e-01 2.8619e-01 3.4268e-02 6.3178e-02
4.6174e-02 6.8006e-02 2.1633e-01 8.7645e-02
8.9290e-02 6.3178e-02 4.4164e-02 3.0801e-01
6.2596e-02 4.4332e-02 2.3360e-02 3.1375e-02
1.0641e-01 6.4889e-03 1.7983e-02 1.5195e-01
1.0141e-01 1.5220e-01 1.8849e-01 1.4980e-01
4.5067e-02 1.5703e-01 1.4060e-02 5.6185e-02
4.4140e-03 4.9627e-03 3.1210e-03 4.2590e-03
1.2144e-01 1.6973e-02 6.5079e-02 3.3259e-02
Columns 5 through 8
1.8741e-01 8.3140e-02 8.6002e-03 3.5213e-02
1.6988e-01 6.4889e-03 1.6519e-02 1.5703e-01
1.7765e-01 3.5688e-02 4.0600e-02 2.7902e-02
1.2022e-01 1.5195e-01 1.6259e-02 5.6185e-02
6.9093e-01 2.0483e-02 1.4940e-03 1.3211e-02
7.8488e-02 3.0017e-01 1.9041e-02 1.6442e-01
5.2744e-02 1.7543e-01 2.6243e-02 9.6965e-02
5.0624e-02 1.6442e-01 1.0524e-02 1.3657e-01
1.3257e-02 4.2675e-03 1.2340e-05 7.6671e-03
6.4331e-02 3.7338e-02 5.8159e-03 1.5402e-01
Columns 9 through 10
9.5138e-02 9.4886e-02
1.3690e-01 1.6973e-02
1.7086e-01 1.2915e-01
1.1749e-01 3.3259e-02
9.5433e-02 1.6788e-02
1.1772e-01 3.7338e-02
3.1363e-03 5.3585e-02
2.1150e-01 1.5402e-01
9.5713e-01 9.1108e-04
2.5133e-02 4.7661e-01
guarda P(5,5) e P(9,9)
epsilon =
1.0000e-01
ans =
Columns 1 through 4
7.1745e-01 2.1661e-12 1.5224e-42 1.7180e-12
1.1258e-01 3.3698e-01 1.1654e-31 6.3178e-02
4.6174e-02 6.8006e-02 2.5693e-01 8.7645e-02
8.9290e-02 6.3178e-02 1.5019e-31 3.6843e-01
4.4369e-49 7.7379e-60 1.3867e-89 5.4763e-60
1.0641e-01 6.4889e-03 6.1157e-32 1.5195e-01
1.0141e-01 1.5220e-01 1.8849e-01 1.4980e-01
4.5067e-02 1.5703e-01 4.7815e-32 5.6185e-02
4.2300e-135 1.1711e-145 2.5047e-175 1.0050e-145
1.2144e-01 1.6973e-02 2.2132e-31 3.3259e-02
Columns 5 through 8
1.8741e-01 2.0473e-12 7.0505e-109 8.6711e-13
1.6988e-01 6.4889e-03 5.4995e-98 1.5703e-01
1.7765e-01 3.5688e-02 3.9745e-68 2.7902e-02
1.2022e-01 1.5195e-01 5.4128e-98 5.6185e-02
9.0457e-01 3.5752e-60 8.6813e-157 2.3059e-60
7.8488e-02 3.3719e-01 6.3390e-98 1.6442e-01
5.2744e-02 1.7543e-01 2.6243e-02 9.6965e-02
5.0624e-02 1.6442e-01 3.5038e-98 1.6116e-01
1.7923e-87 1.0070e-145 9.6949e-244 1.8093e-145
6.4331e-02 3.7338e-02 1.9362e-98 1.5402e-01
Columns 9 through 10
9.5138e-02 2.3365e-12
1.3690e-01 1.6973e-02
1.7086e-01 1.2915e-01
1.1749e-01 3.3259e-02
9.5433e-02 2.9303e-60
1.1772e-01 3.7338e-02
3.1363e-03 5.3585e-02
2.1150e-01 1.5402e-01
1.0000e+00 2.1500e-146
2.5133e-02 5.4751e-01
guarda P(5,5) e P(9,9)
epsilon =
1.0000e-02
ans =
Columns 1 through 4
7.1745e-01 7.8269e-109 0 6.2078e-109
1.1258e-01 3.3698e-01 1.4853e-299 6.3178e-02
4.6174e-02 6.8006e-02 2.5693e-01 8.7645e-02
8.9290e-02 6.3178e-02 1.9142e-299 3.6843e-01
0 0 0 0
1.0641e-01 6.4889e-03 7.7943e-300 1.5195e-01
1.0141e-01 1.5220e-01 1.8849e-01 1.4980e-01
4.5067e-02 1.5703e-01 6.0938e-300 5.6185e-02
0 0 0 0
1.2144e-01 1.6973e-02 2.8207e-299 3.3259e-02
Columns 5 through 8
1.8741e-01 7.3977e-109 0 3.1332e-109
1.6988e-01 6.4889e-03 0 1.5703e-01
1.7765e-01 3.5688e-02 0 2.7902e-02
1.2022e-01 1.5195e-01 0 5.6185e-02
9.0457e-01 0 0 0
7.8488e-02 3.3719e-01 0 1.6442e-01
5.2744e-02 1.7543e-01 2.6243e-02 9.6965e-02
5.0624e-02 1.6442e-01 0 1.6116e-01
0 0 0 0
6.4331e-02 3.7338e-02 0 1.5402e-01
Columns 9 through 10
9.5138e-02 8.4429e-109
1.3690e-01 1.6973e-02
1.7086e-01 1.2915e-01
1.1749e-01 3.3259e-02
9.5433e-02 0
1.1772e-01 3.7338e-02
3.1363e-03 5.3585e-02
2.1150e-01 1.5402e-01
1.0000e+00 0
2.5133e-02 5.4751e-01
guarda P(5,5) e P(9,9)
diary off