code 31253;
procedure GSSWTS(N,ZER,B,C)RESULTS:(W);
value N; integer N;
array ZER,B,C,W;
begin 
    integer J,K; real S; array P[0:N-1];
    for J:=1 step 1 until N do 
    begin 
        ALLORTPOL(N-1,ZER[J],B,C,P);
        S:=0.0;
        for K:=N-1 step -1 until 1 do 
        S:=(S+P[K]**2)/C[K];
        W[J]:=1/(1+S);
    end 
end GSSWTS;
        eop