Apply a digital filter (direct form II transposed): Difference between revisions
Content added Content deleted
m (→version 1: made the displaying of the results into a subroutine to correctly handle the number of (specified) decimal digits.) |
Walterpachl (talk | contribs) (→REXX version 2: increase precision to avoid rounding of s.1=etc.) |
||
Line 465: | Line 465: | ||
{{trans|Julia}} |
{{trans|Julia}} |
||
<lang REXX>/* REXX */ |
<lang REXX>/* REXX */ |
||
Numeric Digits 20 |
|||
acoef = '1.00000000, -2.77555756e-16, 3.33333333e-01, -1.85037171e-17' |
acoef = '1.00000000, -2.77555756e-16, 3.33333333e-01, -1.85037171e-17' |
||
bcoef = '0.16666667, 0.5, 0.5, 0.16666667' |
bcoef = '0.16666667, 0.5, 0.5, 0.16666667' |
||
Line 493: | Line 494: | ||
End |
End |
||
ret.i=temp/a.1 |
ret.i=temp/a.1 |
||
Say format(i,2) format(ret.i,2, |
Say format(i,2) format(ret.i,2,10) |
||
End</lang> |
End</lang> |
||
{{out|output}} |
{{out|output}} |
||
<pre> |
<pre>1 -0.1529739895 |
||
2 -0. |
2 -0.4352578291 |
||
3 -0. |
3 -0.1360433970 |
||
4 0. |
4 0.6975033265 |
||
5 0. |
5 0.6564446925 |
||
6 -0. |
6 -0.4354824533 |
||
7 -1. |
7 -1.0892394612 |
||
8 -0. |
8 -0.5376765496 |
||
9 0. |
9 0.5170499923 |
||
10 1. |
10 1.0522497472 |
||
11 0. |
11 0.9618543004 |
||
12 0. |
12 0.6956900940 |
||
13 0. |
13 0.4243562951 |
||
14 0. |
14 0.1962622318 |
||
15 -0. |
15 -0.0278351245 |
||
16 -0. |
16 -0.2117219155 |
||
17 -0. |
17 -0.1747455622 |
||
18 0. |
18 0.0692584089 |
||
19 0. |
19 0.3854458743 |
||
20 0. |
20 0.6517708388</pre> |
||
=={{header|Sidef}}== |
=={{header|Sidef}}== |