% -------------------------------------------------------------------------------
%  Assigment I, Problem 4
% -------------------------------------------------------------------------------
% 
% Calculating the probability of t-student distribution
% Juan F Rubio-Ramirez
% Atlanta, 06-12-2005


clear all;
close all;
clc
% number of simulations
N=1000000;
mean=0;
var=1;
x(1)=0;
x(2)=0.67;
x(3)=1.28;
x(4)=1.65;
x(5)=2.32;
x(6)=2.58;
x(7)=3.09;
x(8)=3.72;
for d=1:10:100;
    R=random('t',d,1,N);
    test1=zeros(1,8);
    for i=1:8,
        for j=1:N,
            if (R(j)<x(i));
                test1(i)=test1(i)+1;
            end
        end
        test1(i)=test1(i)/N;
    end
    R=random('norm',mean,var,1,N);
    test=zeros(1,8);
    for i=1:8,
        for j=1:N,
            if (R(j)<x(i));
                test(i)=test(i)+1;
            end
        end
        test(i)=test(i)/N;
    end
    test-test1
end