syms f(x) c

f(x)=c*sqrt(x)

df(x)=diff(f(x))

%1-1
syms h
E=2*pi*int(f(x)*sqrt(1+df(x)^2),0,h)

syms c h positive
E=2*pi*int(f(x)*sqrt(1+df(x)^2),0,h)


%1-2
M=int(sqrt(1+df(x)^2),0,h)

simplify(M)
%2
findsym(M)
solve(M==10)

%example for c=1
subs(M,c,1)
solve(subs(M,c,1)==10)

call=[0.8:0.1:1.2];
for k=1:length(call)
    a(k,1)=call(k);
    a(k,2)=double(solve(subs(M,c,call(k))==10));
    a(k,3)=double(subs(E,[h,c],[a(k,2),call(k)]));
end
a

 for k=1:length(call)
    if a(k,2)>9.3
        a(k,2)=0;
        a(k,3)=0;
    end
end
a