## ekmantan.m -- compute and plot the tangent of the ekman angle

function ekmantan
  ratio = [0:0.1:3];
  tangle = ((1 + 4*ratio.^4).^0.5 - 1)/1.414;
  angle = 57.3*atan(tangle);
  plot(ratio,angle);
endfunction
