clear
clc

x=[27 27 39 40 44 44 49 58 71 84 85 98 109 131 128 133 134 161 194]';
%x=[11318 13213 14369 16392 16683 16952 18394 19275 19505 19765 16968 14385 13650 15441 17000 18750 19453 21400 21782 20906 21097 23996 25014 24048 26494 27527 28452 28631 28273 30477 33212 32056 35450 38728 40731 37911 39151 40473 39408 40755 44624 43529 44266 45649 44510 46662 50454 49417 51230 50839 46251 45262 45711]';
%[P,yy]=mcgm1(x,0);
[Ymcgm,Ygm,P1,ERROR]=mcgm(x)