Exame Modelo de AIPE

sim, já percebi onde querias chegar...

Obrigado pela ajuda.

em média a solução mais rápida seria:

syms H R; %Retirei aqui a declaração de "real"
V=2/3*pi*R^3+pi*R^2*H;
H=solve(V-500,H);
C=2*pi*H*R*300+2*pi*R^2*400;
f=diff(C);
solucao=solve(f);

for i=1:3
if abs(solucao(i))==solucao(i) %Outra maneira de verificar se é a solução pretendida
%Também seria possível if imag(solucao(i))==0
resultado=solucao(i);
break
else
continue
end
end

RAIO=double(resultado)

% Não é grande ajuda, mas é uma maneira de optimizar...

Administrador escreveu:Vejam lá se descobrem o que isto faz:

syms H R real
V=2/3*pi*R^3+pi*R^2*H
H=solve(V-500,H)
C=2*pi*H*R*300+2*pi*R^2*400
f=diff(C)
solucao=solve(f)
for i=1:3
if conj(solucao(i))==solucao(i)
resultado=solucao(i);
break
else
continue
end
end
RAIO=double(resultado)

optimizar ainda mais seria retirar o "else" e o "continue" que eu coloquei dentro do ciclo for .