Сопряженные диаметры эллипса

>

restart;

>	with(plots):

Полуоси

>	a:=4: b:=2:

Кривая

>	F:=(x/a)^2+(y/b)^2=1;

>	Эллипс:=implicitplot(F,x=-a..a,y=-b..b,color=blue):

Диаметр

>

k:=1:

>	for c from -5 to 5 do

>	L[c]:=k*x+c:

Сопряженый диаметр

>	L1[c]:=-(b/a)^2*x/k+c:

>

od:

Хорды и диаметр (c=0)

>	g1:=plot({seq(L[c],c=-4..4)},x=-a..a,-b..b,color=black):

>	g2:=plot({seq(L1[c],c=-4..4)},x=-a..a,-b..b,color=red):

>	display(g1,g2,Эллипс,scaling=constrained);