Hi,
I would like to draw field lines and equipotential lines of specific field.
I know that plotdf is designed for ODE.
My field is not based on ODE but maybe it can be.
I will describe this field in my non-technical way :
My field is field of complex potential :
http://en.wikibooks.org/wiki/Fractals/Iterations_in_the_complex_plane/Julia_set#Complex_potential_-_Boettcher_coordinate
I can easly compute radial part of this potential :
/* returns floating point number
full definition
*/
GiveFLogPhi(z0,c,e_r,i_max):=
block(
[z:z0,
logphi:log(cabs(z)),
fac:1/2,
i:0],
while i<i_max and cabs(z)<e_r do
(z:z*z+c,
logphi:logphi+fac*log(cabs(1+c/(z*z))),
i:i+1
),
if i=iMax then logphi:0,
return(float(logphi))
)$
or simpler
/* returns floating point number
simple definition
*/
GiveSLogPhi(z0,c,e_r,i_max):=
block(
[z:z0,
logphi,
fac:1/2,
i:0],
while i<i_max and cabs(z)<e_r do
(z:z*z+c,
fac:fac/2,
i:i+1
),
if i=iMax
then logphi:0
else logphi:fac*log(cabs(z)),
return(float(logphi))
)$
Can I plot equipotential lines (contour lines) of this filed using plotdf ?
Can I plot field lines ( gradient vectors orthogonal to
contour lines ) using plotdf ?
(Or other package )
Adam