subplot(1,2,1)
mesh(t,x,a1)
V = AXIS;
AXIS([0, 21, V(3:6) ] )
title('Angiostatin')
xlabel('Time')
ylabel('Position')
zlabel('Concentration of Angiostatin')

subplot(1,2,2)
mesh(t,x,c1)
V = AXIS;
AXIS([0, 21, V(3:6) ] )
title('Protease')
xlabel('Time')
ylabel('Position')
zlabel('Concentration of Protease')

figure

subplot(1,2,1)
mesh(t,x,eta1)
V = AXIS;
AXIS([0, 21, V(3:6) ] )
title('Endothelial Cells')
xlabel('Time')
ylabel('Position')
zlabel('Concentration of Endothelial Cells')

subplot(1,2,2)
mesh(t,x,f1)
V = AXIS;
AXIS([0, 21, V(3:4), 0, 1 ] )
title('Fibronectin')
xlabel('Time')
ylabel('Position')
zlabel('Concentration of Fibronectin')

figure

subplot(1,2,1)
mesh(t,x,pi1)
V = AXIS;
AXIS([0, 21, V(3:6) ] )
title('Inhibitor')
xlabel('Time')
ylabel('Position')
zlabel('Concentration of Inhibitor')

subplot(1,2,2)
mesh(t,x,rmac1)
V = AXIS;
AXIS([0, 21, V(3:6) ] )
title('Macrophages')
xlabel('Time')
ylabel('Position')
zlabel('Concentration of Macrophage Cells')

figure

subplot(1,2,1)
mesh(t,x,sig1)
V = AXIS;
AXIS([0, 21, V(3:6) ] )
title('Pericytes')
xlabel('Time')
ylabel('Position')
zlabel('Concentration of Pericyte cells')

subplot(1,2,2)
mesh(t,x,u1)
V = AXIS;
AXIS([0, 21, V(3:6) ] )
title('Chemotactic agent')
xlabel('Time')
ylabel('Position')
zlabel('Concentration of Chemotactic agent')

figure

subplot(1,2,1)
mesh(t,x,v1)
V = AXIS;
AXIS([0, 21, V(3:6) ] )
title('Angiogenic factor')
xlabel('Time')
ylabel('Position')
zlabel('Concentration of Angiogenic factor')
