

function [x,z] = LH82functionFF( L,W,  thetad, x0, z0, mu_range)

mu = [-1.33:0.01:1.33];
if (nargin<5)
    x = [];
    z = [];
    return;
end;

if (nargin==6),
    mu = mu_range;
end;


x1 = 3*mu.^2;
z1 = -mu.^3 + 2*mu;
  
x1 = x1*L;
z1 = z1*W;

G = [cosd(thetad) sind(thetad); -sind(thetad) cosd(thetad)];

TN = [x1; z1];
TT = G * TN;;
     
x = TT(1,:) + x0;     
z = TT(2,:) + z0;     

x = x';
z = z';
