
関数をプロットする方法
$f:x\mapsto
\int_0^{x}\frac{1}{\sqrt{1+t^2*ln^2(t)}}\, \textrm{d}t$
TikZで?ありがとう
答え1
上記のリンクで示されているように、この場合、Asymptote が適切な描画ツールです。私は @g.kov のコードを真似しただけです。
// http://asymptote.ualberta.ca/
import graph;
size(300,200,IgnoreAspect);
real f(real t){return (t==0)?1:1/sqrt(1+t^2*(log(t))^2);}
real F(real x){return simpson(f,0,x);}
//real a=20;
real a=1000;
draw(graph(F,0,a,n=200),red+1.5pt);
string noZero(real x) {return (x==0)?"":string(x);}
xaxis(LeftTicks(noZero,Step=a/5));
yaxis(RightTicks(noZero,Step=1));
string s="$F(x)=\displaystyle\int_0^{x}\frac{\textrm{d}t}{\sqrt{1+t^2\ln^2(t)}}$";
label(s,point(E),SW);
shipout(bbox(5mm,invisible));
あなたはLaTeX文書内のasyコード。
説明とコメント積分関数はゼロで定義されていないので、Asymptoteに積分の計算方法を教えるには、次のように書くことができます(このリンク)
\documentclass{article}
\usepackage{amsmath,amssymb}
\begin{document}
Note that $\lim\limits_{t\to 0} t\ln(t)=0$. One may write
\[f(t)=
\begin{cases}
\dfrac{1}{\sqrt{1+t^2\ln^2(t)}},\quad &t\not=0\\
\qquad 1,\quad &t=0.
\end{cases}
\]
\[F(x)=
\displaystyle\int_0^{x} f(t)\;\textrm{d}t, \quad x\geq 0.\]
We see that $F(0)=0$, $F'(x)=f(x)\geq 0$, and $\lim\limits_{t\to +\infty} f(t)=1$.
Therefore $F$ is increasing to infinity \[\lim\limits_{x\to +\infty} F(x)=+\infty.\]
\end{document}
答え2
PGFplotsでプロットし、パッケージpst-ode
(RKF45法)を使用して定積分評価を行う。積分関数は定義されていないため、t=0で積分を開始するt=1e-9.
lualatex
またはlatex
+ dvips
+をps2pdf -DNOSAFER
少なくとも 2 回入力してください。
\documentclass[tikz,border=5pt]{standalone}
\usepackage{pgfplots} \pgfplotsset{compat=newest}
\usepackage{pst-ode}
\pstODEsolve[algebraicAll,saveData]{integral}{t | y[0]}{10^-9}{20}{201}{0}{1/sqrt(1+(t*ln(t))^2)}
\begin{document}
\IfFileExists{integral.dat}{}{dummy text\end{document}}
\begin{tikzpicture}
\begin{axis}[
xlabel=$x$,
ylabel=$F(x)$,
xmin=0,xmax=20,ymin=0
]
\addplot[red,thick] table {integral.dat};
\end{axis}
\end{tikzpicture}
\end{document}