"Lorenz curve con σ=10, r=28, y b=8/3, es el sistema con las ecuaciones siguientes";

Clear[x,y,z]
σ=10; r=28; b=8/3;
eqs11:={x'[t]==σ*(y[t]-x[t]),y'[t]==r*x[t]-y[t]-x[t]*z[t],z'[t]==x[t]*y[t]-b*z[t],x[0]==5,y[0]==0,z[0]==0}
var:={x,y,z}
lorenzsol=NDSolve[eqs11,var,{t,0,60}]

"Ver las diferencias según la duración de la dinámica, tfinal= 0.5,5, 25 y 60"
lorenzgraph11=ParametricPlot3D[Evaluate[{x[t],y[t],z[t]}/.lorenzsol],{t,0,0.5},PlotPoints->4000,PlotRange->All]
lorenzgraph11=ParametricPlot3D[Evaluate[{x[t],y[t],z[t]}/.lorenzsol],{t,0,5},PlotPoints->4000,PlotRange->All]
lorenzgraph11=ParametricPlot3D[Evaluate[{x[t],y[t],z[t]}/.lorenzsol],{t,0,25},PlotPoints->4000,PlotRange->All]
lorenzgraph11=ParametricPlot3D[Evaluate[{x[t],y[t],z[t]}/.lorenzsol],{t,0,60},PlotPoints->4000,PlotRange->All]

"Proyeccciones en los planos, para tfinal=25"
graphxxyy11=ParametricPlot[Evaluate[{x[t],y[t]}/.lorenzsol],{t,0,25},PlotPoints->4000,PlotRange->All,AxesLabel->{"x","y"},PlotStyle->RGBColor[1,0,0]]
graphxxzz11=ParametricPlot[Evaluate[{x[t],z[t]}/.lorenzsol],{t,0,25},PlotPoints->4000,PlotRange->All,AxesLabel->{"x","z",},PlotStyle->RGBColor[0,1,0]]
graphyyzz11=ParametricPlot[Evaluate[{y[t],z[t]}/.lorenzsol],{t,0,25},PlotPoints->4000,PlotRange->All,AxesLabel->{"y","z"},PlotStyle->RGBColor[0,0,1]]