poisson.mw

> with(plots):
 

> implicitplot( (x^2-y^2-7/2)^2=(1-y^2)*(9/2-2*x)^2 ,x=-4..2.5,y=-1..1, gridrefine=4, crossingrefine=2, resolution=2000, thickness=2,scaling=constrained);
 

Plot_2d
 

> solve( (x^2-y^2-7/2)^2-(1-y^2)*(9/2-2*x)^2=0,y);
 

`+`(`*`(`/`(1, 4), `*`(`^`(`+`(`-`(218), `*`(6, `*`(`^`(`+`(`*`(848, `*`(`^`(x, 2))), `-`(`*`(1872, `*`(x))), `-`(`*`(128, `*`(`^`(x, 3)))), 1377), `/`(1, 2)))), `*`(144, `*`(x)), `-`(`*`(16, `*`(`^`(...
`+`(`*`(`/`(1, 4), `*`(`^`(`+`(`-`(218), `*`(6, `*`(`^`(`+`(`*`(848, `*`(`^`(x, 2))), `-`(`*`(1872, `*`(x))), `-`(`*`(128, `*`(`^`(x, 3)))), 1377), `/`(1, 2)))), `*`(144, `*`(x)), `-`(`*`(16, `*`(`^`(...
`+`(`*`(`/`(1, 4), `*`(`^`(`+`(`-`(218), `*`(6, `*`(`^`(`+`(`*`(848, `*`(`^`(x, 2))), `-`(`*`(1872, `*`(x))), `-`(`*`(128, `*`(`^`(x, 3)))), 1377), `/`(1, 2)))), `*`(144, `*`(x)), `-`(`*`(16, `*`(`^`(...
`+`(`*`(`/`(1, 4), `*`(`^`(`+`(`-`(218), `*`(6, `*`(`^`(`+`(`*`(848, `*`(`^`(x, 2))), `-`(`*`(1872, `*`(x))), `-`(`*`(128, `*`(`^`(x, 3)))), 1377), `/`(1, 2)))), `*`(144, `*`(x)), `-`(`*`(16, `*`(`^`(...
(1)
 

> plot([%],x=-4..2.5,thickness=2,scaling=constrained);
 

Plot_2d
 

>