Convert Matlab code to vb or c#

case ('BCC')
switch orientation;
    % Input oriented;
    case ('io')
        Z = zeros(n,n+m+s+1);

        % Objective function of the BCC model: min(0*lambda - epsilon*(s+ + s-) + theta);
        f = [zeros(1,n) -epsilon*ones(1,s+m) 1];

        lblambda = zeros(n,1);                % Lower bounds for (n) lambdas;
        lboutput = zeros(s,1);                % Lower bounds for (s) outputs;
        lbinput = zeros(m,1);                 % Lower bounds for (m) inputs ;
        lb = [lblambda; lboutput; lbinput];   % Lower bounds for lambdas, outputs (s+) and inputs (s-);
        for j=1:n
            Aeq = [Y', -eye(s,s), zeros(s,m+1);
                  -X', zeros(m,s), -eye(m,m) X(j,:)';
                  ones(1,n), zeros(1,s), zeros(1,m+1)];
            beq = [Y(j,:)';zeros(m,1);1];
            z = linprog(f,[],[],Aeq,beq,lb);
            Z(j,:) = z;
        end
        Z

    % Output oriented;
    case ('oo')
        Z = zeros(n,n+m+s+1);

        % Objective function of the BCC_oo model: max(0*lambda + epsilon*(s+ + s-) + theta);
        f = -[zeros(1,n), epsilon*ones(1,s+m), 1];

        lblambda = zeros(n,1);                % Lower bounds for (n) lambdas;
        lboutput = zeros(s,1);                % Lower bounds for (s) outputs;
        lbinput = zeros(m,1);                 % Lower bounds for (m) inputs ;
        lb = [lblambda; lboutput; lbinput];   % Lower bounds for lambdas, outputs (s+) and inputs (s-);
        for j=1:n
            Aeq = [-Y', eye(s,s), zeros(s,m), Y(j,:)';
                    X', zeros(m,s), eye(m,m), zeros(m,1);
                    ones(1,n), zeros(1,s+m+1)];
            beq = [zeros(s,1);X(j,:)';1];
            z = linprog(f,[],[],Aeq,beq,lb);
            Z(j,:) = z;
        end
        Z

end