function [code]=orthogoldseq(user,seq,stage,ptap1, ptap2,regi1 ,regi2)                                                           % number of users
                                                 

%**************************** Preparation part *****************************




%********************** Spreading code initialization **********************

%user  = 2;                                                          % number of users
%seq   = 3;                                                          % 1:M-sequence  2:Gold  3:Orthogonal Gold
%stage = 3;                                                          % number of stages
%ptap1 = [1 3];                                                      % position of taps for 1st
%ptap2 = [2 3];                                                      % position of taps for 2nd
%regi1 = [1 1 1];                                                    % initial value of register for 1st
%regi2 = [1 1 1];                                                    % initial value of register for 2nd

%******************** Generation of the spreading code *********************

switch seq
case 1                                                              % M-sequence
    code = mseq(stage,ptap1,regi1,user);
case 2                                                              % Gold sequence
    m1   = mseq(stage,ptap1,regi1);
    m2   = mseq(stage,ptap2,regi2);
    code = goldseq(m1,m2,user);
case 3                                                              % Orthogonal Gold sequence
    m1   = mseq(stage,ptap1,regi1);
    m2   = mseq(stage,ptap2,regi2);
    code = [goldseq(m1,m2,user),zeros(user,1)];
end
code;

clen = length(code);