function out=oe(t)
a=5;
b=2;
T=0.33;
k=3.4;
l=(b^2)*(cos(t).^2) + (a^2)*(sin(t).^2);
l=sqrt(l);
out = (a + T*b./l).*cos(t)-k;