ECUACIONES EN DIFERENCIAS. EJEMPLO 9.2

PRIMER MÉTODO: SIMULACIÓN ECUACIÓN DIFERENCIAS

(Subscript[y, 0])=100;

For[t=0,t<9,t++,(Subscript[y, t+1])=2(Subscript[y, t])+10]

puntos=ListPlot[{Table[{i,(Subscript[y, i])},{i,9}]}]

f[t_]:=110*(2^t)-10

curva=Plot[f[t],{t,0,8}]

Show[puntos,curva]

SEGUNDO MÉTODO : DIAGRAMA DE COBWEB

f[x_]:=2x+10
g[x_]:=x

NSolve[f[x]==g[x],x]

Cálculo de 25 términos de la órbita para el valor inicial x0=10

NestList[f,10.,25]

orbita=NestList[f,10,25];

a1=ListPlot[Partition[Flatten[Transpose[{orbita,orbita}]],2,1],PlotJoined->True,DisplayFunction->Identity]

a2=Plot[{x,f[x]},{x,0,5*(10^9)},PlotStyle->{RGBColor[1,0,0],RGBColor[0,0,1]},
DisplayFunction->Identity]

Show[a1,a2,AspectRatio->1,DisplayFunction->$DisplayFunction]

grafico=FixedPointList[f,10.,25]

ListPlot[grafico,PlotJoined->False]