cur_time_pulse = [ 0 , 1 , 7 , 8 , 10 ] + 1
cur_current_pulse = [ 0 , 1 , 1 , 0 , 0 ]
cur_x = deepcopy(cur_time_pulse)
cur_y = deepcopy(cur_current_pulse)
append!(cur_x, cur_time_pulse + 1*( maximum(cur_time_pulse) - 1 ) )
append!(cur_y, cur_current_pulse)
append!(cur_x, cur_time_pulse + 2*( maximum(cur_time_pulse) - 1 ) )
append!(cur_y, cur_current_pulse)
cur_plot = plot()
plot!([0,2,NaN,8,12,NaN,18,22,NaN,28,30], [1,1,NaN,1,1,NaN,1,1,NaN,1,1], color=2, label="Steady-State Current", style=:dash)
plot!(cur_x, cur_y, fill=true, fillopacity=0.2, label="Pulsed Current", color=1)
ylims!(0, 1.5)
xlims!(0, 30)
plot!(xticks=([10,20,30],[]),yticks=([0,0.5,1.0, 1.5],["0","",L"I_P",""]))
xlabel!("Time")
title!("Pulsed vs Steady-State Operation")