function parciaal

disp('Ejercicio A=1')
disp('Ejercicio B=2')
disp('Para terminar el programa ingrese 0')
n=input('ingrese el número del ejercicio que se desee realizar:')
switch n
case 0
disp('el programa ha finalizado')
case 1
syms x
k1=40000
k2=40
m=95
g=9.8
h=0.43
f=2*k2*x^(1/2)/5+(1/2)*k1*x^2-2*m*g*h
ezplot(f),grid
a=input('ingrese el valor de a=')
b=input('ingrese el valor de b=')
fa=double(subs(f,a));
fb=double(subs(f,b));
x=(a+b)/2;
T=1;
c=1;
disp('++++++++++++++++++++++++++++++++++')
disp('iteracion a b x fa fx')
disp('++++++++++++++++++++++++++++++++++')
disp([ c a b x])
if fa*fb < 0
while T > 0.001
fx=double(subs(f,x));
if fa*fx < 0
b=x;
else
a=x
end
x1=(a+b)/2;
T=abs((x1-x)/x1);
x=x1;
c=c+1;
disp([c a b x fa fx T])
end
fprintf('El valor de la raiz es %1.3f \n',x)
else
disp('no tiene raices')
end
case 2

end
end