m0= 156; fs0=22050.0; tm0=[ 0.000000E+00, 1.000000E+00, 2.000000E+00, 3.000000E+00, 4.000000E+00, 5.000000E+00, 6.000000E+00, 7.000000E+00, 8.000000E+00, 9.000000E+00, 1.000000E+01, 1.100000E+01, 1.200000E+01, 1.300000E+01, 1.400000E+01, 1.500000E+01, 1.600000E+01, 1.700000E+01, 1.800000E+01, 1.900000E+01, 2.000000E+01, 2.100000E+01, 2.200000E+01, 2.300000E+01, 2.400000E+01, 2.500000E+01, 2.600000E+01, 2.700000E+01, 2.800000E+01, 2.900000E+01, 3.000000E+01, 3.100000E+01, 3.200000E+01, 3.300000E+01, 3.400000E+01, 3.500000E+01, 3.600000E+01, 3.700000E+01, 3.800000E+01, 3.900000E+01, 4.000000E+01, 4.100000E+01, 4.200000E+01, 4.300000E+01, 4.400000E+01, 4.500000E+01, 4.600000E+01, 4.700000E+01, 4.800000E+01, 4.900000E+01, 5.000000E+01, 5.100000E+01, 5.200000E+01, 5.300000E+01, 5.400000E+01, 5.500000E+01, 5.600000E+01, 5.700000E+01, 5.800000E+01, 5.900000E+01, 6.000000E+01, 6.100000E+01, 6.200000E+01, 6.300000E+01, 6.400000E+01, 6.500000E+01, 6.600000E+01, 6.700000E+01, 6.800000E+01, 6.900000E+01, 7.000000E+01, 7.100000E+01, 7.200000E+01, 7.300000E+01, 7.400000E+01, 7.500000E+01, 7.600000E+01, 7.700000E+01, 7.800000E+01, 7.900000E+01, 8.000000E+01, 8.100000E+01, 8.200000E+01, 8.300000E+01, 8.400000E+01, 8.500000E+01, 8.600000E+01, 8.700000E+01, 8.800000E+01, 8.900000E+01, 9.000000E+01, 9.100000E+01, 9.200000E+01, 9.300000E+01, 9.400000E+01, 9.500000E+01, 9.600000E+01, 9.700000E+01, 9.800000E+01, 9.900000E+01, 1.000000E+02, 1.010000E+02, 1.020000E+02, 1.030000E+02, 1.040000E+02, 1.050000E+02, 1.060000E+02, 1.070000E+02, 1.080000E+02, 1.090000E+02, 1.100000E+02, 1.110000E+02, 1.120000E+02, 1.130000E+02, 1.140000E+02, 1.150000E+02, 1.160000E+02, 1.170000E+02, 1.180000E+02, 1.190000E+02, 1.200000E+02, 1.210000E+02, 1.220000E+02, 1.230000E+02, 1.240000E+02, 1.250000E+02, 1.260000E+02, 1.270000E+02, 1.280000E+02, 1.290000E+02, 1.300000E+02, 1.310000E+02, 1.320000E+02, 1.330000E+02, 1.340000E+02, 1.350000E+02, 1.360000E+02, 1.370000E+02, 1.380000E+02, 1.390000E+02, 1.400000E+02, 1.410000E+02, 1.420000E+02, 1.430000E+02, 1.440000E+02, 1.450000E+02, 1.460000E+02, 1.470000E+02, 1.480000E+02, 1.490000E+02, 1.500000E+02, 1.510000E+02, 1.520000E+02, 1.530000E+02, 1.540000E+02, 1.550000E+02]; ym0=[ 4.040000E+02, 9.300000E+01, -2.890000E+02, -7.390000E+02, -1.244000E+03, -1.787000E+03, -2.348000E+03, -2.900000E+03, -3.416000E+03, -3.865000E+03, -4.217000E+03, -4.444000E+03, -4.521000E+03, -4.431000E+03, -4.160000E+03, -3.706000E+03, -3.075000E+03, -2.281000E+03, -1.350000E+03, -3.140000E+02, 7.850000E+02, 1.902000E+03, 2.988000E+03, 3.994000E+03, 4.870000E+03, 5.574000E+03, 6.070000E+03, 6.332000E+03, 6.343000E+03, 6.100000E+03, 5.612000E+03, 4.899000E+03, 3.993000E+03, 2.936000E+03, 1.778000E+03, 5.720000E+02, -6.250000E+02, -1.759000E+03, -2.778000E+03, -3.638000E+03, -4.303000E+03, -4.746000E+03, -4.955000E+03, -4.927000E+03, -4.673000E+03, -4.213000E+03, -3.580000E+03, -2.810000E+03, -1.949000E+03, -1.043000E+03, -1.400000E+02, 7.140000E+02, 1.482000E+03, 2.126000E+03, 2.622000E+03, 2.950000E+03, 3.103000E+03, 3.082000E+03, 2.896000E+03, 2.563000E+03, 2.106000E+03, 1.555000E+03, 9.420000E+02, 3.000000E+02, -3.370000E+02, -9.380000E+02, -1.475000E+03, -1.925000E+03, -2.268000E+03, -2.492000E+03, -2.591000E+03, -2.564000E+03, -2.416000E+03, -2.159000E+03, -1.805000E+03, -1.375000E+03, -8.880000E+02, -3.680000E+02, 1.610000E+02, 6.780000E+02, 1.160000E+03, 1.589000E+03, 1.946000E+03, 2.218000E+03, 2.396000E+03, 2.474000E+03, 2.451000E+03, 2.328000E+03, 2.114000E+03, 1.818000E+03, 1.454000E+03, 1.040000E+03, 5.930000E+02, 1.340000E+02, -3.170000E+02, -7.440000E+02, -1.127000E+03, -1.452000E+03, -1.706000E+03, -1.881000E+03, -1.971000E+03, -1.976000E+03, -1.899000E+03, -1.747000E+03, -1.530000E+03, -1.261000E+03, -9.550000E+02, -6.270000E+02, -2.950000E+02, 2.600000E+01, 3.230000E+02, 5.830000E+02, 7.960000E+02, 9.540000E+02, 1.054000E+03, 1.096000E+03, 1.080000E+03, 1.012000E+03, 8.980000E+02, 7.460000E+02, 5.650000E+02, 3.660000E+02, 1.570000E+02, -5.000000E+01, -2.500000E+02, -4.330000E+02, -5.930000E+02, -7.250000E+02, -8.250000E+02, -8.900000E+02, -9.180000E+02, -9.080000E+02, -8.610000E+02, -7.770000E+02, -6.590000E+02, -5.100000E+02, -3.340000E+02, -1.380000E+02, 7.200000E+01, 2.880000E+02, 5.030000E+02, 7.050000E+02, 8.850000E+02, 1.034000E+03, 1.143000E+03, 1.205000E+03, 1.215000E+03, 1.172000E+03, 1.078000E+03, 9.360000E+02, 7.560000E+02, 5.520000E+02, 3.370000E+02, 1.300000E+02, -5.000000E+01, -1.860000E+02]; wm0=ones(m0,1);