Question: Why the frist row of my matrix is fixed, while I expect it varies?

In the following program why the first row of Matrix P is costant, while I expect it varies?


 

``

> 

``

> 

restart

> 

with(LinearAlgebra):

> 

MCK := Matrix(1, 1, {(1, 1) = 0.1627682387e-16*mu})

> 

P := Matrix(2, 111):

> 

`ρ__∞` := 1.225:

> 

j := 1:

> 

for `U__&infin;` from 333 to 335 do `M__&infin;` := `U__&infin;`/(331.2); mu := `&rho;__&infin;`*`U__&infin;`*(`M__&infin;`^2-2)/sqrt((`M__&infin;`^2-1)^3); P[1 .. 2, j] := `<,>`(`<,>`(MCK), `<,>`(j)); j := j+1 end do:

> 

P[1, 1 .. j-1]

Vector[row]([0.1627682387e-16*mu, 0.1627682387e-16*mu, 0.1627682387e-16*mu])

(1)
> 

``


 

Download Sooal.mw

Please Wait...