function plot_nl(model_params, u_max, V_L, V_R),
if nargin < 1,
disp('USAGE: plotnl(model_params, u_max)');
end;
phi_max = asin(2*model_params.d0/model_params.bm);
phi_max = 9/10*phi_max;
if nargin < 2,
u_max = 200;
end;
if nargin < 3,
V_L = -500;
V_R = 500;
end;
y_vec = linspace(-phi_max, phi_max,100);
u_vec = linspace(-u_max, u_max, 3);
yout = zeros(length(u_vec), length(y_vec));
for h = 1:length(u_vec),
for k = 1:length(y_vec),
yout(h,k) = sfunmoment([],[],[u_vec(h), y_vec(k), V_L, V_R], 3, model_params);
end;
end;
figure; plot(y_vec, yout, 'b');
set(gca, 'FontSize', 16);
grid on; xlabel('Input -->', 'FontSize', 18);
ylabel('Output -->', 'FontSize', 18);
title('Function Plot', 'FontSize', 20);
%% End of File