Despejando tenemos:

Y dando valores a

Con lo que ya tenemos el algoritmo. Lo siguiente es la implementación en sage ploteando juntos el primer punto, los puntos calculados con el método de Euler y la función original (en linea roja) para comparar.
En este caso la función original es . De la que únicamente conocemos la derivada y el punto . Finalmente el algoritmo quedaría:

#Variables
var('h,t,n')
h = 1/10 #paso
puntos = 100 #puntos a calcular
#Condiciones iniciales
f = exp(t)*cos(t) - exp(t)*sin(t) #ecuacion diferencial
x = exp(t)*cos(t) #funcion original para comparar
x0 = 1 #el punto conocido

p = point((0,1))
for i in range(puntos):
    x1 = RR(x0 + h*f(i*h))
    #print x1,RR(x(i*h)),x1-RR(x(i*h))
    p += point((i*h,x1))
    x0 = x1

#graficar la funcion original junto con los puntos hallados
p += plot(x,0,puntos*h,color="red")
p

• No hay artículos relacionados

Supongo que en tiempos de crisis algunos ven positivo cortar las tres libertades de la economía global, sin embargo es en realidad muy peligroso. David analiza muy bien el discurso y concluye que además de la crisis deberíamos temer a quienes la aprovechan para reclamar un mundo más cerrado y autoritario en nombre de una supuesta economía democrática.

Gabilondo dice:

La especulación financiera se hizo global antes de que el mundo construyera alguna estructura global

Pero a mi me parece que las casas no se empiezan por el tejado y que el significado que le da a esa afirmación no tiene pies ni cabeza, porque es una afirmación trivial.

Entiendo poquito de economía, pero tengo alguna idea de los cimientos sobre los que se ha levantado: matemáticas lineales que asumen que el mercado tiene distribución normal. Error, y de los gordos.

No lo sé, y espero investigar pronto cómo ha sido posible desarrollar la ciencia económica si la teoría está equivocada, tal vez no le prestan mucha atención y la práctica empírica se impone a la larga. Pero el caso es que sería muy interesante aplicar la Teoría del caos a la economía.

Quizás me equivoco (una vez y otra) y sean solo impresiones mías pero en un post de Juan Urrutia sobre neurología y neuroeconomía leo entre lineas un poquito de ese orden caótico no lineal que subyace en la economía.

La imparable Neuroeconomía no parece tan útil para la economía como seguramente lo es para la neurología [...] todavía estoy por ver algún experimento que me abra un camino en economía que me parezca interesante recorrer. Y, sin embargo, estoy convencido de que el cogollo de la ciencia económica es una fuente segura de ideas y de caminos a recorrer par los neurólogos

Por el momento no me importa el sentido en el que brotan las ideas en esa fuente, me interesa la implicación y vuelvo a la frase de gabilondo sobre el sistema global no construido. La economía, como todo sistema social es un sistema vivo. Y en vez de ponernos a sangrar al paciente como propone él tal vez deberíamos antes explorar y asimilar estas nuevas matemáticas de la no linealidad y obrar en consecuencia.

• Balance anual 2008

Un día un colega me pegó un telefonazo y me dijo: “¡Ey! No me renuevan el contrato… ¿tienes alguna idea?”. No daría un duro por ninguna de mis ideas, pero tenerlas las tengo a montones, así que sabiendo que él tenía herramientas le propuse montarnos algo por nuestra cuenta. Y ahí está, funcionando, poquito a poco: Telecomunicaciones Gemavi, o la empresa de cómo hacer ñapas con estilo.

Gemavi está superada, fue divertido y enriquecedor, aprendí mucho. Pero no me interesa pasarme la vida grimpando efes y coaxiales… Gemavi me sirve para ir acumulando un pequeñito capital que me está permitiendo desarrollar en los ratos libres otro proyecto, con otro pirado. De momento puedo decir que es un proyecto neoveneciano puro. Es arriesgado, mucho, pero a los osados les sonríe la fortuna y esta vez, de verdad, si que es apasionante :)

• No hay artículos relacionados

En el artículo se habla de la presentación (tras casi diez años de trabajo) del Rosetta disk. Un disco de larga duración que almacena a escala microscópica aproximadamente 13000 páginas en 1500 lenguajes diferentes. Un legado de información para las próximas generaciones.

Toda una proeza.

Sin embargo leo el último artículo de David sobre multilingüísmo y el contraste clama al cielo, el mundo anglo nos ofrece petrificar las lenguas, conservar un fosil para cuando dentro de tres, cinco o diez generaciones toda “la humanidad” hable inglés. Los anglófonos son sordos, ya lo sabemos, pero a mí me parece que además son ciegos: David relata un mundo que ya existe y que cada vez se hará más presente, ciudades-metrópoli como tablero de juego neoveneciano, donde coexistiran en una estructura fractal miles de lenguas distintas y vivas.

