format short
clear all    %just some matlab stuff
close all

myhaste=856;            %you personal hasterating, change this number to get your table
AB=2.5;                 %arcane blast initial cast time
n=5;                    %the number of situations we're looking at
ABcast=zeros(5,n,700);  %matrix that will contain information

% increase in castspeed from different buffs

totem=1.05;             %spellhaste totem
moonkin=1.03;           %the 3% raidhaste
pot=500;                %potion of speed
set2=1.12;              %the t10 2set bonus
BM=250;                 %black magic enchant
BL=1.3;                 %bloodlust
IV=1.2;                 %icy veins
talents=1.06;           %talents

% probabilty table for MBarr procs
mb=zeros(6,1);
mb(1)=0.8704;
mb(2)=0.0519;
mb(3)=0.0311;
mb(4)=0.0186;
mb(5)=0.0112;


%stating what range of hasteratings we want to look at
hastemin=400;   
hastemax=1200;

min1=hastemin-1;

%some empty vectors used for plotting later
cast=zeros(701,1);
castBL=zeros(701,1);
castBLIV=zeros(701,1);
castIVpot=zeros(701,1);
one=ones(hastemax-hastemin+1,1);

%index different situations
factor=[1;BL;(BL*IV);(BL);(IV)];


for HR=hastemin:hastemax 

   
for k=1:n    %loop for looking at different situations                   
   
    HRl=HR;
    if k==4                                 %adding 500 hasterating from pot
        HRl=HR+500;                         %if we are looking at situation
    elseif k==5                             %4 or 5.
        HRl=HR+500;
    end

    gear=((HRl/3279)+1);                                %haste rating haste gains
    basef=1/(talents*gear*totem*moonkin*factor(k));     %total casttime reduction without cd's
    ABbasic=AB*basef;                                   %casttime for AB
    
    basef=1/(talents*gear*totem*moonkin*factor(k)*set2);
    AB2set=AB*basef;                                        %casttime for AB with 2setbonus
    
    gear=(((HRl+250)/3279)+1);                              %haste rating haste gains with black magic
    basef=1/(talents*gear*totem*moonkin*factor(k)*set2);    %total casttime reduction without cd's
    AB2setBM=AB*basef;                                      %casttime for AB with 2setbonus and Black Magic proc

    ABcast(1,k,HR-min1)=ABbasic;             %saving the above values in the matrix
    ABcast(2,k,HR-min1)=AB2set;
    ABcast(3,k,HR-min1)=AB2setBM;

    casts=ceil(5/AB2set)-1;         %the number of AB's you can cast having the buff from 2set bonus
    castsBM=ceil(5/AB2setBM)-1;     %the number of AB's you can cast having the buff from 2set bonus and black magic
        if casts<castsBM            
        change=1;                   %lets us know if there is a difference in number of casts having the 2setbuff   
        end                         %from having both2set and black magic

    ABcast(4,k,HR-min1)=casts;

%creating a vector containing the uptime of buff for the 
%different rotations in the vector 'mb' earlier.
upt=zeros(6,1);
   for l=1:5
       upt(l)=((1+casts)*100)/(l+4);
        if upt(l)>100
        upt(l)=100;
        end
    end

totupt2set=sum(upt.*mb); %multiplying the values of mb and upt vectors
if totupt2set>100        %making sure we don't get above 100% uptime
   totupt2set=100;
end
ABcast(5,k,HR-min1)=totupt2set; %adding uptime information in the matrix ABcast

end
%saving the information from the matrix into vectors for the different

cast(HR-399)=ABcast(3,1,HR-min1);
castBL(HR-399)=ABcast(3,2,HR-min1);
castBLIV(HR-399)=ABcast(3,3,HR-min1);
castIVpot(HR-399)=ABcast(3,5,HR-min1);

castnoBM(HR-399)=ABcast(2,1,HR-min1);
castnoBMBL(HR-399)=ABcast(2,2,HR-min1);
castnoBMBLIV(HR-399)=ABcast(2,3,HR-min1);
castnoBMIVpot(HR-399)=ABcast(2,5,HR-min1);


end


disp('                      no cooldowns   BL   BL+IV   BL+pot    IV+pot')
disp('AB cast')
disp('AB cast,2set')
disp('AB cast,2set,BM')
disp('ABs cast with 2set')
disp('2set uptime')

ABcast(:,:,myhaste-min1)     %shows the table for your haste rating
hasterating=(hastemin:1:hastemax); %making a vector for the values of hasteerating, spaced 1 apart from 400 to 1100

hold on
plot(hasterating,cast)
plot(hasterating,castBL)
plot(hasterating,castBLIV)
plot(hasterating,castIVpot)
plot(hasterating,one,'-r')
axis([hastemin hastemax 0.8 1.8])
xlabel('hasterating')
title('AB and Mbarr cast times with 2setbonus and Black Magic')
        legend('No cooldowns','Bloodlust','Icy veins+speed potion,','Bloodlust+Icy veins');
        
figure
hold on
plot(hasterating,castnoBM)
plot(hasterating,castnoBMBL)
plot(hasterating,castnoBMBLIV)
plot(hasterating,castnoBMIVpot)
plot(hasterating,one,'-r')
axis([hastemin hastemax 0.8 1.8])
xlabel('hasterating')
title('AB and Mbarr cast times with 2setbonus')
        legend('No cooldowns','Bloodlust','Icy veins+speed potion,','Bloodlust+Icy veins');