global lamida 
        lamida=0.6627;
        nmesh=21
        nplot=nmesh;

        x=linspace(-1,1,nmesh)
        solinit=bvpinit(x,@guess)
        sol=bvp4c(@odes,@bcs,solinit)

        x=linspace(-1,1,nplot);
        y=deval(sol,x);
        
        y(1,:)
        plot(x,y(1,:));