La lengua articula el pensamiento, y es la misma lengua la rejilla sobre la cual se cultivan las redes. Por eso si queremos que nuestras redes sean verdaderamente nuestras tenemos que emitir en nuestras lenguas (sean las que sean), solo así nuestra red creará un marco conversacional importante, digno de ser escuchado por otros.

En un mundo transnacional, verdaderamente globalizado, en cada ciudad-metrópoli, un precioso fractal de flores-red se abrirá y su esplendor será como un sonido selvático en el que si escuchamos directamente, sin filtros, lo que oiremos será diversidad, a cada comunidad real emitiendo sus pensamientos, proclamando su modo de vida y defendiendo sus intereses.

• Generador de grafos aleatorio con Sagemaths
• Música libre

Para conseguirlo utilizamos el programa Magpie a través de la extensión Feedimport para Mediawiki, el proceso es muy sencillo:

1. Descargamos Magpie y lo colocamos por ejemplo en ~/lib/Magpie


2. Colocamos el fichero feedimport.php en el directorio de extensiones de Mediawiki:
   ~/mediawiki/extensions


3. Editamos feedimport.php para configurar el directorio en el que se encuentra Magpie:
   # define location of magpie
define('MAGPIE_DIR', '/home/janzo/lib/magpierss/');


4. Y editamos el fichero de configuración de Mediawiki (LocalSettings.php) para incluir la extensión, añadiendo la siguiente linea:
   require_once("$IP/extensions/feedimport.php");


5. Para utilizarlo editamos cualquier página en Mediawiki:
   <rss>http://nonumerable.net/feed|max=5|template=Template:FEED</rss>


6. Y configuramos la apariencia en una plantilla (en este caso FEED) con el formato:
   : [{{{link}}} {{{title}}}]
:: {{{description}}}
:: {{{dc_creator}}} {{{date}}}

Actualización

• La prueba de funcionamiento en este minipost y en la portada de mi contextopedia.
• Como Mediawiki hace los enlaces a su forma particular para que un link aparezca en ambos el formato debe ser la combinación siguiente: <a href="http://micontextopedia/articulo">[[articulo]]</a>

• Lineas de tiempo para mediawiki
• Protocolo ed2k para firefox
• Cómo unir videos con mencoder
• Script para respaldo de bases de datos
• Metaverso Opensim en GNU/Linux
• Disco duro cifrado con Truecrypt
• Instalando Rails con apache en 3 minutos

• ¿Cual es el impacto de la meditación en niños y ancianos? ¿Cual es la dosis óptima de meditación segun edades? Más ampliamente ¿Cómo afecta la meditación en el marco del desarrollo en la vida?
• ¿Qué técnicas se ajustan mejor a cada individuo?
• ¿Cuales son las consecuencias biológicas periféricas? ¿Podemos utilizar la meditación para examinar los mecanismos de interacción mente-cuerpo?

Referencias:
Resumen y artículo en inglés.
Neurociencia y la meditación de atención. Con-ciencia..
Neuroscience of meditation and attention. Mind Hacks.

• No hay artículos relacionados

En mi caso utilizo amule, por lo que primero instalo los útiles de consola amule-utils. Después en firefox abrimos la pestaña de configuración about:config e introducimos dos configuraciones:

network.protocol-handler.external.ed2k boolean true
network.protocol-handler.app.ed2k string /usr/bin/ed2k

Y listo. Si se comete algún error, para borrar una linea de configuración buscaríamos el fichero prefs.js y tras cerrar el navegador haríamos una copia de seguridad y lo editaríamos borrando la entrada erronea.

• Cómo integrar feed en mediawiki
• Cómo unir videos con mencoder
• Script para respaldo de bases de datos
• Lineas de tiempo para mediawiki
• Metaverso Opensim en GNU/Linux
• Disco duro cifrado con Truecrypt
• Instalando Rails con apache en 3 minutos

• No hay artículos relacionados

Un sencillo generador de grafos aleatorios con sagemath

Se llama matriz de adyacencia de un grafo a una matriz cuadrada para representar relaciones entre pares de vértices. Un 1 en la fila i columna j representa una arista entre los vertices i y j.
Es decir, la matriz de adyacencia:

M=Matrix([[1,0,1,1],[0,1,0,1],[1,1,0,1],[0,1,0,1]])
M

[1 0 1 1]
[0 1 0 1]
[1 1 0 1]
[0 1 0 1]

Ejecutada con este código:
M=Matrix([[1,0,1,1],[0,1,0,1],[1,1,0,1],[0,1,0,1]])
G = Graph(M)
G.plot().show()

Representa el siguiente grafo

Modificando ese código se puede construir un sencillo generador de grafos aleatorio.
Para generar un 1 o un 0 aleatoriamente utilizamos: int(round(random()))
y para crear un vector de 10 componentes [i for i in range(10)]. Combinando ambas [int(round(random())) for i in range(10)] se generan vectores como:
[0, 0, 1, 0, 0, 0, 1, 1, 0, 0]