求解非线性方程周期解matlab程序！

function shootingmethod
m=20;h=0.3142;t0=0.0;ep=0.00000001;
x（1）=0.06594;x（2）=0.01;a=1;b=1;
c（1）=0.5*h;c（2）=c（1）;c（5）=c（1）;
c（3）=h;c（4）=h;s（1）=x（1）;s（2）=x（2）;
x（1）=1.0;x1（2）=0.0;
x2（1）=0.0;x2（2）=1.0;
v=1;
while v==1
    t=t0;
    x（1）=s（1）;x（2）=s（2）;
    for i=1:m
        t1=t
        ts=t
        for ii=1:2
            p（ii）=x（ii）;
            w（ii）=x（ii）;
        end
        for jj=1:4
            f（1）=p（2）;
            f（2）=-0.2*p（2）-4.0*p（1）-p（1）.^3+0.03*cos（t）;
            t=ts+c（jj）;
            for ii=1:2
                p（ii）=c（jj）.*f（ii）+x（ii）;
                x（ii）=c（jj+1）.*f（ii）/3.0+x（ii）;
            end
        end
        t=t1;ts=t;
        for ii=1:2
            p（ii）=x1（ii）;
            w（ii）=x1（ii）;
        end
        for jj=1:4
            f（1）=p（2）;
            f（2）=-0.2*p（2）-（4.0+3.0*x（1）.^2）.*p（1）;
            t=ts+c（jj）;
            for ii=1:2
                p（ii）=c（jj）.*f（ii）+w（ii）;
                x1（ii）=c（jj