....
function [W, b] = LS()
        a1 = ann_radbas_run(P, W(1), b(1));

        z = a1;
        z($+1, 1:$) = 1;

        t = T';
        U = z';

        I = eye(size(U, 2), size(U, 2));
        x_star = inv(U'*U + ro*I) * U' * t

        W = x_star(1:$-1)';
        b = x_star($)';
    endfunction

    [W2, b2] = LS();

    W(2) = W2;
    b(2) = b2;