www.pudn.com > Patternrecognition.rar > clusteringtt.m, change:2009-05-16,size:17638b


clear all; 
close all; 
 
x1(1,1)=13.6687;  x1(1,2)=7.4379;  x1(1,3)=6.0827; 
x1(2,1)=12.2050;  x1(2,2)=5.8589;  x1(2,3)=7.8590; 
x1(3,1)=11.5415;  x1(3,2)=5.2403;  x1(3,3)=-3.6579; 
x1(4,1)=10.9466;  x1(4,2)=6.6703;  x1(4,3)=6.3669; 
x1(5,1)=11.9364;  x1(5,2)=5.2574;  x1(5,3)=8.4027; 
x1(6,1)=11.2357;  x1(6,2)=4.7387;  x1(6,3)=2.1749; 
x1(7,1)=13.5203;  x1(7,2)=6.4709;  x1(7,3)=3.1906; 
x1(8,1)=11.5270;  x1(8,2)=5.3449;  x1(8,3)=7.5019; 
x1(9,1)=12.6621;  x1(9,2)=4.4915;  x1(9,3)=1.3567; 
x1(10,1)=12.1226;  x1(10,2)=7.0985;  x1(10,3)=3.7731; 
x1(11,1)=8.8631;  x1(11,2)=6.6462;  x1(11,3)=5.7317; 
x1(12,1)=10.7487;  x1(12,2)=4.6012;  x1(12,3)=7.1788; 
x1(13,1)=11.6167;  x1(13,2)=6.5704;  x1(13,3)=4.7585; 
x1(14,1)=13.1000;  x1(14,2)=6.3040;  x1(14,3)=3.6836; 
x1(15,1)=11.5450;  x1(15,2)=4.5166;  x1(15,3)=5.7479; 
x1(16,1)=12.2808;  x1(16,2)=3.5478;  x1(16,3)=8.5439; 
x1(17,1)=12.5677;  x1(17,2)=7.6336;  x1(17,3)=7.6535; 
x1(18,1)=14.2888;  x1(18,2)=5.8923;  x1(18,3)=7.6080; 
x1(19,1)=12.4203;  x1(19,2)=2.9261;  x1(19,3)=0.5860; 
x1(20,1)=11.2617;  x1(20,2)=10.5130;  x1(20,3)=5.4427; 
x1(21,1)=12.4753;  x1(21,2)=5.1749;  x1(21,3)=3.0941; 
x1(22,1)=12.3204;  x1(22,2)=6.3387;  x1(22,3)=4.3305; 
x1(23,1)=11.7361;  x1(23,2)=8.1901;  x1(23,3)=1.0057; 
x1(24,1)=10.1189;  x1(24,2)=5.4394;  x1(24,3)=4.9840; 
x1(25,1)=12.7332;  x1(25,2)=3.4123;  x1(25,3)=0.0594; 
x1(26,1)=13.2772;  x1(26,2)=2.2310;  x1(26,3)=7.5745; 
x1(27,1)=13.3337;  x1(27,2)=8.3792;  x1(27,3)=10.0253; 
x1(28,1)=13.9313;  x1(28,2)=4.9984;  x1(28,3)=5.0184; 
x1(29,1)=11.5336;  x1(29,2)=4.1162;  x1(29,3)=6.2753; 
x1(30,1)=12.0861;  x1(30,2)=4.7394;  x1(30,3)=2.1361; 
x1(31,1)=12.4535;  x1(31,2)=7.0876;  x1(31,3)=4.6464; 
x1(32,1)=11.2932;  x1(32,2)=6.7306;  x1(32,3)=10.4726; 
x1(33,1)=10.7001;  x1(33,2)=2.4137;  x1(33,3)=2.1650; 
x1(34,1)=11.3208;  x1(34,2)=4.8488;  x1(34,3)=6.5925; 
x1(35,1)=11.4305;  x1(35,2)=6.2956;  x1(35,3)=6.3976; 
x1(36,1)=12.2733;  x1(36,2)=5.8662;  x1(36,3)=8.4433; 
x1(37,1)=12.1302;  x1(37,2)=8.6439;  x1(37,3)=11.8747; 
x1(38,1)=13.2625;  x1(38,2)=5.8481;  x1(38,3)=8.2453; 
x1(39,1)=13.4992;  x1(39,2)=6.9207;  x1(39,3)=4.8383; 
x1(40,1)=11.5195;  x1(40,2)=5.5312;  x1(40,3)=-0.0670; 
x1(41,1)=13.3727;  x1(41,2)=7.4000;  x1(41,3)=7.8219; 
x1(42,1)=12.2184;  x1(42,2)=6.0162;  x1(42,3)=7.4501; 
x1(43,1)=10.4744;  x1(43,2)=8.5130;  x1(43,3)=6.7497; 
x1(44,1)=10.4676;  x1(44,2)=7.2821;  x1(44,3)=6.7668; 
x1(45,1)=10.4011;  x1(45,2)=6.3032;  x1(45,3)=7.0190; 
x1(46,1)=12.0606;  x1(46,2)=6.0003;  x1(46,3)=5.2678; 
x1(47,1)=9.3498;  x1(47,2)=7.8227;  x1(47,3)=7.8611; 
x1(48,1)=12.7020;  x1(48,2)=5.9381;  x1(48,3)=5.3034; 
x1(49,1)=10.3735;  x1(49,2)=5.3045;  x1(49,3)=1.3789; 
x1(50,1)=12.3489;  x1(50,2)=6.2771;  x1(50,3)=1.9048; 
x1(51,1)=13.8009;  x1(51,2)=7.0795;  x1(51,3)=4.7638; 
x1(52,1)=13.9791;  x1(52,2)=12.0935;  x1(52,3)=11.8116; 
x1(53,1)=13.4456;  x1(53,2)=4.5023;  x1(53,3)=3.8819; 
x1(54,1)=12.5760;  x1(54,2)=6.2388;  x1(54,3)=7.0131; 
x1(55,1)=11.7849;  x1(55,2)=7.1185;  x1(55,3)=6.3819; 
x1(56,1)=12.3997;  x1(56,2)=8.1569;  x1(56,3)=4.7003; 
x1(57,1)=13.9923;  x1(57,2)=5.9755;  x1(57,3)=3.8076; 
x1(58,1)=11.5029;  x1(58,2)=4.2965;  x1(58,3)=7.4457; 
x1(59,1)=12.6563;  x1(59,2)=11.0555;  x1(59,3)=13.1794; 
x1(60,1)=11.4296;  x1(60,2)=6.8431;  x1(60,3)=11.2363; 
x1(61,1)=11.0676;  x1(61,2)=7.4283;  x1(61,3)=6.2538; 
x1(62,1)=12.3292;  x1(62,2)=3.7805;  x1(62,3)=1.6234; 
x1(63,1)=11.6633;  x1(63,2)=5.1718;  x1(63,3)=7.5366; 
x1(64,1)=12.0141;  x1(64,2)=7.5826;  x1(64,3)=1.3548; 
x1(65,1)=11.1760;  x1(65,2)=8.0340;  x1(65,3)=8.1898; 
x1(66,1)=12.0125;  x1(66,2)=9.9384;  x1(66,3)=3.4477; 
x1(67,1)=10.3447;  x1(67,2)=8.0100;  x1(67,3)=8.8906; 
x1(68,1)=10.9849;  x1(68,2)=5.4324;  x1(68,3)=3.7603; 
x1(69,1)=13.5125;  x1(69,2)=10.2985;  x1(69,3)=8.0724; 
x1(70,1)=12.0107;  x1(70,2)=5.8313;  x1(70,3)=7.0179; 
x1(71,1)=11.6265;  x1(71,2)=5.7369;  x1(71,3)=5.6393; 
x1(72,1)=12.0946;  x1(72,2)=6.8559;  x1(72,3)=1.3868; 
x1(73,1)=12.6534;  x1(73,2)=7.0153;  x1(73,3)=13.0916; 
x1(74,1)=12.2848;  x1(74,2)=7.6449;  x1(74,3)=10.3351; 
x1(75,1)=11.1052;  x1(75,2)=6.6351;  x1(75,3)=4.7131; 
x1(76,1)=11.5894;  x1(76,2)=8.3472;  x1(76,3)=6.3328; 
x1(77,1)=12.1890;  x1(77,2)=7.0273;  x1(77,3)=3.6144; 
x1(78,1)=12.5933;  x1(78,2)=4.5677;  x1(78,3)=10.5093; 
x1(79,1)=11.8215;  x1(79,2)=3.5661;  x1(79,3)=6.5529; 
x1(80,1)=13.4307;  x1(80,2)=7.0848;  x1(80,3)=4.1261; 
x1(81,1)=11.6796;  x1(81,2)=8.2750;  x1(81,3)=11.3900; 
x1(82,1)=10.2997;  x1(82,2)=6.1464;  x1(82,3)=4.1917; 
x1(83,1)=11.4254;  x1(83,2)=5.6168;  x1(83,3)=9.2822; 
x1(84,1)=11.0403;  x1(84,2)=8.1857;  x1(84,3)=3.2159; 
x1(85,1)=11.4751;  x1(85,2)=8.1787;  x1(85,3)=1.2492; 
x1(86,1)=11.7660;  x1(86,2)=8.9312;  x1(86,3)=9.9028; 
x1(87,1)=12.2176;  x1(87,2)=7.1134;  x1(87,3)=4.7973; 
x1(88,1)=13.3970;  x1(88,2)=2.6870;  x1(88,3)=7.4252; 
x1(89,1)=11.9241;  x1(89,2)=7.6113;  x1(89,3)=6.4029; 
x1(90,1)=12.2923;  x1(90,2)=6.3877;  x1(90,3)=-3.1790; 
x1(91,1)=12.5364;  x1(91,2)=6.8525;  x1(91,3)=4.9206; 
x1(92,1)=11.1807;  x1(92,2)=7.0602;  x1(92,3)=6.7870; 
x1(93,1)=12.4587;  x1(93,2)=11.2299;  x1(93,3)=0.7074; 
x1(94,1)=12.0199;  x1(94,2)=8.7060;  x1(94,3)=-0.4945; 
x1(95,1)=12.4906;  x1(95,2)=6.5278;  x1(95,3)=1.3834; 
x1(96,1)=12.3254;  x1(96,2)=2.1426;  x1(96,3)=5.5241; 
x1(97,1)=12.4678;  x1(97,2)=6.6505;  x1(97,3)=4.7319; 
x1(98,1)=11.8556;  x1(98,2)=7.2343;  x1(98,3)=8.1199; 
x1(99,1)=13.0847;  x1(99,2)=4.3309;  x1(99,3)=9.5755; 
x1(100,1)=12.1327;  x1(100,2)=7.3295;  x1(100,3)=10.2708; 
 
 
x2(1,1)=3.1843;  x2(1,2)=25.4958;  x2(1,3)=6.2499; 
x2(2,1)=3.5201;  x2(2,2)=26.3668;  x2(2,3)=4.9336; 
x2(3,1)=6.0144;  x2(3,2)=20.8143;  x2(3,3)=6.5988; 
x2(4,1)=7.6114;  x2(4,2)=23.2741;  x2(4,3)=6.3866; 
x2(5,1)=5.8053;  x2(5,2)=21.6878;  x2(5,3)=5.3835; 
x2(6,1)=9.0814;  x2(6,2)=23.0674;  x2(6,3)=6.1280; 
x2(7,1)=6.8167;  x2(7,2)=24.7780;  x2(7,3)=6.2330; 
x2(8,1)=3.2865;  x2(8,2)=22.4849;  x2(8,3)=5.1316; 
x2(9,1)=5.1583;  x2(9,2)=26.3514;  x2(9,3)=6.0671; 
x2(10,1)=1.1600;  x2(10,2)=22.2980;  x2(10,3)=5.5749; 
x2(11,1)=8.9414;  x2(11,2)=26.2355;  x2(11,3)=4.7787; 
x2(12,1)=1.2236;  x2(12,2)=22.7244;  x2(12,3)=6.4586; 
x2(13,1)=7.2660;  x2(13,2)=25.7167;  x2(13,3)=4.9934; 
x2(14,1)=4.9821;  x2(14,2)=24.8474;  x2(14,3)=5.9344; 
x2(15,1)=4.0812;  x2(15,2)=26.2629;  x2(15,3)=6.4116; 
x2(16,1)=2.8662;  x2(16,2)=23.8820;  x2(16,3)=5.8242; 
x2(17,1)=6.2467;  x2(17,2)=24.3826;  x2(17,3)=5.8080; 
x2(18,1)=6.8524;  x2(18,2)=20.7153;  x2(18,3)=5.0504; 
x2(19,1)=5.2066;  x2(19,2)=22.1057;  x2(19,3)=6.2690; 
x2(20,1)=7.3841;  x2(20,2)=21.2046;  x2(20,3)=6.2227; 
x2(21,1)=5.8463;  x2(21,2)=22.2282;  x2(21,3)=6.4193; 
x2(22,1)=1.9046;  x2(22,2)=25.1827;  x2(22,3)=7.0073; 
x2(23,1)=8.4363;  x2(23,2)=26.4361;  x2(23,3)=5.5118; 
x2(24,1)=3.1800;  x2(24,2)=22.5305;  x2(24,3)=6.6473; 
x2(25,1)=8.3544;  x2(25,2)=25.0643;  x2(25,3)=6.5561; 
x2(26,1)=7.0058;  x2(26,2)=24.8974;  x2(26,3)=5.1903; 
x2(27,1)=5.0480;  x2(27,2)=25.5560;  x2(27,3)=6.2479; 
x2(28,1)=2.1324;  x2(28,2)=25.6444;  x2(28,3)=4.8277; 
x2(29,1)=5.7128;  x2(29,2)=22.1297;  x2(29,3)=6.5748; 
x2(30,1)=6.4180;  x2(30,2)=21.0087;  x2(30,3)=6.0532; 
x2(31,1)=3.0957;  x2(31,2)=23.8864;  x2(31,3)=5.3477; 
x2(32,1)=5.6951;  x2(32,2)=21.8869;  x2(32,3)=4.5925; 
x2(33,1)=9.2521;  x2(33,2)=25.0839;  x2(33,3)=8.7636; 
x2(34,1)=3.8145;  x2(34,2)=23.1228;  x2(34,3)=4.4204; 
x2(35,1)=3.6982;  x2(35,2)=23.2438;  x2(35,3)=6.5376; 
x2(36,1)=7.3852;  x2(36,2)=24.2771;  x2(36,3)=6.1474; 
x2(37,1)=7.3715;  x2(37,2)=24.2567;  x2(37,3)=6.9647; 
x2(38,1)=1.8430;  x2(38,2)=23.9234;  x2(38,3)=4.6966; 
x2(39,1)=7.9294;  x2(39,2)=22.5550;  x2(39,3)=5.3630; 
x2(40,1)=7.3969;  x2(40,2)=21.3036;  x2(40,3)=5.3393; 
x2(41,1)=4.2299;  x2(41,2)=22.8500;  x2(41,3)=4.9652; 
x2(42,1)=4.3676;  x2(42,2)=23.1789;  x2(42,3)=6.6717; 
x2(43,1)=6.7703;  x2(43,2)=25.4218;  x2(43,3)=7.1126; 
x2(44,1)=9.0604;  x2(44,2)=24.6221;  x2(44,3)=5.8069; 
x2(45,1)=7.1173;  x2(45,2)=22.6033;  x2(45,3)=6.4089; 
x2(46,1)=9.9674;  x2(46,2)=23.7614;  x2(46,3)=3.8907; 
x2(47,1)=7.2099;  x2(47,2)=26.9854;  x2(47,3)=7.0986; 
x2(48,1)=4.6552;  x2(48,2)=24.3528;  x2(48,3)=5.5873; 
x2(49,1)=4.7756;  x2(49,2)=22.0086;  x2(49,3)=6.4575; 
x2(50,1)=4.6647;  x2(50,2)=25.8687;  x2(50,3)=5.5880; 
x2(51,1)=4.6709;  x2(51,2)=24.6630;  x2(51,3)=8.4769; 
x2(52,1)=8.3759;  x2(52,2)=23.9867;  x2(52,3)=5.8562; 
x2(53,1)=4.4335;  x2(53,2)=24.2433;  x2(53,3)=5.4538; 
x2(54,1)=8.1182;  x2(54,2)=21.2482;  x2(54,3)=6.4147; 
x2(55,1)=11.1130;  x2(55,2)=22.5043;  x2(55,3)=6.5645; 
x2(56,1)=6.2114;  x2(56,2)=24.4345;  x2(56,3)=5.0425; 
x2(57,1)=5.4432;  x2(57,2)=22.5287;  x2(57,3)=5.2405; 
x2(58,1)=10.1960;  x2(58,2)=26.0088;  x2(58,3)=6.5210; 
x2(59,1)=9.3881;  x2(59,2)=25.8389;  x2(59,3)=6.0721; 
x2(60,1)=5.6462;  x2(60,2)=25.9087;  x2(60,3)=5.6580; 
x2(61,1)=4.9725;  x2(61,2)=22.9453;  x2(61,3)=5.9013; 
x2(62,1)=4.9183;  x2(62,2)=26.8221;  x2(62,3)=7.8209; 
x2(63,1)=5.4790;  x2(63,2)=23.3171;  x2(63,3)=5.3990; 
x2(64,1)=5.8526;  x2(64,2)=26.2082;  x2(64,3)=6.6361; 
x2(65,1)=7.2550;  x2(65,2)=25.0171;  x2(65,3)=6.9439; 
x2(66,1)=9.6819;  x2(66,2)=21.1287;  x2(66,3)=6.5656; 
x2(67,1)=11.7152;  x2(67,2)=26.0091;  x2(67,3)=5.5865; 
x2(68,1)=3.6672;  x2(68,2)=25.5212;  x2(68,3)=5.0588; 
x2(69,1)=6.9020;  x2(69,2)=24.0887;  x2(69,3)=5.0071; 
x2(70,1)=7.6878;  x2(70,2)=26.8686;  x2(70,3)=6.4273; 
x2(71,1)=5.2503;  x2(71,2)=23.3019;  x2(71,3)=6.0763; 
x2(72,1)=6.0410;  x2(72,2)=24.7043;  x2(72,3)=6.6294; 
x2(73,1)=7.5536;  x2(73,2)=23.5847;  x2(73,3)=4.2668; 
x2(74,1)=7.2037;  x2(74,2)=26.3529;  x2(74,3)=7.1832; 
x2(75,1)=10.2444;  x2(75,2)=26.7344;  x2(75,3)=7.0040; 
x2(76,1)=7.1714;  x2(76,2)=20.5219;  x2(76,3)=6.2873; 
x2(77,1)=8.0821;  x2(77,2)=24.7492;  x2(77,3)=8.1449; 
x2(78,1)=6.9859;  x2(78,2)=23.3574;  x2(78,3)=6.7680; 
x2(79,1)=7.7146;  x2(79,2)=23.2373;  x2(79,3)=6.0963; 
x2(80,1)=11.3263;  x2(80,2)=24.7325;  x2(80,3)=4.8551; 
x2(81,1)=9.7939;  x2(81,2)=25.3192;  x2(81,3)=5.6647; 
x2(82,1)=5.6017;  x2(82,2)=25.2142;  x2(82,3)=5.4410; 
x2(83,1)=-1.1629;  x2(83,2)=22.0094;  x2(83,3)=6.2602; 
x2(84,1)=5.7512;  x2(84,2)=26.4160;  x2(84,3)=5.7084; 
x2(85,1)=10.9030;  x2(85,2)=25.3536;  x2(85,3)=5.8204; 
x2(86,1)=2.1047;  x2(86,2)=22.6004;  x2(86,3)=7.4531; 
x2(87,1)=6.6088;  x2(87,2)=25.8447;  x2(87,3)=3.6004; 
x2(88,1)=5.4351;  x2(88,2)=26.2417;  x2(88,3)=5.3689; 
x2(89,1)=8.7346;  x2(89,2)=20.4857;  x2(89,3)=4.5104; 
x2(90,1)=6.0986;  x2(90,2)=18.2535;  x2(90,3)=5.1711; 
 
 
x3(1,1)=7.0472;  x3(1,2)=4.9858;  x3(1,3)=23.2239; 
x3(2,1)=7.3158;  x3(2,2)=8.8404;  x3(2,3)=25.9957; 
x3(3,1)=10.4216;  x3(3,2)=2.8697;  x3(3,3)=23.5629; 
x3(4,1)=4.2359;  x3(4,2)=3.2804;  x3(4,3)=23.6645; 
x3(5,1)=5.7188;  x3(5,2)=12.0217;  x3(5,3)=22.8479; 
x3(6,1)=8.5223;  x3(6,2)=11.3082;  x3(6,3)=24.1968; 
x3(7,1)=4.0174;  x3(7,2)=4.3157;  x3(7,3)=24.9286; 
x3(8,1)=7.8232;  x3(8,2)=4.0656;  x3(8,3)=22.7443; 
x3(9,1)=6.6261;  x3(9,2)=8.3312;  x3(9,3)=25.0079; 
x3(10,1)=6.5049;  x3(10,2)=4.1503;  x3(10,3)=24.6356; 
x3(11,1)=6.7090;  x3(11,2)=7.4244;  x3(11,3)=22.5042; 
x3(12,1)=3.4950;  x3(12,2)=8.8357;  x3(12,3)=22.9566; 
x3(13,1)=9.1101;  x3(13,2)=4.4902;  x3(13,3)=23.4999; 
x3(14,1)=5.4134;  x3(14,2)=8.3542;  x3(14,3)=24.2257; 
x3(15,1)=6.1286;  x3(15,2)=3.9606;  x3(15,3)=22.7396; 
x3(16,1)=8.2770;  x3(16,2)=7.8868;  x3(16,3)=22.7296; 
x3(17,1)=6.3133;  x3(17,2)=3.0897;  x3(17,3)=24.4597; 
x3(18,1)=6.9705;  x3(18,2)=8.4312;  x3(18,3)=23.6748; 
x3(19,1)=2.0449;  x3(19,2)=8.5232;  x3(19,3)=23.3662; 
x3(20,1)=6.0051;  x3(20,2)=7.2427;  x3(20,3)=23.9693; 
x3(21,1)=5.2971;  x3(21,2)=1.3920;  x3(21,3)=23.6839; 
x3(22,1)=5.6138;  x3(22,2)=-0.0507;  x3(22,3)=23.0104; 
x3(23,1)=6.2404;  x3(23,2)=4.4323;  x3(23,3)=22.6639; 
x3(24,1)=7.1563;  x3(24,2)=4.1706;  x3(24,3)=23.8074; 
x3(25,1)=4.1763;  x3(25,2)=5.3487;  x3(25,3)=23.8818; 
x3(26,1)=7.4088;  x3(26,2)=10.1372;  x3(26,3)=22.3257; 
x3(27,1)=8.8329;  x3(27,2)=8.2794;  x3(27,3)=25.7425; 
x3(28,1)=5.3694;  x3(28,2)=12.9262;  x3(28,3)=22.5570; 
x3(29,1)=4.8695;  x3(29,2)=6.2980;  x3(29,3)=24.1852; 
x3(30,1)=8.2390;  x3(30,2)=4.8479;  x3(30,3)=25.4433; 
x3(31,1)=3.2932;  x3(31,2)=6.2103;  x3(31,3)=24.0331; 
x3(32,1)=5.9466;  x3(32,2)=6.4455;  x3(32,3)=22.4879; 
x3(33,1)=4.2467;  x3(33,2)=6.0882;  x3(33,3)=23.8310; 
x3(34,1)=7.6598;  x3(34,2)=6.8358;  x3(34,3)=23.4306; 
x3(35,1)=5.4593;  x3(35,2)=9.9100;  x3(35,3)=25.0607; 
x3(36,1)=7.8556;  x3(36,2)=10.0791;  x3(36,3)=22.5322; 
x3(37,1)=6.1585;  x3(37,2)=6.9780;  x3(37,3)=22.7033; 
x3(38,1)=5.4251;  x3(38,2)=3.5673;  x3(38,3)=24.6822; 
x3(39,1)=8.4130;  x3(39,2)=3.7417;  x3(39,3)=22.9200; 
x3(40,1)=4.1150;  x3(40,2)=10.1365;  x3(40,3)=23.1791; 
x3(41,1)=10.5809;  x3(41,2)=7.3908;  x3(41,3)=23.8329; 
x3(42,1)=7.6448;  x3(42,2)=8.9141;  x3(42,3)=24.7826; 
x3(43,1)=9.2826;  x3(43,2)=10.9942;  x3(43,3)=23.2796; 
x3(44,1)=9.1321;  x3(44,2)=10.3728;  x3(44,3)=23.3821; 
x3(45,1)=4.8713;  x3(45,2)=9.4258;  x3(45,3)=25.4085; 
x3(46,1)=4.2622;  x3(46,2)=8.1413;  x3(46,3)=22.5201; 
x3(47,1)=4.4076;  x3(47,2)=5.7547;  x3(47,3)=25.0195; 
x3(48,1)=3.4928;  x3(48,2)=8.4803;  x3(48,3)=23.3948; 
x3(49,1)=1.9748;  x3(49,2)=2.8132;  x3(49,3)=22.8140; 
x3(50,1)=8.6446;  x3(50,2)=8.4229;  x3(50,3)=24.2352; 
x3(51,1)=9.5094;  x3(51,2)=3.3529;  x3(51,3)=24.7325; 
x3(52,1)=5.6135;  x3(52,2)=4.6926;  x3(52,3)=23.3019; 
x3(53,1)=5.3485;  x3(53,2)=6.1470;  x3(53,3)=23.4457; 
x3(54,1)=9.1105;  x3(54,2)=9.7038;  x3(54,3)=22.9437; 
x3(55,1)=6.9126;  x3(55,2)=6.4601;  x3(55,3)=24.1031; 
x3(56,1)=4.9743;  x3(56,2)=6.5621;  x3(56,3)=24.3424; 
x3(57,1)=6.4120;  x3(57,2)=6.0273;  x3(57,3)=24.3392; 
x3(58,1)=3.1418;  x3(58,2)=5.2375;  x3(58,3)=21.8816; 
x3(59,1)=7.0881;  x3(59,2)=11.4931;  x3(59,3)=22.4968; 
x3(60,1)=7.8305;  x3(60,2)=3.5696;  x3(60,3)=23.7339; 
x3(61,1)=7.0774;  x3(61,2)=7.4042;  x3(61,3)=23.9355; 
x3(62,1)=5.6762;  x3(62,2)=10.9617;  x3(62,3)=26.9641; 
x3(63,1)=4.7352;  x3(63,2)=3.0980;  x3(63,3)=23.8276; 
x3(64,1)=7.7643;  x3(64,2)=10.8227;  x3(64,3)=24.0413; 
x3(65,1)=9.5063;  x3(65,2)=1.6303;  x3(65,3)=25.1990; 
x3(66,1)=4.8809;  x3(66,2)=4.1454;  x3(66,3)=25.4158; 
x3(67,1)=6.2552;  x3(67,2)=5.2467;  x3(67,3)=23.2922; 
x3(68,1)=8.0987;  x3(68,2)=8.2072;  x3(68,3)=22.7545; 
x3(69,1)=7.7719;  x3(69,2)=7.3520;  x3(69,3)=24.2627; 
x3(70,1)=4.1214;  x3(70,2)=4.6192;  x3(70,3)=24.5777; 
x3(71,1)=4.8790;  x3(71,2)=4.2724;  x3(71,3)=22.8663; 
x3(72,1)=9.7280;  x3(72,2)=5.4827;  x3(72,3)=24.2064; 
x3(73,1)=7.2929;  x3(73,2)=4.7312;  x3(73,3)=23.5626; 
x3(74,1)=11.2079;  x3(74,2)=9.0513;  x3(74,3)=27.3125; 
x3(75,1)=6.9272;  x3(75,2)=12.6744;  x3(75,3)=23.4386; 
x3(76,1)=7.3920;  x3(76,2)=7.6244;  x3(76,3)=23.1026; 
x3(77,1)=4.0842;  x3(77,2)=12.5081;  x3(77,3)=24.9780; 
x3(78,1)=5.7820;  x3(78,2)=8.7764;  x3(78,3)=23.5453; 
x3(79,1)=8.3525;  x3(79,2)=10.0578;  x3(79,3)=23.8787; 
x3(80,1)=6.9984;  x3(80,2)=1.5655;  x3(80,3)=25.6115; 
 
