x = linspace(0,3*pi,1000000);
y1 = sin(x+0.5*pi)
y2 = sin(x)
y3 = -sin(x)
tiledlayout(2,1)
ax1 = nexttile;
plot(ax1,x,y1,x,y2,'LineWidth',2)
xlim([0 3*pi])
ylim([-1.1 1.1])
title(ax1,'Capacitor V - I phase relationship', 'FontSize', 14)
xlabel(ax1,'Time')
legend('Current','Voltage')
ax2 = nexttile;
plot(ax2,x,y1,x,y3,'LineWidth',2)
xlim([0 3*pi])
ylim([-1.1 1.1])
title(ax2,'Inductor V - I phase relationship', 'FontSize', 14)
xlabel(ax2,'Time')
legend('Current','Voltage')