درخواست اجرای کد در یک سیستم پرسرعت!

[ASedJavad]

سلام
یه کد نوشتم که رو سیستم من هنگ میکنه و اجرا نمیشه.
بیزحمت اگه کسی سیستم خوبی داره، یه تست بکنه و اگه تونست جواب بگیره، خروجی اون رو برام بفرسته (فایل ASedJavad.txt)
کد رو در ادامه مینویسم، اما m فایل مربوطه رو هم ضمیمه میکنم.

clear all
syms k57 k58 k59 k60 k61 k62 k63 k64 k65 k66 k67 k68 k69;
syms k70 k71 k72 k73 k74 k75 k76 k77 k78 k79 k80 k81 k82 ;
syms k83 k84 k85 k86 k87 k88 k89 k90 k91 k92 k93 k94 k95 ;
syms k96 k97 k98 k99 k100 k101 k102 k103 k104 k105 k106;
syms x g1 g2 g3 g4 g5 g6 g7 g8 g9 g10;
g = [g5 , g4 , g3 , g2 , g1, 0 , 0 , 0
g10, g9 , g8 , g7 , g6, 0 , 0 , 0
0 , g5 , g4 , g3 , g2, g1, 0 , 0
0 , g10, g9 , g8 , g7, g6, 0 , 0
0 , 0 , g5 , g4 , g3, g2, g1, 0
0 , 0 , g10, g9 , g8, g7, g6, 0
0 , 0 , 0 , g5 , g4, g3, g2, g1
0 , 0 , 0 , g10, g9, g8, g7, g6];
g = det(g);
g0 = [k57 * x^4 + k58 * x^3 + k59 * x^2 + k60 * x + k61
k62 * x^4 + k63 * x^3 + k64 * x^2 + k65 * x + k66
k67 * x^4 + k68 * x^3 + k69 * x^2 + k70 * x + k71
k72 * x^4 + k73 * x^3 + k74 * x^2 + k75 * x + k76
k77 * x^4 + k78 * x^3 + k79 * x^2 + k80 * x + k81
k82 * x^4 + k83 * x^3 + k84 * x^2 + k85 * x + k86
k87 * x^4 + k88 * x^3 + k89 * x^2 + k90 * x + k91
k92 * x^4 + k93 * x^3 + k94 * x^2 + k95 * x + k96
k97 * x^4 + k98 * x^3 + k99 * x^2 + k100 * x + k101
k102 * x^4 + k103 * x^3 + k104 * x^2 + k105 * x + k106];
g = subs(g,[g1;g2;g3;g4;g5;g6;g7;g8;g9;g10],g0);
g0 = collect(g,x);
fid = fopen('ASedJavad.txt', 'wt');
fprintf(fid, '%s\n', char(g0));
fclose(fid);

[rahnema1]

یه نرم افزار هست به اسم Maxima که رایگان هم هست و کم حجم کد شما به زبان maxima اینجور میشه ممکنه محاسبه طول بکشه اما سیستم هنگ نمی کنه

g : matrix([g5 , g4 , g3 , g2 , g1, 0 , 0 , 0],
[g10, g9 , g8 , g7 , g6, 0 , 0 , 0],
[0 , g5 , g4 , g3 , g2, g1, 0 , 0],
[0 , g10, g9 , g8 , g7, g6, 0 , 0],
[0 , 0 , g5 , g4 , g3, g2, g1, 0],
[0 , 0 , g10, g9 , g8, g7, g6, 0],
[0 , 0 , 0 , g5 , g4, g3, g2, g1],
[0 , 0 , 0 , g10, g9, g8, g7, g6]
)$
gd:determinant(g)$
gs: psubst( [g1=k57 * x^4 + k58 * x^3 + k59 * x^2 + k60 * x + k61 ,
g2=k62 * x^4 + k63 * x^3 + k64 * x^2 + k65 * x + k66 ,
g3=k67 * x^4 + k68 * x^3 + k69 * x^2 + k70 * x + k71 ,
g4=k72 * x^4 + k73 * x^3 + k74 * x^2 + k75 * x + k76 ,
g5=k77 * x^4 + k78 * x^3 + k79 * x^2 + k80 * x + k81 ,
g6=k82 * x^4 + k83 * x^3 + k84 * x^2 + k85 * x + k86 ,
g7=k87 * x^4 + k88 * x^3 + k89 * x^2 + k90 * x + k91 ,
g8=k92 * x^4 + k93 * x^3 + k94 * x^2 + k95 * x + k96 ,
g9=k97 * x^4 + k98 * x^3 + k99 * x^2 + k100 * x + k101,
g10=k102 * x^4 + k103 * x^3 + k104 * x^2 + k105 * x + k106] ,gd)$


gc:facsum (expand(gs),x)$
stringout ("c:\\ASedJavad.txt",gc)