close all, clear all
%Form two vectors of the stats in (HP, speed, Atk, SpAtk, Def, SpDef)
P1=[1 1 1 1 1 0];
P2=[1 1 1 1 1 0];
%Items (have it equal to the stat you want to keep, 7 for Everstone, 8 for
%knot)
item1 = 7;
item2 = 8;


% P1 = zeros(1,6);
% P2 = zeros(1,6);
% P1(P1t)=1;
% P2(P2t)=1;
P1=P1*31;
P2=P2*31;

size=100000;
PIVs = [P1;P2];
IVs = zeros(size,6);
items = [item1 item2];
if sum(items == 8)>0
    parenttraits = 5;
else
    parenttraits = 3;
end

for i = 1:size
    traitskept = randperm(6);
    loopl=6;
    parenttraitsl = parenttraits;
    if sum(items<=6)>0
        IVs(i,items(items<=6))= PIVs(items<=6,items(items<=6));
        parenttraitsl = parenttraits-1;
        traitskept(traitskept==items(items<=6))=[];
        loopl=5;
    end
    for j = 1:parenttraitsl
        IVs(i,traitskept(j))=PIVs(randi(2),traitskept(j));
    end
    for j=(parenttraitsl+1):loopl
        IVs(i,traitskept(j))=randi(32)-1;
    end
end
PerfectIVs = (IVs==31);
PerfectIVs = sum(PerfectIVs,2);
for i=1:6
    perfectIVprob(i) = sum(PerfectIVs>=i)/size;
end
round(perfectIVprob*100)/100