clc;
clear;
close all;

%% Load Data

load HWMar26data;
X = A;
clear A B C

%% Create and Train Neural Gas Network

params.N = 20;

params.MaxIt = 30;

params.tmax = 9000;

params.epsilon_initial = 0.2; 
params.epsilon_final = 0.001; 

params.lambda_initial = 0.75;  
params.lambda_final = 0.1;  


params.T_initial = 5;
params.T_final = 30;

net = NeuralGasNetwork(X, params);
PlotResults(X,net.w,net.C);