clear
sim('gangguan_matahari.slx')
N = length(time);
for i = 1:N
    time_hour(i)= 6 + (time(i)/3600);
end
plot(time_hour,gangguan_matahari), grid on, xlabel('time[hour]'), ylabel('gangguan matahari [C]')