Sponsored links: Algebra eBooks
 

Related

linel-linsolve-string

linel: 500;

eq1:a*x_1^5 + b*x_1^4...

eq2:a*x_2^5 + b*x_2^4...

Calculate

linel-linsolve-string

linel: 200;

eq1:a*x_1^5 + b*x_1^4...

eq2:a*x_2^5 + b*x_2^4...

Calculate

linel-linsolve-string

linel: 500;

eq1:a*x_1^5 + b*x_1^4...

eq2:a*x_2^5 + b*x_2^4...

Calculate

linel-linsolve-string

linel: 200;

eq1:a*x_1^5 + b*x_1^4...

eq2:a*x_2^5 + b*x_2^4...

Calculate

linel

Run Example
(%i1)linel: 200;
(%o1)                                                                                             200
(%i2) eq1:a*x_1^5 + b*x_1^4+ c*x_1^3+ d*x_1^2  + 1 = y_1;
                                                                                   5        4        3        2
(%o2)                                                                         a x_1  + b x_1  + c x_1  + d x_1  + 1 = y_1
(%i3) eq2:a*x_2^5 + b*x_2^4+ c*x_2^3+ d*x_2^2 + 1 = y_2;
                                                                                   5        4        3        2
(%o3)                                                                         a x_2  + b x_2  + c x_2  + d x_2  + 1 = y_2
(%i4) eq3:5*a*x_2^4 + 4*b*x_2^3+ 3*c*x_2^2+ 2*d*x_2  = 0;
                                                                                     4          3          2
(%o4)                                                                         5 a x_2  + 4 b x_2  + 3 c x_2  + 2 d x_2 = 0
(%i5) eq4:a*x_3^5 + b*x_3^4+ c*x_3^3+ d*x_3^2 + 1 = y_3;
                                                                                   5        4        3        2
(%o5)                                                                         a x_3  + b x_3  + c x_3  + d x_3  + 1 = y_3
(%i6) string(linsolve([eq1,eq2,eq3,eq4],[a,b,c,d]));
(%o6) [a = (x_1^3*(x_2^4*(2*y_3-2)+(2*x_3^4-4*x_2^2*x_3^2)*y_2-2*x_3^4+4*x_2^2*x_3^2)+x_1^2*(x_2^5*(1-y_3)+(4*x_2^2*x_3^3-3*x_2*x_3^4)*y_2+3*x_2*x_3^4-4*x_2^2*x_3^3)+x_1^4*(x_2^3*(1-y_3)+(3*x_2*x_3^2\
-2*x_3^3)*y_2+2*x_3^3-3*x_2*x_3^2)+(x_2^3*x_3^4-2*x_2^4*x_3^3+x_2^5*x_3^2)*y_1-x_2^3*x_3^4+2*x_2^4*x_3^3-x_2^5*x_3^2)/(x_1^2*(-x_2^5*x_3^5+2*x_2^6*x_3^4-x_2^7*x_3^3)+x_1^3*(2*x_2^4*x_3^5-3*x_2^5*x_3^\
4+x_2^7*x_3^2)+x_1^4*(-x_2^3*x_3^5+3*x_2^5*x_3^3-2*x_2^6*x_3^2)+x_1^5*(x_2^3*x_3^4-2*x_2^4*x_3^3+x_2^5*x_3^2)),b = -(x_1^3*(x_2^5*(3*y_3-3)+(2*x_3^5-5*x_2^3*x_3^2)*y_2-2*x_3^5+5*x_2^3*x_3^2)+x_1^5*(x\
_2^3*(1-y_3)+(3*x_2*x_3^2-2*x_3^3)*y_2+2*x_3^3-3*x_2*x_3^2)+x_1^2*(x_2^6*(2-2*y_3)+(5*x_2^3*x_3^3-3*x_2*x_3^5)*y_2+3*x_2*x_3^5-5*x_2^3*x_3^3)+(x_2^3*x_3^5-3*x_2^5*x_3^3+2*x_2^6*x_3^2)*y_1-x_2^3*x_3^5\
+3*x_2^5*x_3^3-2*x_2^6*x_3^2)/(x_1^2*(-x_2^5*x_3^5+2*x_2^6*x_3^4-x_2^7*x_3^3)+x_1^3*(2*x_2^4*x_3^5-3*x_2^5*x_3^4+x_2^7*x_3^2)+x_1^4*(-x_2^3*x_3^5+3*x_2^5*x_3^3-2*x_2^6*x_3^2)+x_1^5*(x_2^3*x_3^4-2*x_2\
^4*x_3^3+x_2^5*x_3^2)),c = (x_1^4*(x_2^5*(3*y_3-3)+(2*x_3^5-5*x_2^3*x_3^2)*y_2-2*x_3^5+5*x_2^3*x_3^2)+x_1^2*(x_2^7*(1-y_3)+(5*x_2^3*x_3^4-4*x_2^2*x_3^5)*y_2+4*x_2^2*x_3^5-5*x_2^3*x_3^4)+x_1^5*(x_2^4*\
(2-2*y_3)+(4*x_2^2*x_3^2-2*x_3^4)*y_2+2*x_3^4-4*x_2^2*x_3^2)+(2*x_2^4*x_3^5-3*x_2^5*x_3^4+x_2^7*x_3^2)*y_1-2*x_2^4*x_3^5+3*x_2^5*x_3^4-x_2^7*x_3^2)/(x_1^2*(-x_2^5*x_3^5+2*x_2^6*x_3^4-x_2^7*x_3^3)+x_1\
^3*(2*x_2^4*x_3^5-3*x_2^5*x_3^4+x_2^7*x_3^2)+x_1^4*(-x_2^3*x_3^5+3*x_2^5*x_3^3-2*x_2^6*x_3^2)+x_1^5*(x_2^3*x_3^4-2*x_2^4*x_3^3+x_2^5*x_3^2)),d = -(x_1^4*(x_2^5*(2*y_3-2)+(3*x_3^5-5*x_2^2*x_3^3)*y_2-3\
*x_3^5+5*x_2^2*x_3^3)+x_1^3*(x_2^6*(1-y_3)+(5*x_2^2*x_3^4-4*x_2*x_3^5)*y_2+4*x_2*x_3^5-5*x_2^2*x_3^4)+x_1^5*(x_2^4*(1-y_3)+(4*x_2*x_3^3-3*x_3^4)*y_2+3*x_3^4-4*x_2*x_3^3)+(x_2^4*x_3^5-2*x_2^5*x_3^4+x_\
2^6*x_3^3)*y_1-x_2^4*x_3^5+2*x_2^5*x_3^4-x_2^6*x_3^3)/(x_1^2*(-x_2^4*x_3^5+2*x_2^5*x_3^4-x_2^6*x_3^3)+x_1^3*(2*x_2^3*x_3^5-3*x_2^4*x_3^4+x_2^6*x_3^2)+x_1^4*(-x_2^2*x_3^5+3*x_2^4*x_3^3-2*x_2^5*x_3^2)+\
x_1^5*(x_2^2*x_3^4-2*x_2^3*x_3^3+x_2^4*x_3^2))]
(%i7) 
Run Example
linel: 500;
(%o1)                                                                                                                                                                                                                                                   500
(%i2) eq1:a*x_1^5 + b*x_1^4+ c*x_1^3+ d*x_1^2  + 1 = y_1;
                                                                                                                                                                                                                                         5        4        3        2
