使用迭代方法绘制隐式函数

Question

这是我用来在‍ 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，但这并不是实际的曲线。