How can I solve this system of equations?

I'm trying to solve a system of equations of two equations with two unknowns but Matlab has been running for an hour and it is still showing busy and I hence I have no results. Is Matlab able to solve systems like this?

syms c_d c_e V_d;

ec1='0.9259259259e-2*sqrt(1/2+(1/2)*sqrt(1+2328.331561*c_d^2))+(sqrt(1/2+(1/2)*sqrt(1+c_e^2))*sinh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2)))-sqrt(-1/2+(1/2)*sqrt(1+c_e^2))*sin(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2))))/(cosh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2)))+cos(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2))))+(.2000000000*(2.521200000*10^10*A*sech(A)^2+2.489400000*10^10*tanh(A)))*(-0.9259259259e-2*sqrt(-1/2+(1/2)*sqrt(1+2328.331561*c_d^2))+(-sqrt(1/2+(1/2)*sqrt(1+c_e^2))*sin(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2)))-sqrt(-1/2+(1/2)*sqrt(1+c_e^2))*sinh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2))))/(cosh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2)))+cos(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_d^2)))))/(3.750000000*10^10*A*sech(A)^2+3.750000000*10^10*tanh(A)+1.590000000*10^8*tanh(A)^2)=1.528755555*10^7*V_d*((2.500000000*(A*sech(A)^2+tanh(A)))/sqrt(1/2+(1/2)*sqrt(1+2328.331561*c_d^2))-tanh(A)*(-sqrt(1/2+(1/2)*sqrt(1+c_e^2))*sin(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2)))-sqrt(-1/2+(1/2)*sqrt(1+c_e^2))*sinh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2))))/(sqrt(1/2+(1/2)*sqrt(1+c_e^2))*sqrt(-1/2+(1/2)*sqrt(1+c_e^2))*(cosh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2)))+cos(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2))))))/(3.750000000*10^10*A*sech(A)^2+3.750000000*10^10*tanh(A)+1.590000000*10^8*tanh(A)^2)';

ec2='-0.9259259259e-2*sqrt(-1/2+(1/2)*sqrt(1+2328.331561*c_d^2))+(-sqrt(1/2+(1/2)*sqrt(1+c_e^2))*sin(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2)))-sqrt(-1/2+(1/2)*sqrt(1+c_e^2))*sinh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2))))/(cosh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2)))+cos(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2))))-(.2000000000*(2.521200000*10^10*A*sech(A)^2+2.489400000*10^10*tanh(A)))*(0.9259259259e-2*sqrt(1/2+(1/2)*sqrt(1+2328.331561*c_d^2))+(sqrt(1/2+(1/2)*sqrt(1+c_e^2))*sinh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2)))-sqrt(-1/2+(1/2)*sqrt(1+c_e^2))*sin(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2))))/(cosh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2)))+cos(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2)))))/(3.750000000*10^10*A*sech(A)^2+3.750000000*10^10*tanh(A)+1.590000000*10^8*tanh(A)^2) = -1.528755555*10^7*V_d*(-2.500000000*sqrt(-1/2+(1/2)*sqrt(1+2328.331561*c_d^2))*(A*sech(A)^2+tanh(A))/(sqrt(1/2+(1/2)*sqrt(1+2328.331561*c_d^2))*(sqrt(1/2+(1/2)*sqrt(1+2328.331561*c_d^2))+1))-tanh(A)*(sqrt(1/2+(1/2)*sqrt(1+c_e^2))*sinh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2)))-sqrt(-1/2+(1/2)*sqrt(1+c_e^2))*sin(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2)))-tanh(A)*(cosh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2)))+cos(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2)))))/(sqrt(1/2+(1/2)*sqrt(1+c_e^2))*sqrt(-1/2+(1/2)*sqrt(1+c_e^2))*(cosh(2*A*sqrt(1/2+(1/2)*sqrt(1+c_e^2)))-cos(2*A*sqrt(-1/2+(1/2)*sqrt(1+c_e^2))))))/(3.750000000*10^10*A*sech(A)^2+3.750000000*10^10*tanh(A)+1.590000000*10^8*tanh(A)^2)';

[c1 c2]=solve(ec1,ec2,c_d,c_e)

Thank you!

Answer 1

[Torsten]

Maybe numerically, but not symbolically.

Try fsolve.

Best wishes

Torsten.

[Lorena Ortiz]

Thank you Torsten. I'll try.