function [T,Y]=cancer(parameters,orders,TSim,Y0)
h=0.005;
n=round(TSim/h);
q1=orders(1);q2=orders(2);q3=orders(3);
b1=parameters(1);b2=parameters(2);b3=parameters(3);b4=parameters(4);b5=parameters(5);
b6=parameters(6);
cp1=1; cp2=1; cp3=1;
for j=1:n
    c1(j)=(1-(1+q1)/j)*cp1;
    c2(j)=(1-(1+q2)/j)*cp2;
    c3(j)=(1-(1+q3)/j)*cp3;
    cp1=c1(j); cp2=c2(j); cp3=c3(j);
end
x(1)=Y0(1); y(1)=Y0(2); z(1)=Y0(3);
for i=2:n
    x(i)=(1+b1*x(i-1)*(1-x(i-1))-x(i-1)*y(i-1))*h^q1-memo(x,c1,i);
    y(i)=(b2*y(i-1)*z(i-1)-b3*y(i-1))*h^q2-memo(y,c2,i);
    z(i)=(b4*z(i-1)*(1-z(i-1))-b5*y(i)*z(i-1)-b6*y(i))*h^q3-memo(z,c3,i);
end;
for j=1:n
    Y(j,1)=x(j);
    Y(j,2)=y(j);
    Y(j,3)=z(j);
end
T=h:h:TSim;