clear?all;
close?all;
global?C?M0?F

ts=0.01;
T=30;
TimeSet=[0:ts:T];

para=[];
options=odeset（‘RelTol‘1e-3‘AbsTol‘[1e-3?1e-3?1e-3?1e-3]）;
%options=[];
x0=[0.50.300];
[txout]=ode45（‘chap2_3eq‘TimeSetx0optionspara）;
x1=xout（:1）;
x2=xout（:2）;
x3=xout（:3）;
x4=xout（:4）;

s=C（1）*x1+C（2）*x2+C（3）*x3+C（4）*x4;

if?F==1
???M0=40;
???u=-M0*sign（s）;
elseif?F==2
???beta=3