close all
clear all

g = @(u,v) u.^2-v.^2*sin(u).^2;

n_max = 1000;
u_max = pi;

v = 2;

for n = 1:n_max
    u = (n-1)/(n_max-1)*u_max;
    x(n) = u;
    lhs)n) = u.^2;
    rhs(n) = g(u,v);
end

plot(x,lhs,'b',x,rhs,'r','LineWidth',2)