close all;
clc;
clear;

filename2 = 'results.txt';
path2 = '/Users/sergeymakovkin/Documents/Kuramoto_network/Ver_2/Solver5_8_8_N_50_p_rew_00_p_inh_00';
Data = importdata( [ path2 '/' filename2 ] );
time = Data(:, 1);
theta = Data(:, 2:end);

omega_end(:) = ( bsxfun(@minus, theta(end, :), theta(round(end/2), :)) ) / ( time(end) - time(round(end/2)) );

Delta_wired_glial = sqrt(var(omega_end(1:end/2), 1));
Delta_random_neuron = sqrt(var(omega_end(end/2+1:end), 1));

rx_wired_glial = (1/(size(theta(:, 1:end/2), 2)))*sum(cos(theta(:, 1:end/2)')); % Sum of mean x-part of theta
ry_wired_glial = (1/(size(theta(:, 1:end/2), 2)))*sum(sin(theta(:, 1:end/2)')); % Sum of mean y-part of theta
r_wired_glial = sqrt(rx_wired_glial.*rx_wired_glial + ry_wired_glial.*ry_wired_glial);
r_wired_glial_mean = mean(r_wired_glial(:, round(size(r_wired_glial, 2)*0.5):end));

rx_random_neuron = (1/(size(theta(:, end/2+1:end), 2)))*sum(cos(theta(:, end/2+1:end)')); % Sum of mean x-part of theta
ry_random_neuron = (1/(size(theta(:, end/2+1:end), 2)))*sum(sin(theta(:, end/2+1:end)')); % Sum of mean y-part of theta
r_random_neuron = sqrt(rx_random_neuron.*rx_random_neuron + ry_random_neuron.*ry_random_neuron);
r_random_neuron_mean = mean(r_random_neuron(:, round(size(r_random_neuron, 2)*0.5):end));

fid = fopen('s_G_s_N_Delta_rho.txt', 'at');
fprintf(fid, '%f\t %f\t %f\t %f \n', [Delta_wired_glial Delta_random_neuron r_wired_glial_mean r_random_neuron_mean]');
fclose(fid);

clear;