tcolc.m

syms x ;
syms a1;

cc1=0.25;
u=(x-x^2);
u1=a1*u;
x1=0.5;
L(x)=diff(u,2)-4*u-4*x;
E(x)=-2-4*(x-x^2)*a1-4*x;