function [] = solutions_Final_p2c(t, theta1, L1, L2, h)

    [x, theta2] = solutions_Final_p2b(theta1, L2, L2, h);
    
    v = zeros(size(x));
    v(1:end-1) = (x(2:end) - x(1:end-1))./(t(2:end) - t(1:end-1));
    v(end) = (x(end) - x(end-1))/(t(end) - t(end-1));
    
    a = zeros(size(x));
    a(1:end-1) = (v(2:end) - v(1:end-1))./(t(2:end) - t(1:end-1));
    a(end) = (v(end) - v(end-1))/(t(end) - t(end-1));
    
    figure;
    subplot(2,3,1);
    plot(t, theta1);
    title('$\theta_1$ vs. $t$', 'interpreter', 'latex');
    xlabel('t');
    ylabel('$\theta_1$', 'interpreter', 'latex');
    
    subplot(2,3,2);
    plot(t, theta2);
    title('$\theta_2$ vs. $t$', 'interpreter', 'latex');
    xlabel('t');
    ylabel('$\theta_2$', 'interpreter', 'latex');
    
    subplot(2,3,3);
    plot(t, x);
    title('$x$ vs. $t$', 'interpreter', 'latex');
    xlabel('t');
    ylabel('$x$', 'interpreter', 'latex');
    
    subplot(2,3,4);
    plot(t, v);
    title('$v(t)$ vs. $t$', 'interpreter', 'latex');
    xlabel('t');
    ylabel('$v(t)$', 'interpreter', 'latex');
    
    subplot(2,3,5);
    plot(t, a);
    title('$a(t)$ vs. $t$', 'interpreter', 'latex');
    xlabel('t');
    ylabel('$a(t)$', 'interpreter', 'latex');

end

