% EX2_ISTARAL
%   Be sure to append the file fxy_test.m to the routine
%   Istaral.m
%
global N alpha beta
N=50;
al=-.9:.02:.9; be=-.9:.02:.9;
for i=1:size(al,2)
  for j=1:size(be,2)
    alpha=al(i); beta=al(j);
    if beta<=alpha
      Z(i,j)=acc_Istaral(alpha,beta);
    end
  end
end
meshc(al,be,Z)
axis([-1 1 -1 1 1 16.5])
xlabel('\alpha','FontSize',14)
ylabel('\beta','FontSize',14)
zlabel('acc','Fontsize',14)
