Hello marcin,
Actually, there is an influence on your Iamp() function. You changed the variable "i" and used it as a loop variable. This way you've lost the imaginary unit. Advice: Do not use variable "i" as a loop variable if using complex numbers.
Numerical result
[MATH=eng]Iamp(225)=1.06*10^{-2}[/MATH]
Symbolical result
[MATH=eng]Iamp(225)—abs({10*(500000+225000*i*π)*(1000+225*i*π)}/{sqrt(2)*(225000*i*π*(500000+225000*i*π)+500000000*(1000+225*i*π))})[/MATH]
Before for loop
[MATH=eng]i—i[/MATH]
[MATH=eng]for(i,range(1,n_fx),el(tab_Iamp,i)←Iamp(el(fx,i)))[/MATH]
After for loop
[MATH=eng]i—101[/MATH]
[MATH=eng]Iamp(225)=7.12*10^{-3}[/MATH]
[MATH=eng]Iamp(225)—abs({10*(500000+22725000*π)*(1000+22725*π)}/{sqrt(2)*(22725000*π*(500000+22725000*π)+500000000*(1000+22725*π))})[/MATH]
Regards,
Radovan
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"