> with(plots);with(plottools);

> k := t-> 6/(9*sin(t)^2+4*cos(t)^2)^(3/2);

> OCE:=proc(TT)
    s:= TT ->if(evalf(Pi/2)<=TT and TT<evalf(3*Pi/2))  then -1 else 1 fi;
    display(circle([3*cos(TT) -s(TT)*(1/k(TT))*cos(arctan(3/2*tan(TT))) ,2*sin(TT)-s(TT)*(1/k(TT))*sin(arctan(3/2*tan(TT)))],1/k(TT),color=green, thickness=2,numpoints=1000));
end:

> animate(OCE,[TT], TT=0..2*Pi,background= plot([3*cos(t),2*sin(t),t=0..2*Pi],view=[-5..5,-8..8],thickness=2) ,scaling=constrained,numpoints=10000);

[Plot]

>

>