
function [result] = simple_linear(X_l, Y_l, X_t, Y_t, X_eval, Y_eval, bw_list)
    
    X=sparse([X_l;X_t]);
    X=[ones(size(X,1),1),X];
    warning('off','MATLAB:rankDeficientMatrix');
    beta=spqr_solve (X'*X,X'*[Y_l;Y_t],struct ('solution','min2norm'));
    warning('on','MATLAB:rankDeficientMatrix');
    %beta=regress([Y_l;Y_t],[X_l;X_t]);       
    Y_calc=max([ones(size(X_eval,1),1),X_eval]*beta,0);

    result.mse=mean((Y_calc-Y_eval).^2);
    result.res=Y_calc-Y_eval;
    result.param=bw_list(1);
        
end

