?function?example11_2
clear
clc
????function?dx=lorenz（tx）
??????
%????????dx/dt?=?-a（x-y）
%????????dy/dt?=?x（r-z）-y
%????????dz/dt?=?xy-bz
%????a=y（4）
%????r=y（5）
%????b=y（6）
dx=zeros（61）;
dx（1）=-x（4）*（x（1）-x（2））;
dx（2）=x（1）*（x（5）-x（3））-x（2）;
dx（3）=x（1）*x（2）-x（6）*x（3）;
dx（4）=0;
dx（5）=0;
dx（6）=0;
%?dx（11）=-x（4）*x（1）+x（2）*（x（3）+x（5））;
%?dx（21）=-x（4）*x（2）+x（1）*（x（3）-x（5））;
%?dx（31）=x（3）-x（1）*x（2）;
%?dx（41）=0;
%?dx（51）=0;

????end
%?a=10;b=8/3;
%?
%?tspan1=[0?1];
%?tspan2=[0?50];
%?x0=[1?1?0];
z=[];
?for?r=linspace（15001000）??????????????????%linspace（x0xnn）?????
????[T?Y]=ode45（@lorenz[0?1][1;1;1;16;r;4]）;?????
????[T?Y]=ode45（@lorenz[0?50]Y（length（Y）:））;
?????Y（:1）=Y（:2）-Y（:1）;???????%去掉開始一段（暫態過程），保留剩余迭代點。
????%?對計算結果進行判斷，如果滿足x=y，則取點
?????for?k=2:length（Y）
????????f=k-1;
????????if?Y（k1）<0
????????????if?Y（f1）>0
????????????????y=Y（k2）-Y（k1）*（Y（f2）-Y（k2））/（Y（f1）-Y（k1））;
???????????????