(%o2)                                                                                                                                                                                                                               a x_1  + b x_1  + c x_1  + d x_1  + 1 = y_1
(%i3) eq2:a*x_2^5 + b*x_2^4+ c*x_2^3+ d*x_2^2 + 1 = y_2;
                                                                                                                                                                                                                                         5        4        3        2
(%o3)                                                                                                                                                                                                                               a x_2  + b x_2  + c x_2  + d x_2  + 1 = y_2
(%i4) eq3:5*a*x_2^4 + 4*b*x_2^3+ 3*c*x_2^2+ 2*d*x_2  = 0;
                                                                                                                                                                                                                                           4          3          2
(%o4)                                                                                                                                                                                                                               5 a x_2  + 4 b x_2  + 3 c x_2  + 2 d x_2 = 0
(%i5) eq4:a*x_3^5 + b*x_3^4+ c*x_3^3+ d*x_3^2 + 1 = y_3;
                                                                                                                                                                                                                                         5        4        3        2
(%o5)                                                                                                                                                                                                                               a x_3  + b x_3  + c x_3  + d x_3  + 1 = y_3
(%i6) string(linsolve([eq1,eq2,eq3,eq4],[a,b,c,d]));
(%o6) [a = (x_1^3*(x_2^4*(2*y_3-2)+(2*x_3^4-4*x_2^2*x_3^2)*y_2-2*x_3^4+4*x_2^2*x_3^2)+x_1^2*(x_2^5*(1-y_3)+(4*x_2^2*x_3^3-3*x_2*x_3^4)*y_2+3*x_2*x_3^4-4*x_2^2*x_3^3)+x_1^4*(x_2^3*(1-y_3)+(3*x_2*x_3^2-2*x_3^3)*y_2+2*x_3^3-3*x_2*x_3^2)+(x_2^3*x_3^4-2*x_2^4*x_3^3+x_2^5*x_3^2)*y_1-x_2^3*x_3^4+2*x_2^4*x_3^3-x_2^5*x_3^2)/(x_1^2*(-x_2^5*x_3^5+2*x_2^6*x_3^4-x_2^7*x_3^3)+x_1^3*(2*x_2^4*x_3^5-3*x_2^5*x_3^4+x_2^7*x_3^2)+x_1^4*(-x_2^3*x_3^5+3*x_2^5*x_3^3-2*x_2^6*x_3^2)+x_1^5*(x_2^3*x_3^4-2*x_2^4*x_3^3+x_2^\
5*x_3^2)),b = -(x_1^3*(x_2^5*(3*y_3-3)+(2*x_3^5-5*x_2^3*x_3^2)*y_2-2*x_3^5+5*x_2^3*x_3^2)+x_1^5*(x_2^3*(1-y_3)+(3*x_2*x_3^2-2*x_3^3)*y_2+2*x_3^3-3*x_2*x_3^2)+x_1^2*(x_2^6*(2-2*y_3)+(5*x_2^3*x_3^3-3*x_2*x_3^5)*y_2+3*x_2*x_3^5-5*x_2^3*x_3^3)+(x_2^3*x_3^5-3*x_2^5*x_3^3+2*x_2^6*x_3^2)*y_1-x_2^3*x_3^5+3*x_2^5*x_3^3-2*x_2^6*x_3^2)/(x_1^2*(-x_2^5*x_3^5+2*x_2^6*x_3^4-x_2^7*x_3^3)+x_1^3*(2*x_2^4*x_3^5-3*x_2^5*x_3^4+x_2^7*x_3^2)+x_1^4*(-x_2^3*x_3^5+3*x_2^5*x_3^3-2*x_2^6*x_3^2)+x_1^5*(x_2^3*x_3^4-2*x_2^4*\
x_3^3+x_2^5*x_3^2)),c = (x_1^4*(x_2^5*(3*y_3-3)+(2*x_3^5-5*x_2^3*x_3^2)*y_2-2*x_3^5+5*x_2^3*x_3^2)+x_1^2*(x_2^7*(1-y_3)+(5*x_2^3*x_3^4-4*x_2^2*x_3^5)*y_2+4*x_2^2*x_3^5-5*x_2^3*x_3^4)+x_1^5*(x_2^4*(2-2*y_3)+(4*x_2^2*x_3^2-2*x_3^4)*y_2+2*x_3^4-4*x_2^2*x_3^2)+(2*x_2^4*x_3^5-3*x_2^5*x_3^4+x_2^7*x_3^2)*y_1-2*x_2^4*x_3^5+3*x_2^5*x_3^4-x_2^7*x_3^2)/(x_1^2*(-x_2^5*x_3^5+2*x_2^6*x_3^4-x_2^7*x_3^3)+x_1^3*(2*x_2^4*x_3^5-3*x_2^5*x_3^4+x_2^7*x_3^2)+x_1^4*(-x_2^3*x_3^5+3*x_2^5*x_3^3-2*x_2^6*x_3^2)+x_1^5*(x_2\
^3*x_3^4-2*x_2^4*x_3^3+x_2^5*x_3^2)),d = -(x_1^4*(x_2^5*(2*y_3-2)+(3*x_3^5-5*x_2^2*x_3^3)*y_2-3*x_3^5+5*x_2^2*x_3^3)+x_1^3*(x_2^6*(1-y_3)+(5*x_2^2*x_3^4-4*x_2*x_3^5)*y_2+4*x_2*x_3^5-5*x_2^2*x_3^4)+x_1^5*(x_2^4*(1-y_3)+(4*x_2*x_3^3-3*x_3^4)*y_2+3*x_3^4-4*x_2*x_3^3)+(x_2^4*x_3^5-2*x_2^5*x_3^4+x_2^6*x_3^3)*y_1-x_2^4*x_3^5+2*x_2^5*x_3^4-x_2^6*x_3^3)/(x_1^2*(-x_2^4*x_3^5+2*x_2^5*x_3^4-x_2^6*x_3^3)+x_1^3*(2*x_2^3*x_3^5-3*x_2^4*x_3^4+x_2^6*x_3^2)+x_1^4*(-x_2^2*x_3^5+3*x_2^4*x_3^3-2*x_2^5*x_3^2)+x_1^5*\
(x_2^2*x_3^4-2*x_2^3*x_3^3+x_2^4*x_3^2))]
(%i7) 
Run Example
linel:4096;
(%o1)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         4096
(%i2) stardisp:true;
(%o2)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         true
(%i3) display2d:false;

(%o3) false
(%i4) 
linsolve([ay1*t1+by1=ay2*t2+by2,ax1*t1+bx1=ax2*t2+bx2],[t1,t2]);

(%o4) [t1 = -(ax2*by2-ax2*by1+ay2*(bx1-bx2))/(ax1*ay2-ax2*ay1),t2 = -(ax1*by2-ax1*by1+ay1*(bx1-bx2))/(ax1*ay2-ax2*ay1)]
(%i5) 

Related Help

Help for Linel