how to calculate Z S from the equation

a(1)=0;

b(1)=1;

c(1)=0;

k3=2.8;

k1=0.33;

k2=1;

u=0;

v=0;

w=0;

S=0;

Z=0;

T=0.1;

t=0:0.01:0.1;% t from 0 to 1 interval 0.01

tau = t ./ T;

for i=1:20

for p=1:i-1

a(i+1)=(1/(i+1)).*(a(i)-b(i))*k2;

S=S+a(p)*c(i-p);

Z=Z+a(p)*b(i-p);

b(i+1)=(1/(i+1)).*(a(i).*k3-T*S-T*b(i));

c(i+1)=(1/(i+1))*(T*Z-k1*c(i));

u=u+(a(i).*(t/T).^i);

v=v+(b(i).*(t/T).^i);

w=w+(c(i).*(t/T).^i);

disp([u v w])

end

end

sir any one you are req to write the Z S calculation from equation of a(i),b(i),c(i) how to calculate the value

the phii=i+1;

I have puzzle about calculation of Z ,S

[Jan]

You have been asked repeatedly to use a proper code formatting. Please read an consider this: https://www.mathworks.com/matlabcentral/answers/help/rtc#rtc_summary . Thanks.