Matlab - Alguien me puede ayudar con un algoritmo en matlab que no me sale bien las graficas

function H=HCircunferencia

x1=input('valor real de Z1:')

y1=input('valor imaginario de Z1:')

p=sqrt(x1^2+y1^2)

x2=input('valor real de Z2:')

y2=input('valor imaginario de Z2:')

r=sqrt(x2^2+y2^2)

t=atan(y2/x2)

n=-2:1:2

Dn=(t+2*pi.*n)

h=[log(r)]/[2.*log(p)]

k=[Dn]/[2.*log(p)]

R=sqrt([[log(r)]^2+[Dn].^2]/[4.*[log(p)]^2])

x=R.*cos(t)+h

y=R.*sin(t)+k

axis equal

grid on

plot(x,y)