global N
N=8
N_t_divisions = 8;
[T,Y] = ode45('example_7p4_func',[0,5], [0;0;0;0;0;0;0]);
stride = floor(size(T,1)/N_t_divisions);
Yplot = zeros(N_t_divisions,N-1);
Yplot(1,:) = Y(1,:);
Xsurf(1) = T(1);
for i=2:N_t_divisions
    Yplot(i,:) = Y(1+i*stride,:);
    Xsurf(i) = T(1+i*stride);
end
%Xsurf = T';
for i=1:N-1;Ysurf(i) = i/(N);end;
% Ysurf=Ysurf';
surf(Ysurf,Xsurf,Yplot)
size(T)
