Rosetta Code

Arithmetic/Complex: Difference between revisions

Page Discussion

Arithmetic/Complex (view source)

Revision as of 21:41, 18 January 2021

398 bytes removed ,  3 years ago
→‎{{header|Phix}}: use builtin
(→‎{{header|Phix}}: use builtin)
Line 3,780:
 
=={{header|Phix}}==
<lang Phix>constant-- REAL = 1,demo\rosetta\ArithComplex.exw
include complex.e
IMAG = 2
complex a = complex_new(1,1),
 
b = complex_new(3.14159,1.25),
type complex(sequence s)
sc = sprintfcomplex_new("%g"1,ar0),▼
return length(s)=2 and atom(s[REAL]) and atom(s[IMAG])
d = complex_new(0,1)
end type
printf(1,"a = %s\n",{scomplexcomplex_sprint(a)})▼
 
printf(1,"b = %s\n",{scomplexcomplex_sprint(b)})▼
function add(complex a, complex b)
printf(1,"1/ac = %s\n",{scomplexcomplex_sprint(inv(a)c)})▼
return sq_add(a,b)
printf(1,"a+bd = %s\n",{scomplexcomplex_sprint(add(a,b)d)})▼
end function
printf(1,"a*+b = %s\n",{scomplexcomplex_sprint(multcomplex_add(a,b))})▼
 
printf(1,"-a*b = %s\n",{scomplexcomplex_sprint(negcomplex_mul(a,b))})</lang>▼
function mult(complex a, complex b)
printf(1,"1/a = %s\n",{complex_sprint(complex_inv(a))})
return {a[REAL] * b[REAL] - a[IMAG] * b[IMAG],
printf(1,"c/a = %s\n",{complex_sprint(complex_div(c,a))})
a[REAL] * b[IMAG] + a[IMAG] * b[REAL]}
printf(1,"c-a = %s\n",{complex_sprint(complex_sub(c,a))})
end function
printf(1,"d-a = %s\n",{complex_sprint(complex_sub(d,a))})
 
printf(1,"-a = %s\n",{complex_sprint(complex_neg(a))})
function inv(complex a)
printf(1,"conj a = %s\n",{complex_sprint(complex_conjugate(a))})</lang>
atom denom
denom = a[REAL] * a[REAL] + a[IMAG] * a[IMAG]
return {a[REAL] / denom, -a[IMAG] / denom}
end function
 
function neg(complex a)
return sq_uminus(a)
end function
 
function scomplex(complex a)
sequence s = ""
atom ar, ai
{ar, ai} = a
if ar!=0 then
▲ s = sprintf("%g",ar)
end if
if ai!=0 then
if ai=1 then
s &= "+i"
elsif ai=-1 then
s &= "-i"
else
s &= sprintf("%+gi",ai)
end if
end if
if length(s)=0 then
return "0"
end if
return s
end function
 
complex a, b
a = { 1.0, 1.0 }
b = { 3.14159, 1.2 }▼
▲printf(1,"a = %s\n",{scomplex(a)})
▲printf(1,"b = %s\n",{scomplex(b)})
▲printf(1,"a+b = %s\n",{scomplex(add(a,b))})
▲printf(1,"a*b = %s\n",{scomplex(mult(a,b))})
▲printf(1,"1/a = %s\n",{scomplex(inv(a))})
▲printf(1,"-a = %s\n",{scomplex(neg(a))})</lang>
{{out}}
<pre>
a = 1+i
b = 3.14159+1.2i25i
c = 1
a+b = 4.14159+2.2i
d = i
a*b = 1.94159+4.34159i▼
▲a+b = { 34.14159, 1+2.2 }25i
▲a*b = 1.9415989159+4.34159i39159i
1/a = 0.5-0.5i
c/a = 0.5-0.5i
c-a = -i
d-a = -1
-a = -1-i
conj a = 1-i
</pre>
 
Petelomax
7,806

edits

Retrieved from "https://rosettacode.org/wiki/Special:MobileDiff/63667"