clc;close all;clear all;

tic;

%n=9;

syms u1 u2 u3 u4 u5 u6 u7 u8 u9 u10  a0

u(1,1)=u1;

u(1,2)=u2;

u(1,3)=u3;

u(1,4)=u4;

u(1,5)=u5;

u(1,6)=u6;

u(1,7)=u7;

u(1,8)=u8;

u(1,9)=u9;

u(1,10)=u10;

t1=ones(10,10);

t2=eye(10);

y=t1-t2;

%x=zeros(1,10);

z0=100;

for i=1:11

    z(1,i)=a0;

end

%z(1,1)=a0

z(1,1)=subs(z(1,1),a0,z0);

%z=zeros(1,10)

for m=1:10

    for i=1:10

        if i<=m

            x(1,i)=u(1,i)*z(1,i);

            %xx(i)=u(1,i)*z(1,i)

        else

            x(1,i)=z(1,i);

        end

        if m==10 && i==m

                x(1,i)=0;

        end

        z(1,i+1)=z(1,i)-x(1,i)+2*x(1,i)*y(m,i);    

    end 

    w(1,m)=z(1,i+1);

    %w(1,m)

end

for m=1:10

    for i=1:10

        if i>m

            w(1,m)=subs(w(1,m),u(1,i),1);

        end

    end

    %w(1,m)

end

            

A=solve(w(1,1)-w(1,2),w(1,1)-w(1,3),w(1,1)-w(1,4),w(1,1)-w(1,5),w(1,1)-w(1,6),w(1,1)-w(1,7),w(1,1)-w(1,8),w(1,1)-w(1,9),w(1,1)-w(1,10));

u(1,1)=A.u1;

u(1,2)=A.u2;

u(1,3)=A.u3;

u(1,4)=A.u4;

u(1,5)=A.u5;

u(1,6)=A.u6;

u(1,7)=A.u7;

u(1,8)=A.u8;

u(1,9)=A.u9;

u(1,10)=subs(u(1,10),0/2);

r=u

sum=subs(w(1,1),u1,u(1,1))

toc;