phi0[1]=7^3*(L+e3)^3/(3^2 *2 * (T/2)^7);
phi0[2]=(T-4)/2 * phi0[1];
phi[2]=phi0[2]-phi0[1]*t;
phi[1]=phi0[1]*ones(points);
m = (L+e3)/(6/7 * (T/2)^(7/3));
c = (T - 4)/2;
x[2]= 3 * m/4 * ((2 + c)^(4/3)*ones(points) - ((2+c)*ones(points)-t).^(4/3));
x[1]= 2*x[2] + 3*m/4 * (2+c)^(4/3)*t + 9*m/28 * (((2+c)*ones(points)-t).^(7/3)- (2 + c)^(7/3)*ones(points));
j = m^4 * (T/2)^(7/3)*3/7 * (3^(7/3) - 1);
потому что я не умею писать красивый код. да и сильно это все не улучшить.