今回描く楕円曲線は
\begin{align}
y^2=x^3-x
\end{align}
解は複素数になるときもあるが、今回は実平面との交点のみを描く。
ルートの中身
\begin{align}
x^3-x
\end{align}
が正になるときだけ描画すればいい
以下コード
x=linspace(-2,2,10000);
count=1;
for i=1:1:length(x)
if((x(1,i)^3-x(1,i))>0)
xp(1,count)=x(1,i);
yp(1,count)=sqrt(x(1,i)^3-x(1,i));
ym(1,count)=-sqrt(x(1,i)^3-x(1,i));
count=count+1;
end
end
figure
plot(xp,yp,"--")
hold on
plot(xp,ym,"--")
コメント