clc
clear all

%%
data = 100;                % MB
call = 0:500;            % minute

option_rate = 0.9;

% over data charge 50¿ø/1MB
% over call charge 108¿ø/1ºÐ

slim = zeros(length(call),1);
slim_op = zeros(length(call),1);
light = zeros(length(call),1);
light_op = zeros(length(call),1);
talk = zeros(length(call),1);
talk_op = zeros(length(call),1);
medium = zeros(length(call),1);
medium_op = zeros(length(call),1);

for c=1:length(call)
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    over_d = data - 100;
    if over_d < 0
        over_d = 0;
    end
    over_c = call(c) - 150;
    if over_c < 0
        over_c = 0;
    end

    slim(c,1) = (35000 + 108*over_c + 50*over_d)*1.1 + (22000-5500);
    slim_op(c,1) = (35000 + 108/2*over_c*option_rate + 108*over_c*(1-option_rate) + 50*over_d + 2500)*1.1 + (22000-5500);

    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    over_d = data - 500;
    if over_d < 0
        over_d = 0;
    end
    over_c = call(c) - 200;
    if over_c < 0
        over_c = 0;
    end

    light(c,1) = (45000 + 108*over_c + 50*over_d)*1.1 + (22000-8800);
    light_op(c,1) = (45000 + 108/2*over_c*option_rate + 108*over_c*(1-option_rate) + 50*over_d + 2500)*1.1 + (22000-8800);

    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    over_d = data - 100;
    if over_d < 0
        over_d = 0;
    end
    over_c = call(c) - 250;
    if over_c < 0
        over_c = 0;
    end

    talk(c,1) = (45000 + 108*over_c + 50*over_d)*1.1 + (22000-8800);
    talk_op(c,1) = (45000 + 108/2*over_c*option_rate + 108*over_c*(1-option_rate) + 50*over_d + 2500)*1.1 + (22000-8800);
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    over_c = call(c) - 300;
    if over_c < 0
        over_c = 0;
    end

    medium(c,1) = (55000 + 108*over_c)*1.1 + (22000-14300);
    medium_op(c,1) = (55000 + 108/2*over_c*option_rate + 108*over_c*(1-option_rate))*1.1 + (22000-14300);

end

plot(call, slim, call, slim_op, call, light, call, light_op, call,talk, call, talk_op, call, medium, call, medium_op, 'LineWidth', 1.5);
legend1 = legend('Slim', 'Slim op', 'Light', 'Light op', 'Talk', 'Talk op', 'Medium', 'Medium op', 8);
set(legend1,'Location','NorthWest');

% slim = 35000 + (22000-5500);
% light = 45000 + (22000-8800);
% talk = 45000 + (22000-11000);
% medium = 55000 + (22000-14300);