How can I solve this system of equations?
2 views (last 30 days)
Show older comments
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!
0 Comments
Answers (1)
See Also
Categories
Find more on Equation Solving in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!