function f=pred_VAD(azimuth)
[na,nr]=size(azimuth);
if nr==1
    azimuth=azimuth';
end
na=length(azimuth);
az=azimuth;

f(1,:)=ones(1,na);
f(2,:)=sind(az);
f(3,:)=cosd(az);
f(4,:)=sind(2*az);
f(5,:)=cosd(2*az);
f(6,:)=sind(3*az);
f(7,:)=cosd(3*az);

return
end