clear all
close all
clc


N = 100000;
a0 = 5;
a1 = 3;
a2 = -6;

sig2 = 100;

x1 = rand(N,1);
x2 = -2+4*rand(N,1);
eps = sqrt(sig2)*randn(N,1);

y = a0+a1*x1+a2*x2+eps;

yy = y;
XX = [ones(size(yy)) x1 x2];

result = ols(yy,XX);
vnames = strvcat('y','konst','x1','x2');
prt_reg(result,vnames)
plt_reg(result)