syms y
h=[1.5902    1.7614    1.9259    2.085    2.2393    2.3897    2.5366    2.6804    2.8216    2.9603    3.0969    3.2315    3.3644    3.4957    3.6256    3.7542    3.8816    4.0079    4.1332    4.2576    4.3812    4.504    4.6262    4.7478    4.8689    4.9895    5.1097    5.2295    5.349    5.4682    5.5872    5.7061    5.8248    5.9435    6.0622    6.1809    6.2996    6.4185    6.5375    6.6567    6.7763    6.7854    6.9053    6.9082    7.0285    7.1491    7.2703    7.3919    7.5142    7.637    7.6416    7.7653    7.8899    8.0154    8.1419    8.2695    8.3983    8.5284    8.66    8.7932    8.9282    8.9284    9.0653    9.2045    9.3461    9.4905    9.638    9.7891    9.9443    10.1043    10.2699    10.4425    10.6237    10.8159    11.0233    11.2532    11.5236    11.9349
]
A=double(24.5*17.8*int(sqrt(1-y^2./60^2),y,-60,h-60))
A

解答之一是：
clear;clc

syms y h

a=24.5*17.8*int(sqrt(1-y^2/60^2),y,-60,h-60);

H=[...];

A=subs(a,H);

plot(H,A,'r*')