sábado, 23 de julio de 2011

Gráficas con Sage

Ejemplos

L = [[cos(pi*i/100),sin(pi*i/100)] for i in range(200)]
p = polygon(L, rgbcolor=(1,1,0))
show(p)



L=[[-1+cos(pi*i/100)*(1+cos(pi*i/100)),2*sin(pi*i/100)*(1-cos(pi*i/100))] for i in range(200)]
p = polygon(L, rgbcolor=(1/8,3/4,1/2))
show(p)
L = [[cos(pi*i/100)^3,sin(pi*i/100)] for i in range(200)]
p = line(L, rgbcolor=(1/4,1/8,3/4))
t = text("a bulb", (-1.7, 0.5))
x = text("x axis", (2,-0.2))
y = text("y axis", (0.6,1.3))
g = p+t+x+y
show(g, xmin=-1.5, xmax=2, ymin=-1, ymax=1.3)
L = [[sin(5*pi*i/100)^2*cos(pi*i/100)^3,sin(5*pi*i/100)^2*sin(pi*i/100)] for i in range(200)]
p = polygon(L, rgbcolor=(1/3,1/2,3/5))
show(p)
v=[(i,floor(i)) for i in range(-5,5)]
plot_step_function(v, vertical_lines=False)



g(x)=1/(x-1)
plot(g,(x,-2,2),ymin=-10, ymax=10,detect_poles=True,color='red')+ line([(1,-10), (1,10)],color='green',linestyle='--')



g1(x)=tan(x)

plot(g1,(x,-1.5*pi,1.5*pi),ymin=-10, ymax=10,detect_poles=True,color='red')+ line([(1.5*pi,-10), (1.5*pi,10)],color='green',linestyle='--')+ line([(-1.5*pi,-10), (-1.5*pi,10)],color='green',linestyle='--')+ line([(0.5*pi,-10), (0.5*pi,10)],color='green',linestyle='--')+ line([(-0.5*pi,-10), (-0.5*pi,10)],color='green',linestyle='--')



No hay comentarios:

Publicar un comentario