Programming with SCILAB (pdf)

or,
--> surf(x,y,z’)
9.7
0.0
-9.7
-2.50
0.00
2.50 10
5
0
Nested loops can be used for other operations with matrices. For example, a function
that multiplies two matrices is shown next:
function [C] = matrixmult(A,B)
//|--------------------------------------|
//| This function calculates the matrix |
//| C(nxm) = A(nxp)*B(pxm) |
//|--------------------------------------|
//First, check matrices compatibility
[nrA,ncA] = size(A)
[nrB,ncB] = size(B)
if ncA nrB then
error('matrixmult - incompatible matrices A and B')
abort
end
//Calculate matrix product
C = zeros(nrA,ncB)
for i = 1:nrA
for j = 1:ncB
for k = 1:ncA
C(i,j) = C(i,j) + A(i,k)*B(k,j)
end
end
end
An application of function matrixmult is shown next. Here we multiply A 3x4 with B 4x6
resulting in C 3x6 . However, the product B 4x6 times A 3x4 does not work.
-->A = int(10*rand(3,4))
A =
! 2. 3. 8. 0. !
! 7. 6. 6. 5. !
! 0. 6. 8. 6. !
-->B = int(10*rand(4,6))
B =
! 7. 2. 3. 3. 3. 2. !
! 1. 2. 9. 2. 5. 6. !
! 5. 8. 2. 5. 5. 4. !
! 2. 6. 3. 4. 4. 9. !
-->C = matrixmult(A,B)
19