Question: How do polarplot ?


 

restart

_local(gamma)

Warning, A new binding for the name `gamma` has been created. The global instance of this name is still accessible using the :- prefix, :-`gamma`.  See ?protect for details.

 

with(plots)

NULL

sol3 := sqrt(2)*sqrt(-tau*gamma)*tanh(x-tau*t^alpha/alpha)*exp(I*gamma*(x+((2*gamma^3*tau-4*gamma*tau+8*tau^2)/(2*gamma^2)-tau*gamma)*t^alpha/((gamma-2*tau)*alpha)))/gamma

NULL

lprint(indets(sol3, name))

{alpha, gamma, t, tau, x}

 

NULL

P :=   [ alpha=1, gamma=-2,  tau=3]

[alpha = 1, gamma = -2, tau = 3]

(1)

PP := convert(sol3, polar)

polar(2^(1/2)*abs(tau*gamma)^(1/2)*exp(-Im(gamma*(x+((1/2)*(2*gamma^3*tau-4*gamma*tau+8*tau^2)/gamma^2-tau*gamma)*t^alpha/((gamma-2*tau)*alpha))))*abs(tanh(x-tau*t^alpha/alpha)/gamma), argument((-tau*gamma)^(1/2)*tanh(x-tau*t^alpha/alpha)*exp(I*gamma*(x+((1/2)*(2*gamma^3*tau-4*gamma*tau+8*tau^2)/gamma^2-tau*gamma)*t^alpha/((gamma-2*tau)*alpha)))/gamma))

(2)

polarplot(sol3, x = -20 .. 20, t = 0 .. 10, axis[radial] = [color = "Blue"])

NULL

Download polar.mw

Please Wait...