Hello, I'm trying to solve vector polygon question. It's quite complicated,so i wrote it to matlab, and i'm getting error. Can anybody help me?

I'm trying to solve vector polygon question. It's quite complicated,so i wrote it to matlab, and i'm getting error.

Program:

"

x = sym('x');
fi1 = sym('fi1');
fi3 = sym('fi3');
i3 = sym('i3');
i1 = sym('33'); 
i2 = sym('30'); 
i7=sym('18');
y = sym('y');

fi3=2*pi+2*(atan((i3-sqrt(2)* sqrt(-162+i3^2))/(18+i3)))
%fi3=2*pi+2*(atan((sqrt(2)* sqrt(-162+i3^2))/(18+i3)))
fi1=2*pi+2*atan((((2178-y^2)^0.5)-33)/(x-33))
%fi1=2*pi+2*atan(((-(2178-y^2)^0.5)-33)/(x-33))

d=solve(i1^2+i3^2-i2^2+i1*i3*cos(fi1-fi3)+x*[i1*cos(fi1)+i3*cos(fi3)]+x^2,i3)

   d2=((20718*x^2*y^2 - 817638624*x - 324*x^3*y^2 - 3*x^2*y^4 + 6*x^4*y^2 + 883872*x^2*(2178 - y^2)^(1/2) - 15048*x^3*(2178 - y^2)^(1/2) + 264*x^4*(2178 - y^2)^(1/2) + 266112*y^2*(2178 - y^2)^(1/2) + x*(319519672582252032*x + 15908494561344*x^2*y^2 - 222521904000*x^3*y^2 - 4667000976*x^2*y^4 + 5889852144*x^4*y^2 + 41675040*x^3*y^4 - 41884128*x^5*y^2 + 463068*x^2*y^6 - 1044666*x^4*y^4 + 697500*x^6*y^2 - 1032*x^3*y^6 + 2340*x^5*y^4 - 2088*x^7*y^2 - 7*x^2*y^8 + 28*x^4*y^6 - 42*x^6*y^4 + 28*x^8*y^2 + 372213278908416*x^2*(2178 - y^2)^(1/2) - 6535809865728*x^3*(2178 - y^2)^(1/2) + 178286619456*x^4*(2178 - y^2)^(1/2) - 1809343008*x^5*(2178 - y^2)^(1/2) + 32997888*x^6*(2178 - y^2)^(1/2) - 220176*x^7*(2178 - y^2)^(1/2) + 2640*x^8*(2178 - y^2)^(1/2) - 7574107419648*y^2*(2178 - y^2)^(1/2) - 2877259968*y^4*(2178 - y^2)^(1/2) + 247104*y^6*(2178 - y^2)^(1/2) - 274110554234880*x*y^2 + 75823045056*x*y^4 - 7688736*x*y^6 + 126*x*y^8 - 6849698145725952*x*(2178 - y^2)^(1/2) - 17341459944956928*x^2 + 305128327477248*x^3 - 8410551732000*x^4 + 86586895296*x^5 - 1487713392*x^6 + 7897824*x^7 - 116559*x^8 + 654*x^9 - 7*x^10 + 788518683152640*y^2 + 106220572128*y^4 - 38733552*y^6 + 657*y^8 + 41241014899983360*(2178 - y^2)^(1/2) - 255558644352*x^2*y^2*(2178 - y^2)^(1/2) + 3232291392*x^3*y^2*(2178 - y^2)^(1/2) + 50903424*x^2*y^4*(2178 - y^2)^(1/2) - 84148416*x^4*y^2*(2178 - y^2)^(1/2) - 312048*x^3*y^4*(2178 - y^2)^(1/2) + 486288*x^5*y^2*(2178 - y^2)^(1/2) - 2640*x^2*y^6*(2178 - y^2)^(1/2) + 7920*x^4*y^4*(2178 - y^2)^(1/2) - 7920*x^6*y^2*(2178 - y^2)^(1/2) + 4317131934720*x*y^2*(2178 - y^2)^(1/2) - 815756832*x*y^4*(2178 - y^2)^(1/2) + 45936*x*y^6*(2178 - y^2)^(1/2) - 1928017446574222080)^(1/2) - 33*(319519672582252032*x + 15908494561344*x^2*y^2 - 222521904000*x^3*y^2 - 4667000976*x^2*y^4 + 5889852144*x^4*y^2 + 41675040*x^3*y^4 - 41884128*x^5*y^2 + 463068*x^2*y^6 - 1044666*x^4*y^4 + 697500*x^6*y^2 - 1032*x^3*y^6 + 2340*x^5*y^4 - 2088*x^7*y^2 - 7*x^2*y^8 + 28*x^4*y^6 - 42*x^6*y^4 + 28*x^8*y^2 + 372213278908416*x^2*(2178 - y^2)^(1/2) - 6535809865728*x^3*(2178 - y^2)^(1/2) + 178286619456*x^4*(2178 - y^2)^(1/2) - 1809343008*x^5*(2178 - y^2)^(1/2) + 32997888*x^6*(2178 - y^2)^(1/2) - 220176*x^7*(2178 - y^2)^(1/2) + 2640*x^8*(2178 - y^2)^(1/2) - 7574107419648*y^2*(2178 - y^2)^(1/2) - 2877259968*y^4*(2178 - y^2)^(1/2) + 247104*y^6*(2178 - y^2)^(1/2) - 274110554234880*x*y^2 + 75823045056*x*y^4 - 7688736*x*y^6 + 126*x*y^8 - 6849698145725952*x*(2178 - y^2)^(1/2) - 17341459944956928*x^2 + 305128327477248*x^3 - 8410551732000*x^4 + 86586895296*x^5 - 1487713392*x^6 + 7897824*x^7 - 116559*x^8 + 654*x^9 - 7*x^10 + 788518683152640*y^2 + 106220572128*y^4 - 38733552*y^6 + 657*y^8 + 41241014899983360*(2178 - y^2)^(1/2) - 255558644352*x^2*y^2*(2178 - y^2)^(1/2) + 3232291392*x^3*y^2*(2178 - y^2)^(1/2) + 50903424*x^2*y^4*(2178 - y^2)^(1/2) - 84148416*x^4*y^2*(2178 - y^2)^(1/2) - 312048*x^3*y^4*(2178 - y^2)^(1/2) + 486288*x^5*y^2*(2178 - y^2)^(1/2) - 2640*x^2*y^6*(2178 - y^2)^(1/2) + 7920*x^4*y^4*(2178 - y^2)^(1/2) - 7920*x^6*y^2*(2178 - y^2)^(1/2) + 4317131934720*x*y^2*(2178 - y^2)^(1/2) - 815756832*x*y^4*(2178 - y^2)^(1/2) + 45936*x*y^6*(2178 - y^2)^(1/2) - 1928017446574222080)^(1/2) + 422928*x*y^2 + 30*x*y^4 + 12388464*x*(2178 - y^2)^(1/2) - 33166584*x^2 + 727056*x^3 - 22239*x^4 + 294*x^5 - 3*x^6 - 23757624*y^2 + 1521*y^4 - 817638624*(2178 - y^2)^(1/2) - 264*x^2*y^2*(2178 - y^2)^(1/2) - 2376*x*y^2*(2178 - y^2)^(1/2) + 40473111888)/(4*(132*y^2*(2178 - y^2)^(1/2) - 2*x^2*y^2 - 132*x^2*(2178 - y^2)^(1/2) - 574992*x + 132*x*y^2 + 8712*x*(2178 - y^2)^(1/2) + 13068*x^2 - 132*x^3 + x^4 - 13068*y^2 + y^4 - 574992*(2178 - y^2)^(1/2) + 28462104)))^(1/2)

%result  y(x)

d3=solve(y^2+(x-(i2^2-y^2)^0.5)^2-i1^2-d2^2-i1*d2*cos(fi1-fi3),y)

%{d3 should be my y(x) function, but matlab shows error:

d3
Warning: Explicit solution could not be found. 

> In solve at 179

    In Untitled at 36 

d3 =

[ empty sym ]

}

[pfb]

can you format the code (there is a button for that, in the editor) and provide more information?

Answer 1

[Wojciech Krawczyk]

I modified it a little