function compPlot() 

global T      %time
global U      %parameter
global DYDP   %gradient

U = parameterfile;
y0 = initvaluefile;

a1 = 3;                                 % creates a figure
a2 = 4;
a3 = a1*a2;
figure();

    for ii = 1:length(y0)               % iteration over the component
        
        pos = mod(ii,a3);
            if(pos == 0)
                pos = a3;
            end
        subplot(a1,a2,pos);        
        result = zeros(length(T),length(U));
        
        for jj = 1:length(T)            % iteration over the time
            for k = 1:length(U)         % iteration over the parameter
            result(jj,k) = abs(DYDP(jj,ii,k));
            end
        end
        
        mesh(result)                    % plot of the result
        xlabel('Parameter')
        ylabel('Time [s]')
        zlabel('Gradient')
        title(['Component ' num2str(ii)])
        
        if(pos == a3) && (ii < length(y0))
        figure();  
        end
    end
end