load('BreastCancer');

%% 1a

residual = zeros(10, 1);
what = zeros(10,1);
yhat = zeros(295,1);
for i = 1:10
    X_train = X(1:100,:);
    y_train = y(1:100,:);
    w = ista_solve(X_train, y_train, 0.1 * i.^2 );
    residual(i) = norm(X_train*w - y_train);
    what(i) = norm(w,1);

    %%y = numel(find(y_train~=X_train*w))/100;
end
plot(what, residual);