
как построить график функции
$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.
Пояснения и комментарииПодынтегральное выражение не определено в нуле, поэтому, чтобы сообщить 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). Поскольку подынтегральное выражение не определено вт=0, начинаем интеграцию ст=1e-9.
Набирается с помощью lualatex
или latex
+ dvips
+ ps2pdf -DNOSAFER
, как минимум два раза:
\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}