Thanks works fine, in fact I use sagemath(python,matplotlib)
and as I have a lot of programs to translate in ijulia I often fail.
Without your help I couldn't imagine all these modifications.
Julia list is terrific :) Quick and valid answers most the time :)
My code with your suggestions :
using P
I imagine you are using something like sympy that is doing a mathematica
like plot over omega (the var statements)
In which case you need to change the line
x = linspace(0, 2pi)
to
omega = linspace(0, 2*w) # which seems to be the true range from the python
code, not 2pi
and then you need to c