Vector products: Difference between revisions
Content added Content deleted
No edit summary |
(→{{header|Julia}}: A new entry for Julia) |
||
Line 1,422: | Line 1,422: | ||
"a . (b x c) = 6 )" |
"a . (b x c) = 6 )" |
||
"a x (b x c) = [-267, 204, -3]"</lang> |
"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}}== |
=={{header|Liberty BASIC}}== |