Finding the determinant without det(a)
21 次查看(过去 30 天)
显示 更早的评论
I need to physically computate the determinant of 3x3 matrices and lower without using the det(a) function, rather loops. What's causing this error?
1 个评论
DGM
2022-1-27
编辑:DGM
2022-1-27
Like Max, I'm guessing that matrix a1 isn't the expected size. In the future, it's best to post the code as formatted code instead of pictures so that people can troubleshoot it.
I'm just going to throw this out there. You can reuse the code for calculating the 2x2 case.
A = rand(3);
DA1 = det(A)
DA2 = det33(A)
function D = det22(A)
D = prod(A([1 4])) - prod(A([2 3]));
end
function D = det33(A)
D = A(1,1)*det22(A(2:3,2:3)) - A(1,2)*det22(A(2:3,[1 3])) + A(1,3)*det22(A(2:3,1:2));
end
采纳的回答
Max Heimann
2022-1-27
编辑:Max Heimann
2022-1-27
According to the error you are trying to access an index which is not part of the variable.
Are you sure a1 is a 3x3 matrix?
It looks like you are calling the function with a 2x3 matrix [1 2 3; 4 5 6]
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!