
% The multiplicative constants for the 41 point FFT

I = sqrt(-1);
u = [
       -1.025000000000000
       -0.160078105935821
       -0.408082396550413
        0.196134539603996
        0.105973928473208
       -0.511082416460775*I
       -0.400250368444121*I
        0.130270939623715*I
        1.034912896131478*I
        0.384121990474013*I
       -0.764933467311073*I
       -0.261915239835351*I
        0.008064188985054*I
        0.317331263843679*I
       -0.512922550835320
       -1.250469324642098
        0.718168584594119
       -0.089594585231855
       -0.852676772254726
        0.621895325279833
        0.416649291521034
        0.948950031569013
       -0.541132646592800
        1.315793252424125
        0.804209603467736
       -0.935334383092475
        0.406453915324166
       -0.135342021986711
       -0.275329135699426
       -0.655788005031077
       -0.524663225406337
        0.478222949845863
       -0.371953329468526
       -0.662829362092971
        0.461451095112636
        0.021261353975737
       -0.010663232232787
        0.180116913004599
        0.090619147360488
        0.291997414340824
       -0.204836913963709
       -1.081062514774166
       -0.306605243045309
       -0.589619949308502
        3.182817569877233
        1.416657149301860
       -2.286592377523421
       -0.615909913440753
        0.280315278913059
        0.642361392271924
        0.726507922121346
        0.484717302569140
        0.064084427097933
       -1.602039461926485
       -0.702996958534537
        1.053237732259411
        0.262645383040132
       -0.286156346626941
       -0.218762239154107
       -2.408397981403094*I
        1.817024222455438*I
        0.835629620705266*I
       -2.013969090497713*I
        3.192884608494123*I
       -0.638684752662991*I
        0.934083608034837*I
       -1.718570235125866*I
        0.117508351809751*I
        5.088944923797523*I
        2.628311859930334*I
       -2.849901872022901*I
        0.792292033624343*I
        0.308955210722705*I
       -0.570702021531776*I
       -2.809745073306398*I
       -2.588155061213830*I
        1.763700805614981*I
       -2.598954393538631*I
       -2.868728025023146*I
        1.981542856173918*I
        0.849931144122134*I
        0.402438100505788*I
        0.037254024727094*I
        0.654665562091807*I
        1.541850730941035*I
       -0.843062634786771*I
       -3.916720937127244*I
        1.488044138823897*I
        0.426779666932540*I
       -1.960675086316631*I
       -1.240040358626114*I
        0.045851280560194*I
        3.535792550754898*I
        1.620968744998461*I
       -1.125878448649219*I
       -6.434082781810076*I
        0.195924171163376*I
        2.375593127962740*I
       -1.482489854945733*I
        3.066675661602289*I
       -1.089027444180384*I
        2.969462209666952*I
        0.397944910540835*I
       -0.930794560798029*I
        1.426730406868352*I
       -1.247627973731084*I
        0.441579348916443*I
        0.318451908570154*I
        0.528647047397621*I
        0.487596716244488*I
       -1.380713039540308*I
       -0.574664414725666*I
        0.205240277821009*I
        3.162559459265169*I
       -1.652534180639667*I
       -0.631204643818903*I
        1.987322088407172*I
       -0.976422124934004*I
        0.296416736051399*I
       -2.234938079394868*I
        0.048800745063703*I
        0.504185048419734*I
        0.672568929006278*I
       -1.412118015546855*I
        0.237154372030080*I
        0.345098910660695*I
       -1.687815436162498*I
        0.829864732856080*I
       -0.079858568180277*I
        1.095105075336497*I
       -0.416453122408476*I
        0.586111993335139*I
        2.058177999377115*I
       -1.211561102545180*I
       -0.584191526346144*I
       -0.465542573951704*I
       -0.262425370485791*I
        0.363023738724250*I
       -0.483593158107685*I
        0.318911178482881*I
];


% The input permutation for the 41 point FFT

ip = [
   1
   2
   19
   38
   11
   17
   15
   7
   27
   18
   20
   33
   3
   37
   34
   21
   39
   29
   13
   12
   35
   41
   24
   5
   32
   26
   28
   36
   16
   25
   23
   10
   40
   6
   9
   22
   4
   14
   30
   31
   8
];


% The output permutation for the 41 point FFT

op = [
   1
   41
   7
   36
   18
   28
   2
   35
   29
   31
   39
   14
   13
   32
   6
   23
   40
   4
   37
   5
   10
   30
   25
   17
   24
   20
   3
   26
   12
   33
   34
   19
   11
   9
   15
   22
   8
   38
   16
   27
   21
];

