
i=0;
for T=1:0.1:400
    i=i+1;
    x1(i)=T;
    y1(i)=RollGeskeWhaley(50, 100, 50, 0.04, 20, 96, T);
    
    x2(i)=T;
    y2(i)=
end

plot(x,y);