On Wed, 2016-12-14 at 13:46 +0100, Patrick Dupre wrote:
> Thank very good
>
> However, there is something that I do not understand
>
> I define
> W1=0.2
> z(x) = (x + W1*{0,1})
> w(x) = faddeeva (z (x))
> and
> Voi (x) = voigt (x, W1)
>
> Voi (x) and real (w(x)) are exactly the same
> Now,
> d_w (x) = - 2 * z (x) * w (z (x)) + 2*{0,1}/sqrt(pi) " which is the
> derivative
> D_voi(x,del) = (Voi (x+del) - Voi (x-del)) / 2 / del
>
> now
> plot [-4:4] D_voi(x,0.001)
> and
> plot [-4:4] real (d_w (x))
>
> are different except for W1=0
>
> What am I doing wrong?
> It seems that the issue is with the faddeeva, or when I my calls!
>
> Thank.
>
> =====================================================================
> ======
> Patrick DUPRÉ | | email: pdupre@gmx.
> com
> Laboratoire de Physico-Chimie de l'Atmosphère | |
> Université du Littoral-Côte d'Opale | |
> Tel. (33)-(0)3 28 23 76 12 | | Fax: 03 28 65 82
> 44
> 189A, avenue Maurice Schumann | | 59140 Dunkerque,
> France
> =====================================================================
> ======
> _______________________________________________
> users mailing list -- users@xxxxxxxxxxxxxxxxxxxxxxx
> To unsubscribe send an email to users-leave@xxxxxxxxxxxxxxxxxxxxxxx
Hi, Patrick,
I am not a mathematician, but I have programmed extensively
with complex numbers. W(x)=faddeeva(z(x)) will return a complex,
right?
The magnitude of a complex is its full value, so you would need
the sqrt(real^2+imag^2) would you not?
Regards,
Les H
_______________________________________________
users mailing list -- users@xxxxxxxxxxxxxxxxxxxxxxx
To unsubscribe send an email to users-leave@xxxxxxxxxxxxxxxxxxxxxxx
