function [X,f,g,h] = ecu1(N,L,H1)
%function [X,f] = ecu1(L,H1)

format long;
i=1;k=1;
for x=0:0.01:20
    f(i) = x*sin(x*L);
    g(i) = H1*cos(x*L);
    h(i) = x*sin(x*L)/cos(x*L) - H1;
    i=i+1;
end

X=0:0.01:20;

subplot(1,2,1),plot(X,f,X,g)
subplot(1,2,2),plot(X,h)

