clc
clc;
clear all;
close all;

c = [1.1 1.2]
H = [0.01 -0.9999*0.1*0.25;-0.9999*0.1*0.25 0.0625]
H_inv = inv(H)

e= ones(size(c))
alpha = e*H_inv*transpose(c)
sigma = c*H_inv*transpose(c)
delta = e*H_inv*transpose(e)

R_bar = alpha/delta
variance = 1/delta
std= sqrt(variance)

w= variance*H_inv*transpose(e)