t_alpha.m

我认为很不错的语音处理的ｍａｔｌａｂ源代码

        case 29;
            alpha(set_t)=1./2+1./pi.*atan(1./29.*t.*29.^(1./2))+1./pi.*t.*29.^(1./2)./(29+t.^2).*(1+2./3./(1+1./29.*t.^2)+8./15./(1+1./29.*t.^2).^2+16./35./(1+1./29.*t.^2).^3+128./315./(1+1./29.*t.^2).^4+256./693./(1+1./29.*t.^2).^5+1024./3003./(1+1./29.*t.^2).^6+2048./6435./(1+1./29.*t.^2).^7+32768./109395./(1+1./29.*t.^2).^8+65536./230945./(1+1./29.*t.^2).^9+262144./969969./(1+1./29.*t.^2).^10+524288./2028117./(1+1./29.*t.^2).^11+4194304./16900975./(1+1./29.*t.^2).^12+8388608./35102025./(1+1./29.*t.^2).^13);

        case 30;
            alpha(set_t)=1./2+1./2.*t./(30+t.^2).^(1./2).*(1+1./2./(1+1./30.*t.^2)+3./8./(1+1./30.*t.^2).^2+5./16./(1+1./30.*t.^2).^3+35./128./(1+1./30.*t.^2).^4+63./256./(1+1./30.*t.^2).^5+231./1024./(1+1./30.*t.^2).^6+429./2048./(1+1./30.*t.^2).^7+6435./32768./(1+1./30.*t.^2).^8+12155./65536./(1+1./30.*t.^2).^9+46189./262144./(1+1./30.*t.^2).^10+88179./524288./(1+1./30.*t.^2).^11+676039./4194304./(1+1./30.*t.^2).^12+1300075./8388608./(1+1./30.*t.^2).^13+5014575./33554432./(1+1./30.*t.^2).^14);

        case 31;
            alpha(set_t)=1./2+1./pi.*atan(1./31.*t.*31.^(1./2))+1./pi.*t.*31.^(1./2)./(31+t.^2).*(1+2./3./(1+1./31.*t.^2)+8./15./(1+1./31.*t.^2).^2+16./35./(1+1./31.*t.^2).^3+128./315./(1+1./31.*t.^2).^4+256./693./(1+1./31.*t.^2).^5+1024./3003./(1+1./31.*t.^2).^6+2048./6435./(1+1./31.*t.^2).^7+32768./109395./(1+1./31.*t.^2).^8+65536./230945./(1+1./31.*t.^2).^9+262144./969969./(1+1./31.*t.^2).^10+524288./2028117./(1+1./31.*t.^2).^11+4194304./16900975./(1+1./31.*t.^2).^12+8388608./35102025./(1+1./31.*t.^2).^13+33554432./145422675./(1+1./31.*t.^2).^14);

        case 32;
            alpha(set_t)=1./2+1./2.*t./(32+t.^2).^(1./2).*(1+1./2./(1+1./32.*t.^2)+3./8./(1+1./32.*t.^2).^2+5./16./(1+1./32.*t.^2).^3+35./128./(1+1./32.*t.^2).^4+63./256./(1+1./32.*t.^2).^5+231./1024./(1+1./32.*t.^2).^6+429./2048./(1+1./32.*t.^2).^7+6435./32768./(1+1./32.*t.^2).^8+12155./65536./(1+1./32.*t.^2).^9+46189./262144./(1+1./32.*t.^2).^10+88179./524288./(1+1./32.*t.^2).^11+676039./4194304./(1+1./32.*t.^2).^12+1300075./8388608./(1+1./32.*t.^2).^13+5014575./33554432./(1+1./32.*t.^2).^14+9694845./67108864./(1+1./32.*t.^2).^15);

        case 33;
            alpha(set_t)=1./2+1./pi.*atan(1./33.*t.*33.^(1./2))+1./pi.*t.*33.^(1./2)./(33+t.^2).*(1+2./3./(1+1./33.*t.^2)+8./15./(1+1./33.*t.^2).^2+16./35./(1+1./33.*t.^2).^3+128./315./(1+1./33.*t.^2).^4+256./693./(1+1./33.*t.^2).^5+1024./3003./(1+1./33.*t.^2).^6+2048./6435./(1+1./33.*t.^2).^7+32768./109395./(1+1./33.*t.^2).^8+65536./230945./(1+1./33.*t.^2).^9+262144./969969./(1+1./33.*t.^2).^10+524288./2028117./(1+1./33.*t.^2).^11+4194304./16900975./(1+1./33.*t.^2).^12+8388608./35102025./(1+1./33.*t.^2).^13+33554432./145422675./(1+1./33.*t.^2).^14+67108864./300540195./(1+1./33.*t.^2).^15);

        case 34;
            alpha(set_t)=1./2+1./2.*t./(34+t.^2).^(1./2).*(1+1./2./(1+1./34.*t.^2)+3./8./(1+1./34.*t.^2).^2+5./16./(1+1./34.*t.^2).^3+35./128./(1+1./34.*t.^2).^4+63./256./(1+1./34.*t.^2).^5+231./1024./(1+1./34.*t.^2).^6+429./2048./(1+1./34.*t.^2).^7+6435./32768./(1+1./34.*t.^2).^8+12155./65536./(1+1./34.*t.^2).^9+46189./262144./(1+1./34.*t.^2).^10+88179./524288./(1+1./34.*t.^2).^11+676039./4194304./(1+1./34.*t.^2).^12+1300075./8388608./(1+1./34.*t.^2).^13+5014575./33554432./(1+1./34.*t.^2).^14+9694845./67108864./(1+1./34.*t.^2).^15+300540195./2147483648./(1+1./34.*t.^2).^16);

        case 35;
            alpha(set_t)=1./2+1./pi.*atan(1./35.*t.*35.^(1./2))+1./pi.*t.*35.^(1./2)./(35+t.^2).*(1+2./3./(1+1./35.*t.^2)+8./15./(1+1./35.*t.^2).^2+16./35./(1+1./35.*t.^2).^3+128./315./(1+1./35.*t.^2).^4+256./693./(1+1./35.*t.^2).^5+1024./3003./(1+1./35.*t.^2).^6+2048./6435./(1+1./35.*t.^2).^7+32768./109395./(1+1./35.*t.^2).^8+65536./230945./(1+1./35.*t.^2).^9+262144./969969./(1+1./35.*t.^2).^10+524288./2028117./(1+1./35.*t.^2).^11+4194304./16900975./(1+1./35.*t.^2).^12+8388608./35102025./(1+1./35.*t.^2).^13+33554432./145422675./(1+1./35.*t.^2).^14+67108864./300540195./(1+1./35.*t.^2).^15+2147483648./9917826435./(1+1./35.*t.^2).^16);

        case 36;
            alpha(set_t)=1./2+1./2.*t./(36+t.^2).^(1./2).*(1+1./2./(1+1./36.*t.^2)+3./8./(1+1./36.*t.^2).^2+5./16./(1+1./36.*t.^2).^3+35./128./(1+1./36.*t.^2).^4+63./256./(1+1./36.*t.^2).^5+231./1024./(1+1./36.*t.^2).^6+429./2048./(1+1./36.*t.^2).^7+6435./32768./(1+1./36.*t.^2).^8+12155./65536./(1+1./36.*t.^2).^9+46189./262144./(1+1./36.*t.^2).^10+88179./524288./(1+1./36.*t.^2).^11+676039./4194304./(1+1./36.*t.^2).^12+1300075./8388608./(1+1./36.*t.^2).^13+5014575./33554432./(1+1./36.*t.^2).^14+9694845./67108864./(1+1./36.*t.^2).^15+300540195./2147483648./(1+1./36.*t.^2).^16+583401555./4294967296./(1+1./36.*t.^2).^17);

        case 37;
            alpha(set_t)=1./2+1./pi.*atan(1./37.*t.*37.^(1./2))+1./pi.*t.*37.^(1./2)./(37+t.^2).*(1+2./3./(1+1./37.*t.^2)+8./15./(1+1./37.*t.^2).^2+16./35./(1+1./37.*t.^2).^3+128./315./(1+1./37.*t.^2).^4+256./693./(1+1./37.*t.^2).^5+1024./3003./(1+1./37.*t.^2).^6+2048./6435./(1+1./37.*t.^2).^7+32768./109395./(1+1./37.*t.^2).^8+65536./230945./(1+1./37.*t.^2).^9+262144./969969./(1+1./37.*t.^2).^10+524288./2028117./(1+1./37.*t.^2).^11+4194304./16900975./(1+1./37.*t.^2).^12+8388608./35102025./(1+1./37.*t.^2).^13+33554432./145422675./(1+1./37.*t.^2).^14+67108864./300540195./(1+1./37.*t.^2).^15+2147483648./9917826435./(1+1./37.*t.^2).^16+4294967296./20419054425./(1+1./37.*t.^2).^17);

        case 38;
            alpha(set_t)=1./2+1./2.*t./(38+t.^2).^(1./2).*(1+1./2./(1+1./38.*t.^2)+3./8./(1+1./38.*t.^2).^2+5./16./(1+1./38.*t.^2).^3+35./128./(1+1./38.*t.^2).^4+63./256./(1+1./38.*t.^2).^5+231./1024./(1+1./38.*t.^2).^6+429./2048./(1+1./38.*t.^2).^7+6435./32768./(1+1./38.*t.^2).^8+12155./65536./(1+1./38.*t.^2).^9+46189./262144./(1+1./38.*t.^2).^10+88179./524288./(1+1./38.*t.^2).^11+676039./4194304./(1+1./38.*t.^2).^12+1300075./8388608./(1+1./38.*t.^2).^13+5014575./33554432./(1+1./38.*t.^2).^14+9694845./67108864./(1+1./38.*t.^2).^15+300540195./2147483648./(1+1./38.*t.^2).^16+583401555./4294967296./(1+1./38.*t.^2).^17+2268783825./17179869184./(1+1./38.*t.^2).^18);

        case 39;
            alpha(set_t)=1./2+1./pi.*atan(1./39.*t.*39.^(1./2))+1./pi.*t.*39.^(1./2)./(39+t.^2).*(1+2./3./(1+1./39.*t.^2)+8./15./(1+1./39.*t.^2).^2+16./35./(1+1./39.*t.^2).^3+128./315./(1+1./39.*t.^2).^4+256./693./(1+1./39.*t.^2).^5+1024./3003./(1+1./39.*t.^2).^6+2048./6435./(1+1./39.*t.^2).^7+32768./109395./(1+1./39.*t.^2).^8+65536./230945./(1+1./39.*t.^2).^9+262144./969969./(1+1./39.*t.^2).^10+524288./2028117./(1+1./39.*t.^2).^11+4194304./16900975./(1+1./39.*t.^2).^12+8388608./35102025./(1+1./39.*t.^2).^13+33554432./145422675./(1+1./39.*t.^2).^14+67108864./300540195./(1+1./39.*t.^2).^15+2147483648./9917826435./(1+1./39.*t.^2).^16+4294967296./20419054425./(1+1./39.*t.^2).^17+17179869184./83945001525./(1+1./39.*t.^2).^18);

        case 40;
            alpha(set_t)=1./2+1./2.*t./(40+t.^2).^(1./2).*(1+1./2./(1+1./40.*t.^2)+3./8./(1+1./40.*t.^2).^2+5./16./(1+1./40.*t.^2).^3+35./128./(1+1./40.*t.^2).^4+63./256./(1+1./40.*t.^2).^5+231./1024./(1+1./40.*t.^2).^6+429./2048./(1+1./40.*t.^2).^7+6435./32768./(1+1./40.*t.^2).^8+12155./65536./(1+1./40.*t.^2).^9+46189./262144./(1+1./40.*t.^2).^10+88179./524288./(1+1./40.*t.^2).^11+676039./4194304./(1+1./40.*t.^2).^12+1300075./8388608./(1+1./40.*t.^2).^13+5014575./33554432./(1+1./40.*t.^2).^14+9694845./67108864./(1+1./40.*t.^2).^15+300540195./2147483648./(1+1./40.*t.^2).^16+583401555./4294967296./(1+1./40.*t.^2).^17+2268783825./17179869184./(1+1./40.*t.^2).^18+4418157975./34359738368./(1+1./40.*t.^2).^19);

        case 41;
            alpha(set_t)=1./2+1./pi.*atan(1./41.*t.*41.^(1./2))+1./pi.*t.*41.^(1./2)./(41+t.^2).*(1+2./3./(1+1./41.*t.^2)+8./15./(1+1./41.*t.^2).^2+16./35./(1+1./41.*t.^2).^3+128./315./(1+1./41.*t.^2).^4+256./693./(1+1./41.*t.^2).^5+1024./3003./(1+1./41.*t.^2).^6+2048./6435./(1+1./41.*t.^2).^7+32768./109395./(1+1./41.*t.^2).^8+65536./230945./(1+1./41.*t.^2).^9+262144./969969./(1+1./41.*t.^2).^10+524288./2028117./(1+1./41.*t.^2).^11+4194304./16900975./(1+1./41.*t.^2).^12+8388608./35102025./(1+1./41.*t.^2).^13+33554432./145422675./(1+1./41.*t.^2).^14+67108864./300540195./(1+1./41.*t.^2).^15+2147483648./9917826435./(1+1./41.*t.^2).^16+4294967296./20419054425./(1+1./41.*t.^2).^17+17179869184./83945001525./(1+1./41.*t.^2).^18+34359738368./172308161025./(1+1./41.*t.^2).^19);

        case 42;
            alpha(set_t)=1./2+1./2.*t./(42+t.^2).^(1./2).*(1+1./2./(1+1./42.*t.^2)+3./8./(1+1./42.*t.^2).^2+5./16./(1+1./42.*t.^2).^3+35./128./(1+1./42.*t.^2).^4+63./256./(1+1./42.*t.^2).^5+231./1024./(1+1./42.*t.^2).^6+429./2048./(1+1./42.*t.^2).^7+6435./32768./(1+1./42.*t.^2).^8+12155./65536./(1+1./42.*t.^2).^9+46189./262144./(1+1./42.*t.^2).^10+88179./524288./(1+1./42.*t.^2).^11+676039./4194304./(1+1./42.*t.^2).^12+1300075./8388608./(1+1./42.*t.^2).^13+5014575./33554432./(1+1./42.*t.^2).^14+9694845./67108864./(1+1./42.*t.^2).^15+300540195./2147483648./(1+1./42.*t.^2).^16+583401555./4294967296./(1+1./42.*t.^2).^17+2268783825./17179869184./(1+1./42.*t.^2).^18+4418157975./34359738368./(1+1./42.*t.^2).^19+34461632205./274877906944./(1+1./42.*t.^2).^20);

        case 43;
            alpha(set_t)=1./2+1./pi.*atan(1./43.*t.*43.^(1./2))+1./pi.*t.*43.^(1./2)./(43+t.^2).*(1+2./3./(1+1./43.*t.^2)+8./15./(1+1./43.*t.^2).^2+16./35./(1+1./43.*t.^2).^3+128./315./(1+1./43.*t.^2).^4+256./693./(1+1./43.*t.^2).^5+1024./3003./(1+1./43.*t.^2).^6+2048./6435./(1+1./43.*t.^2).^7+32768./109395./(1+1./43.*t.^2).^8+65536./230945./(1+1./43.*t.^2).^9+262144./969969./(1+1./43.*t.^2).^10+524288./2028117./(1+1./43.*t.^2).^11+4194304./16900975./(1+1./43.*t.^2).^12+8388608./35102025./(1+1./43.*t.^2).^13+33554432./145422675./(1+1./43.*t.^2).^14+67108864./300540195./(1+1./43.*t.^2).^15+2147483648./9917826435./(1+1./43.*t.^2).^16+4294967296./20419054425./(1+1./43.*t.^2).^17+17179869184./83945001525./(1+1./43.*t.^2).^18+34359738368./172308161025./(1+1./43.*t.^2).^19+274877906944./1412926920405./(1+1./43.*t.^2).^20);

        case 44;
            alpha(set_t)=1./2+1./2.*t./(44+t.^2).^(1./2).*(1+1./2./(1+1./44.*t.^2)+3./8./(1+1./44.*t.^2).^2+5./16./(1+1./44.*t.^2).^3+35./128./(1+1./44.*t.^2).^4+63./256./(1+1./44.*t.^2).^5+231./1024./(1+1./44.*t.^2).^6+429./2048./(1+1./44.*t.^2).^7+6435./32768./(1+1./44.*t.^2).^8+12155./65536./(1+1./44.*t.^2).^9+46189./262144./(1+1./44.*t.^2).^10+88179./524288./(1+1./44.*t.^2).^11+676039./4194304./(1+1./44.*t.^2).^12+1300075./8388608./(1+1./44.*t.^2).^13+5014575./33554432./(1+1./44.*t.^2).^14+9694845./67108864./(1+1./44.*t.^2).^15+300540195./2147483648./(1+1./44.*t.^2).^16+583401555./4294967296./(1+1./44.*t.^2).^17+2268783825./17179869184./(1+1./44.*t.^2).^18+4418157975./34359738368./(1+1./44.*t.^2).^19+34461632205./274877906944./(1+1./44.*t.^2).^20+67282234305./549755813888./(1+1./44.*t.^2).^21);

        case 45;
            alpha(set_t)=1./2+1./pi.*atan(1./45.*t.*45.^(1./2))+1./pi.*t.*45.^(1./2)./(45+t.^2).*(1+2./3./(1+1./45.*t.^2)+8./15./(1+1./45.*t.^2).^2+16./35./(1+1./45.*t.^2).^3+128./315./(1+1./45.*t.^2).^4+256./693./(1+1./45.*t.^2).^5+1024./3003./(1+1./45.*t.^2).^6+2048./6435./(1+1./45.*t.^2).^7+32768./109395./(1+1./45.*t.^2).^8+65536./230945./(1+1./45.*t.^2).^9+262144./969969./(1+1./45.*t.^2).^10+524288./2028117./(1+1./45.*t.^2).^11+4194304./16900975./(1+1./45.*t.^2).^12+8388608./35102025./(1+1./45.*t.^2).^13+33554432./145422675./(1+1./45.*t.^2).^14+67108864./300540195./(1+1./45.*t.^2).^15+2147483648./9917826435./(1+1./45.*t.^2).^16+4294967296./20419054425./(1+1./45.*t.^2).^17+17179869184./83945001525./(1+1./45.*t.^2).^18+34359738368./172308161025./(1+1./45.*t.^2).^19+274877906944./1412926920405./(1+1./45.*t.^2).^20+549755813888./2893136075115./(1+1./45.*t.^2).^21);

        case 46;
            alpha(set_t)=1./2+1./2.*t./(46+t.^2).^(1./2).*(1+1./2./(1+1./46.*t.^2)+3./8./(1+1./46.*t.^2).^2+5./16./(1+1./46.*t.^2).^3+35./128./(1+1./46.*t.^2).^4+63./256./(1+1./46.*t.^2).^5+231./1024./(1+1./46.*t.^2).^6+429./2048./(1+1./46.*t.^2).^7+6435./32768./(1+1./46.*t.^2).^8+12155./65536./(1+1./46.*t.^2).^9+46189./262144./(1+1./46.*t.^2).^10+88179./524288./(1+1./46.*t.^2).^11+676039./4194304./(1+1./46.*t.^2).^12+1300075./8388608./(1+1./46.*t.^2).^13+5014575./33554432./(1+1./46.*t.^2).^14+9694845./67108864./(1+1./46.*t.^2).^15+300540195./2147483648./(1+1./46.*t.^2).^16+583401555./4294967296./(1+1./46.*t.^2).^17+2268783825./17179869184./(1+1./46.*t.^2).^18+4418157975./34359738368./(1+1./46.*t.^2).^19+34461632205./274877906944./(1+1./46.*t.^2).^20+67282234305./549755813888./(1+1./46.*t.^2).^21+263012370465./2199023255552./(1+1./46.*t.^2).^22);

        case 47;
            alpha(set_t)=1./2+1./pi.*atan(1./47.*t.*47.^(1./2))+1./pi.*t.*47.^(1./2)./(47+t.^2).*(1+2./3./(1+1./47.*t.^2)+8./15./(1+1./47.*t.^2).^2+16./35./(1+1./47.*t.^2).^3+128./315./(1+1./47.*t.^2).^4+256./693./(1+1./47.*t.^2).^5+1024./3003./(1+1./47.*t.^2).^6+2048./6435./(1+1./47.*t.^2).^7+32768./109395./(1+1./47.*t.^2).^8+65536./230945./(1+1./47.*t.^2).^9+262144./969969./(1+1./47.*t.^2).^10+524288./2028117./(1+1./47.*t.^2).^11+4194304./16900975./(1+1./47.*t.^2).^12+8388608./35102025./(1+1./47.*t.^2).^13+33554432./145422675./(1+1./47.*t.^2).^14+67108864./300540195./(1+1./47.*t.^2).^15+2147483648./9917826435./(1+1./47.*t.^2).^16+4294967296./20419054425./(1+1./47.*t.^2).^17+17179869184./83945001525./(1+1./47.*t.^2).^18+34359738368./172308161025./(1+1./47.*t.^2).^19+274877906944./1412926920405./(1+1./47.*t.^2).^20+549755813888./2893136075115./(1+1./47.*t.^2).^21+2199023255552./11835556670925./(1+1./47.*t.^2).^22);

        case 48;
            alpha(set_t)=1./2+1./2.*t./(48+t.^2).^(1./2).*(1+300540195./2147483648./(1+1./48.*t.^2).^16+583401555./4294967296./(1+1./48.*t.^2).^17+2268783825./17179869184./(1+1./48.*t.^2).^18+4418157975./34359738368./(1+1./48.*t.^2).^19+34461632205./274877906944./(1+1./48.*t.^2).^20+67282234305./549755813888./(1+1./48.*t.^2).^21+263012370465./2199023255552./(1+1./48.*t.^2).^22+514589420475./4398046511104./(1+1./48.*t.^2).^23+1./2./(1+1./48.*t.^2)+3./8./(1+1./48.*t.^2).^2+5./16./(1+1./48.*t.^2).^3+35./128./(1+1./48.*t.^2).^4+63./256./(1+1./48.*t.^2).^5+231./1024./(1+1./48.*t.^2).^6+429./2048./(1+1./48.*t.^2).^7+12155./65536./(1+1./48.*t.^2).^9+6435./32768./(1+1./48.*t.^2).^8+46189./262144./(1+1./48.*t.^2).^10+88179./524288./(1+1./48.*t.^2).^11+1300075./8388608./(1+1./48.*t.^2).^13+676039./4194304./(1+1./48.*t.^2).^12+5014575./33554432./(1+1./48.*t.^2).^14+9694845./67108864./(1+1./48.*t.^2).^15);

        case 49;
            alpha(set_t)=1./2+1./pi.*atan(1./49.*t.*49.^(1./2))+1./pi.*t.*49.^(1./2)./(49+t.^2).*(1+274877906944./1412926920405./(1+1./49.*t.^2).^20+549755813888./2893136075115./(1+1./49.*t.^2).^21+262144./969969./(1+1./49.*t.^2).^10+524288./2028117./(1+1./49.*t.^2).^11+8./15./(1+1./49.*t.^2).^2+16./35./(1+1./49.*t.^2).^3+2199023255552./11835556670925./(1+1./49.*t.^2).^22+4194304./16900975./(1+1./49.*t.^2).^12+8388608./35102025./(1+1./49.*t.^2).^13+4398046511104./24185702762325./(1+1./49.*t.^2).^23+33554432./145422675./(1+1./49.*t.^2).^14+128./315./(1+1./49.*t.^2).^4+67108864./300540195./(1+1./49.*t.^2).^15+2147483648./9917826435./(1+1./49.*t.^2).^16+4294967296./20419054425./(1+1./49.*t.^2).^17+256./693./(1+1./49.*t.^2).^5+1024./3003./(1+1./49.*t.^2).^6+2048./6435./(1+1./49.*t.^2).^7+17179869184./83945001525./(1+1./49.*t.^2).^18+34359738368./172308161025./(1+1./49.*t.^2).^19+32768./109395./(1+1./49.*t.^2).^8+65536./230945./(1+1./49.*t.^2).^9+2./3./(1+1./49.*t.^2));

        case 50;
            alpha(set_t)=1./2+1./2.*t./(50+t.^2).^(1./2).*(1+1./2./(1+1./50.*t.^2)+3./8./(1+1./50.*t.^2).^2+5./16./(1+1./50.*t.^2).^3+63./256./(1+1./50.*t.^2).^5+35./128./(1+1./50.*t.^2).^4+231./1024./(1+1./50.*t.^2).^6+6435./32768./(1+1./50.*t.^2).^8+429./2048./(1+1./50.*t.^2).^7+46189./262144./(1+1./50.*t.^2).^10+12155./65536./(1+1./50.*t.^2).^9+88179./524288./(1+1./50.*t.^2).^11+676039./4194304./(1+1./50.*t.^2).^12+1300075./8388608./(1+1./50.*t.^2).^13+9694845./67108864./(1+1./50.*t.^2).^15+5014575./33554432./(1+1./50.*t.^2).^14+300540195./2147483648./(1+1./50.*t.^2).^16+2268783825./17179869184./(1+1./50.*t.^2).^18+583401555./4294967296./(1+1./50.*t.^2).^17+4418157975./34359738368./(1+1./50.*t.^2).^19+34461632205./274877906944./(1+1./50.*t.^2).^20+67282234305./549755813888./(1+1./50.*t.^2).^21+263012370465./2199023255552./(1+1./50.*t.^2).^22+514589420475./4398046511104./(1+1./50.*t.^2).^23+8061900920775./70368744177664./(1+1./50.*t.^2).^24);

        case 51;
            alpha(set_t)=1./2+1./pi.*atan(1./51.*t.*51.^(1./2))+1./pi.*t.*51.^(1./2)./(51+t.^2).*(1+2./3./(1+1./51.*t.^2)+8./15./(1+1./51.*t.^2).^2+16./35./(1+1./51.*t.^2).^3+128./315./(1+1./51.*t.^2).^4+256./693./(1+1./51.*t.^2).^5+1024./3003./(1+1./51.*t.^2).^6+2048./6435./(1+1./51.*t.^2).^7+32768./109395./(1+1./51.*t.^2).^8+65536./230945./(1+1./51.*t.^2).^9+262144./969969./(1+1./51.*t.^2).^10+524288./2028117./(1+1./51.*t.^2).^11+4194304./16900975./(1+1./51.*t.^2).^12+8388608./35102025./(1+1./51.*t.^2).^13+33554432./145422675./(1+1./51.*t.^2).^14+67108864./300540195./(1+1./51.*t.^2).^15+2147483648./9917826435./(1+1./51.*t.^2).^16+4294967296./20419054425./(1+1./51.*t.^2).^17+17179869184./83945001525./(1+1./51.*t.^2).^18+34359738368./172308161025./(1+1./51.*t.^2).^19+274877906944./1412926920405./(1+1./51.*t.^2).^20+549755813888./2893136075115./(1+1./51.*t.^2).^21+2199023255552./11835556670925./(1+1./51.*t.^2).^22+4398046511104./24185702762325./(1+1./51.*t.^2).^23+70368744177664./395033145117975./(1+1./51.*t.^2).^24);

        case 52;
            alpha(set_t)=1./2+1./2.*t./(52+t.^2).^(1./2).*(1+514589420475./4398046511104./(1+1./52.*t.^2).^23+8061900920775./70368744177664./(1+1./52.*t.^2).^24+15801325804719./140737488355328./(1+1./52.*t.^2).^25+4418157975./34359738368./(1+1./52.*t.^2).^19+1./2./(1+1./52.*t.^2)+3./8./(1+1./52.*t.^2).^2+5./16./(1+1./52.*t.^2).^3+35./128./(1+1./52.*t.^2).^4+63./256./(1+1./52.*t.^2).^5+231./1024./(1+1./52.*t.^2).^6+429./2048./(1+1./52.*t.^2).^7+34461632205./274877906944./(1+1./52.*t.^2).^20+6435./32768./(1+1./52.*t.^2).^8+12155./65536./(1+1./52.*t.^2).^9+46189./262144./(1+1./52.*t.^2).^10+88179./524288./(1+1./52.*t.^2).^11+1300075./8388608./(1+1./52.*t.^2).^13+676039./4194304./(1+1./52.*t.^2).^12+5014575./33554432./(1+1./52.*t.^2).^14+67282234305./549755813888./(1+1./52.*t.^2).^21+263012370465./2199023255552./(1+1./52.*t.^2).^22+9694845./67108864./(1+1./52.*t.^2).^15+300540195./2147483648./(1+1./52.*t.^2).^16+583401555./4294967296./(1+1./52.*t.^2).^17+2268783825./17179869184./(1+1./52.*t.^2).^18);

        case 53;
            alpha(set_t)=1./2+1./pi.*atan(1./53.*t.*53.^(1./2))+1./pi.*t.*53.^(1./2)./(53+t.^2).*(1+17179869184./83945001525./(1+1./53.*t.^2).^18+70368744177664./395033145117975./(1+1./53.*t.^2).^24+140737488355328./805867616040669./(1+1./53.*t.^2).^25+4294967296./20419054425./(1+1./53.*t.^2).^17+34359738368./172308161025./(1+1./53.*t.^2).^19+274877906944./1412926920405./(1+1./53.*t.^2).^20+2./3./(1+1./53.*t.^2)+8./15./(1+1./53.*t.^2).^2+16./35./(1+1./53.*t.^2).^3+128./315./(1+1./53.*t.^2).^4+256./693./(1+1./53.*t.^2).^5+1024./3003./(1+1./53.*t.^2).^6+549755813888./2893136075115./(1+1./53.*t.^2).^21+65536./230945./(1+1./53.*t.^2).^9+2048./6435./(1+1./53.*t.^2).^7+32768./109395./(1+1./53.*t.^2).^8+262144./969969./(1+1./53.*t.^2).^10+524288./2028117./(1+1./53.*t.^2).^11+4194304./16900975./(1+1./53.*t.^2).^12+2199023255552./11835556670925./(1+1./53.*t.^2).^22+4398046511104./24185702762325./(1+1./53.*t.^2).^23+33554432./145422675./(1+1./53.*t.^2).^14+8388608./35102025./(1+1./53.*t.^2).^13+2147483648./9917826435./(1+1./53.*t.^2).^16+67108864./300540195./(1+1./53.*t.^2).^15);