Построить неявную функцию, используя итеративный подход

Question

Это код, который я использую для построения графика между 100 L и h.

syms x y L h;
p(1) = (L.*h.*x^2)/2;
for i = 2:25
    integrand = @(x,y) ((h*L).*p(i-1));
    u(i) = int(int(integrand(x,y),0,y),0,x);
    p(i) = (1+h).*p(i-1) + u(i);
end
s = sum(p(i));
s = vpa(s, 5);
q = subs(s,x,1);
% q = simplify(q, 'Steps',70);
q = vpa(q, 5);   

Из этого я получаю значения q, и оттуда я пытаюсь построить L и h, но это не фактическая кривая.