for i=1:100 r1(i)=x1(i,1);end; 
for i=1:100 r2(i)=x1(i,2);end; 
for i=1:100 r3(i)=x1(i,3);end; 
 
for i=1:90 r4(i)=x2(i,1);end; 
for i=1:90 r5(i)=x2(i,2);end; 
for i=1:90 r6(i)=x2(i,3);end; 
 
for i=1:80 r7(i)=x3(i,1);end; 
for i=1:80 r8(i)=x3(i,2);end; 
for i=1:80 r9(i)=x3(i,3);end; 
 
xnew=zeros(270,3); 
for i=1:100 xnew(i,1)=x1(i,1);end; 
for i=1:100 xnew(i,2)=x1(i,2);end; 
for i=1:100 xnew(i,3)=x1(i,3);end; 
 
for i=101:190 xnew(i,1)=x2(i-100,1);end; 
for i=101:190 xnew(i,2)=x2(i-100,2);end; 
for i=101:190 xnew(i,3)=x2(i-100,3);end; 
 
for i=191:270 xnew(i,1)=x3(i-190,1);end; 
for i=191:270 xnew(i,2)=x3(i-190,2);end; 
for i=191:270 xnew(i,3)=x3(i-190,3);end; 
 
figure(1); 
plot3(r1,r2,r3,'*',r4,r5,r6,'o',r7,r8,r9,'.'); 
axis([-5,15,-10,30,-10,30]); 
grid on; 
 
