Código para doble recurrencia en R

a=2333
b=1233333333333
alpha=vector()
beta=vector()
e=vector()
z=function(a,b)
{
return(c((a+b)/2,2*a*b/(a+b)))
}

for (i in 1:49){
e=z(a,b)
a=e[1]
b=e[2]
print(e)
alpha[i]=e[1]
beta[i]=e[2]
}

plot(alpha,beta,col=’red’)

Anuncios

Responder

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s