Anonymous user
Vector products: Difference between revisions
→{{header|Julia}}: A new entry for Julia
No edit summary |
(→{{header|Julia}}: A new entry for Julia) |
||
Line 1,422:
"a . (b x c) = 6 )"
"a x (b x c) = [-267, 204, -3]"</lang>
=={{header|Julia}}==
Julia provides dot and cross products as built-ins. It's easy enough to use these to construct the triple products.
<lang Julia>
function scltrip{T<:Number}(a::AbstractArray{T,1},
b::AbstractArray{T,1},
c::AbstractArray{T,1})
dot(a, cross(b, c))
end
function vectrip{T<:Number}(a::AbstractArray{T,1},
b::AbstractArray{T,1},
c::AbstractArray{T,1})
cross(a, cross(b, c))
end
a = [3, 4, 5]
b = [4, 3, 5]
c = [-5, -12, -13]
println("Test Vectors:")
println(" a = ", a)
println(" b = ", a)
println(" c = ", a)
println("\nVector Products:")
println(" a dot b = ", dot(a, b))
println(" a cross b = ", cross(a, b))
println(" a dot b cross c = ", scltrip(a, b, c))
println(" a cross b cross c = ", vectrip(a, b, c))
</lang>
{{out}}
<pre>
Test Vectors:
a = [3,4,5]
b = [3,4,5]
c = [3,4,5]
Vector Products:
a dot b = 49
a cross b = [5,5,-7]
a dot b cross c = 6
a cross b cross c = [-267,204,-3]
</pre>
=={{header|Liberty BASIC}}==
| |||