c1=zeros(270,1); c2=zeros(270,1); c3=zeros(270,1); 
 
c1(1,1)=1;  c2(1,1)=101;  c3(1,1)=191; 
k1=1;       k2=1;         k3=1; 
FailCount=0;   FailIndex=zeros(200,1); 
ctemp1=[xnew(1,1),xnew(1,2),xnew(1,3)]; 
ctemp2=[xnew(101,1),xnew(101,2),xnew(101,3)]; 
ctemp3=[xnew(191,1),xnew(191,2),xnew(191,3)]; 
 
for i=1:270 
    xtemp=[xnew(i,1),xnew(i,2),xnew(i,3)]; 
    c1dis=dis(xtemp,ctemp1); 
    c2dis=dis(xtemp,ctemp2); 
    c3dis=dis(xtemp,ctemp3); 
    if((c1dis<c2dis)&(c1dis<c3dis)) 
        k1=k1+1; 
        c1(k1,1)=i; 
        ctemp1(1)=ctemp1(1)+(xtemp(1)-ctemp1(1))/k1; 
        ctemp1(2)=ctemp1(2)+(xtemp(2)-ctemp1(2))/k1; 
        ctemp1(3)=ctemp1(3)+(xtemp(3)-ctemp1(3))/k1; 
    end; 
    if((c2dis<c1dis)&(c2dis<c3dis)) 
        k2=k2+1; 
        c2(k2,1)=i; 
        ctemp2(1)=ctemp2(1)+(xtemp(1)-ctemp2(1))/k2; 
        ctemp2(2)=ctemp2(2)+(xtemp(2)-ctemp2(2))/k2; 
        ctemp2(3)=ctemp2(3)+(xtemp(3)-ctemp2(3))/k2;  
    end; 
    if((c3dis<c1dis)&(c3dis<c2dis)) 
        k3=k3+1; 
        c3(k3,1)=i; 
        ctemp3(1)=ctemp3(1)+(xtemp(1)-ctemp3(1))/k3; 
        ctemp3(2)=ctemp3(2)+(xtemp(2)-ctemp3(2))/k3; 
        ctemp3(3)=ctemp3(3)+(xtemp(3)-ctemp3(3))/k3;         
    end; 
end; 
i=1; 
while(i<=k1) 
    if(c1(i,1)>100) 
        FailCount=FailCount+1; 
        FailIndex(FailCount,1)=c1(i,1); 
    end; 
    i=i+1; 
end; 
i=1; 
while(i<=k2) 
    if((c2(i,1)>190)|(c2(i,1)<=100)) 
        FailCount=FailCount+1; 
        FailIndex(FailCount,1)=c2(i,1); 
    end; 
    i=i+1; 
end; 
i=1; 
while(i<=k3) 
    if(c3(i,1)<=190) 
        FailCount=FailCount+1; 
        FailIndex(FailCount,1)=c3(i,1); 
    end; 
    i=i+1; 
end; 
 
e=FailCount/270 
 
if(e>0) 
      for i=1:FailCount r10(i)=xnew(FailIndex(i,1),1);end; 
      for i=1:FailCount r11(i)=xnew(FailIndex(i,1),2);end; 
      for i=1:FailCount r12(i)=xnew(FailIndex(i,1),3);end; 
      figure(2); 
      plot3(r1,r2,r3,'*',r4,r5,r6,'o',r7,r8,r9,'.',r10,r11,r12,'bd'); 
      axis([-5,15,-10,30,-10,30]); 
      grid on; 
  end;