Магнитные силовые линии

[Roman007]

Всем доброго дня, вот такая просьба у меня:
как нарисовать магнитную силовую линию dx/dy=By/Bz , B=mu0*I/2pi*sqrt(x*x+y*y) - это магнитная индукция проводника с током. Вот кусок моей программы
function pr=mg(x,y)
B=1/(sqrt(x*x+y*y))
Bx=-B*y/(sqrt(x*x+y*y));
By=B*x/(sqrt(x*x+y*y));
pr=By/Bx
но программа строит только 1/4 окружности (это логично так как тут только одно из 2-х уравнений магнитной силовой линии), вопрос как дописать программу, чтобы она строила всю окружность. Я думаю, что туда надо встроить условия на 0 и перекидывать её на второе уравнение. Подскажите пожалуйста.

[AKM]

Тема перемещена в карантин. Почему это произошло, можно понять, прочитав тему
Что такое карантин, и что нужно делать, чтобы там оказаться
Там же описано, как исправлять ситуацию.

(не оформлены формулы)