chap2_3.m

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