Matlab-uppgift, polära koordinater

Får fel svar på följande uppgift, nån som har lust att titta på min kod? Vore guld i så fall.

"Rita den slutna kurva som i polära koordinater ges av

r(θ) =[2 + sin(3θ)]/sqrt(1 + exp(cos(θ))

För vilka vinklar skär kurvan enhetscirkeln?"

r=@(theta)(2+sin(3*theta))/sqrt(1+exp(cos(theta)));
theta=linspace(0,2*pi)
x=r(theta).*sin(theta);
y=r(theta).*cos(theta);
plot(x,y);
hold on
%% Rita enhetscirkeln
x=cos(theta);
y=sin(theta);
plot(x,y)
hold off

Men jag får en större cirkel utanför enhetscirkeln? Någon som kan förklara?

På den här raden har du missat en punkt.

r = @(theta) (2 + sin(3*theta))./sqrt(1 + exp(cos(theta)));

bör det vara.

Ah, så enkelt var det! Tusen tack.

Ok har stött på problem igen... Har plottat figuren i samma koordinatsystem som enhetscirkeln, men ska nu ta reda på i vilka punkter kurvorna skär varandra. Tänker att man kan sätta ekvationen=1 och sen allt på ett led, använda fzero på något sätt? Nån som har tips?