Puede verse claramente que hay una línea vertical en x = 1.
Después de encontrar en la Web varias opciones que no me resultaban satisfactorias, tales como actualizar a la versión 6 o dividir el gráfico en segmentos con dominios donde la función es continua (estaba trabajando con una función con 4 polos), decidí acotar el valor de la función y reemplazar los valores superiores con valores indeterminados.
Esta solución aplicada al problema de mostrar correctamente al gráfico anterior quedaría como:
f[x_]:=1/(x-1);
fp[x_]:=Module[{},
ret=f[x];
If[Abs[ret]>30,Indeterminate,ret]];
Off[Plot::plnr];
Plot[fp[x],{x,-3,3}];
On[Plot::plnr];
Después de encontrar en la Web varias opciones que no me resultaban satisfactorias, tales como actualizar a la versión 6 o dividir el gráfico en segmentos con dominios donde la función es continua (estaba trabajando con una función con 4 polos), decidí acotar el valor de la función y reemplazar los valores superiores con valores indeterminados.
Esta solución aplicada al problema de mostrar correctamente al gráfico anterior quedaría como:
f[x_]:=1/(x-1);
fp[x_]:=Module[{},
ret=f[x];
If[Abs[ret]>30,Indeterminate,ret]];
Off[Plot::plnr];
Plot[fp[x],{x,-3,3}];
On[Plot::plnr];
Los comandos Off[] y On[] desactivan y activan respectivamente la emisión del error Plot::plnr, producido por la función Plot[] al encontrarse con valores no numéricos.
En otros temas, más allá de las funciones elípticas, dos interesantes fotos de objetos de oro sólido (particularmente apropiados con el oro a más de $100000 el kg :-)
- Una pepita de oro de más de 4 kg.
- Un ladrillo de 220 kg de oro (más de $22000000 en la actualidad).
2 comentarios:
Error corregido en la segunda aparición de "$22000000".
Mariano, gracias por publicar esta ayuda para graficar las discontinuas.
Aquí, te copio otro intento para graficar. Podrías decirme cómo hago para desactivar el primer gráfico que me aparece?
Saludos! Laura.
f[x_] := 1/(x - 1);
fp[x_] := Module[{},
ret = f[x];
If[Abs[ret] > 30, Indeterminate, ret]];
Off[Plot::plnr];
xmin := -3;
xmax := 3;
grafica = Plot[fp[x], {x, xmin, xmax},
PlotStyle -> {{Thickness[0.01], Hue[.71]}}];
Off[Plot::plnr];
Clear[x];
Clear[y];
Show[grafica,
AspectRatio -> Automatic,
GridLines -> Automatic,
Axes -> Automatic,
AxesLabel -> {x, y},
AxesOrigin -> Automatic,
AxesStyle -> Automatic,
PlotRange -> {{xmin, xmax}, {-6, 4}}]
On[Plot::plnr];
Publicar un comentario