蟒蛇（.m 文件） 你能按照程序代碼做這個功能嗎？ use python （.m file） Can you follow the procedure code to do this function?
Posted: Thu May 26, 2022 9:15 am
蟒蛇（.m 文件）
你能按照程序代碼做這個功能嗎？
use python （.m file）
Can you follow the procedure code to do this function?
Procedure Gauss(n, a, b) • real: (a)nxn, (b)nx1, (x)nx1, rmax, smax, xmult integer: n, lnx1, P for i from 1 to n do • l; <i endfor for k=1 to n-1 do ●rmax < 0 • for ik to n do @ smax (0 • for j = k to n do smax • endfor @ rmax +|a,k/smax| • pri • endif max(smax, |ae, j) if rmax < |ae,,k/smax| •endfor
% Elinimation Step for ik + 1 to n do .xmultae,,k/ alkk • alik 0 for j=k+1 to n do • ajaj - (xmult)ax.j • endfor • be = be, - (xmult) bex • endfor endfor • Xn = ben/ann • for i= n-1 to 1 step -1 do • sum <- be for j =i+1 to n do • sumsum - aj X Xj • endfor • X; (sum)/ai endfor
Procedure Gauss(n, a, b) e real: (a)nxn, (b)nx1, (x)nx1, rmax, smax, xmult integer: n, lnx1, P • for i from 1 to n do • endfor • for k=1 to n-1 do ●rmax < 0 for ik to n do a smax (0 • for j = k to n do • smax + max(smax, la,j) • endfor ← if rmax < |a.k/smax| ●rmax lat,,k/smax • pri endif • endfor
% Elinimation Step for ik 1 to n do xmultak/ak • 31₁ k 40 for j k + 1 to n do = • ajaj(xmult) ak • endfor • bebe (xmult) bex endfor endfor • Xn= ben/an.n for in-1 to 1 step -1 do • sum <- be for j =i+1 to n do sum-sum - aj x xj • endfor X; (sum)/a, endfor
- 1 (22.14 KiB) Viewed 12 times
- 2 (17.57 KiB) Viewed 12 times
use python （.m file）
- 3 (17.88 KiB) Viewed 12 times
- 4 (14.3 KiB) Viewed 12 times
Procedure Gauss(n, a, b) • real: (a)nxn, (b)nx1, (x)nx1, rmax, smax, xmult integer: n, lnx1, P for i from 1 to n do • l; <i endfor for k=1 to n-1 do ●rmax < 0 • for ik to n do @ smax (0 • for j = k to n do smax • endfor @ rmax +|a,k/smax| • pri • endif max(smax, |ae, j) if rmax < |ae,,k/smax| •endfor
% Elinimation Step for ik + 1 to n do .xmultae,,k/ alkk • alik 0 for j=k+1 to n do • ajaj - (xmult)ax.j • endfor • be = be, - (xmult) bex • endfor endfor • Xn = ben/ann • for i= n-1 to 1 step -1 do • sum <- be for j =i+1 to n do • sumsum - aj X Xj • endfor • X; (sum)/ai endfor
Procedure Gauss(n, a, b) e real: (a)nxn, (b)nx1, (x)nx1, rmax, smax, xmult integer: n, lnx1, P • for i from 1 to n do • endfor • for k=1 to n-1 do ●rmax < 0 for ik to n do a smax (0 • for j = k to n do • smax + max(smax, la,j) • endfor ← if rmax < |a.k/smax| ●rmax lat,,k/smax • pri endif • endfor
% Elinimation Step for ik 1 to n do xmultak/ak • 31₁ k 40 for j k + 1 to n do = • ajaj(xmult) ak • endfor • bebe (xmult) bex endfor endfor • Xn= ben/an.n for in-1 to 1 step -1 do • sum <- be for j =i+1 to n do sum-sum - aj x xj • endfor X; (sum)/a, endfor