run; 
runanalysis;
expansion = getresult("Exp","expansion for out");
f=expansion.f;
SP = sourcepower(f);

# these quantites should be equal
?getattribute(expansion,"T_forward");
?abs(expansion.a)^2*expansion.N / ( (1)^2 * SP );