function sy=sigmay(d,esta)%d es la distancia en kilometros y esta es la estabilidad ('A','B','C','D','E' o 'F')
sigy=[3	2	3	5	7	3
0.922861206	0.921091487	0.946394891	0.92646207	0.920802021	0.946394891
0.82718243	0.83931366	0.854322775	0.864652422	0.865874364	0.854322775
217.6880324	158.4327966	106.0659868	67.53860053	49.99807359	35.35532894
241.8159433	167.6728195	117.3561086	74.60280342	55.63813025	39.11870287];
if esta=='A'
  j=1;
elseif esta=='B'
  j=2;
elseif esta=='C'
  j=3;
elseif esta=='D'
  j=4;
elseif esta=='E'
  j=5;
elseif esta=='F'
  j=6;
end
if d<sigy(1,j)
  sy=sigy(4,j)*(d^sigy(2,j));
else
  sy=sigy(5,j)*(d^sigy(3,j));
end
end