package org.matheclipse.core.reflection.system;

import static org.matheclipse.core.expression.F.*;

import org.matheclipse.core.eval.interfaces.AbstractFunctionEvaluator;
import org.matheclipse.core.eval.interfaces.AbstractTrigArg1;
import org.matheclipse.core.eval.interfaces.INumeric;
import org.matheclipse.core.expression.ComplexUtils;
import org.matheclipse.core.expression.ComplexNum;
import org.matheclipse.core.expression.Num;
import org.matheclipse.core.expression.F;
import org.matheclipse.core.interfaces.IAST;
import org.matheclipse.core.interfaces.IExpr;
import org.matheclipse.core.interfaces.ISymbol;
import org.matheclipse.parser.client.SyntaxError;

/**
 * Arccosine
 * 
 * See 
 * Inverse_trigonometric functions
 */
public class ArcCos extends AbstractTrigArg1 implements INumeric {
	// final static String[] RULES = {
	// "ArcCos[0]=Pi/2", "ArcCos[1]=0", "ArcCos[I]=Pi/2-I*ArcSinh[1]",
	// "ArcCos[2^(1/2)/2]=Pi/4", "ArcCos[3^(1/2)/2]=Pi/6",
	// "ArcCos[DirectedInfinity[1]]=I*Infinity",
	// "ArcCos[1/2*(2+2^(1/2))^(1/2)]=1/8*Pi",
	// "ArcCos[1/2*(2-2^(1/2))^(1/2)]=3/8*Pi",
	// "ArcCos[1/4*5^(1/2)+1/4]=1/5*Pi",
	// "ArcCos[1/4*6^(1/2)*(1+1/3*3^(1/2))]=1/12*Pi",
	// "ArcCos[1/4*2^(1/2)*(5-5^(1/2))^(1/2)]=3/10*Pi",
	// "ArcCos[1/4*6^(1/2)*(1-1/3*3^(1/2))]=5/12*Pi", "ArcCos[Cosh[1]]=I",
	// "ArcCos[1/4*5^(1/2)-1/4]=2/5*Pi",
	// "ArcCos[1/4*2^(1/2)*(5+5^(1/2))^(1/2)]=1/10*Pi", "ArcCos[1/2]=1/3*Pi",
	// };

	/**
	 * 
	 * ArcCos[1/2*2^(1/2)]=1/4*Pi,
	 *   ArcCos[1/2*3^(1/2)]=1/6*Pi,
	 *   ArcCos[1/2*(2+2^(1/2))^(1/2)]=1/8*Pi,
	 *   ArcCos[1/2*(2-2^(1/2))^(1/2)]=3/8*Pi,
	 *   ArcCos[(1+1/3*3^(1/2))*1/4*6^(1/2)]=1/12*Pi,
	 *   ArcCos[(1-1/3*3^(1/2))*1/4*6^(1/2)]=5/12*Pi,
	 *   ArcCos[-1/4+1/4*5^(1/2)]=2/5*Pi,
	 *   ArcCos[1]=0,
	 *   ArcCos[0]=1/2*Pi,
	 *   ArcCos[1/2]=1/3*Pi,
	 *   ArcCos[I]=1/2*Pi-I*ArcSinh[1],
	 *   ArcCos[DirectedInfinity[1]]=I*Infinity,
	 *   ArcCos[(5+5^(1/2))^(1/2)*1/4*2^(1/2)]=1/10*Pi,
	 *   ArcCos[(5-5^(1/2))^(1/2)*1/4*2^(1/2)]=3/10*Pi,
	 *   ArcCos[Cosh[1]]=I,
	 *   ArcCos[1/4+1/4*5^(1/2)]=1/5*Pi,
	 * 
*/ final static IAST RULES = List( Set(ArcCos(Times(Times(C1D4, Power(C2, C1D2)), Power(Plus(Times(CN1, Power(C5, C1D2)), C5), C1D2))), Times(fraction(3L, 10L), Pi)), Set(ArcCos(Times(Times(C1D4, Plus(Times(C1D3, Power(C3, C1D2)), C1)), Power(integer(6L), C1D2))), Times(fraction( 1L, 12L), Pi)), Set(ArcCos(Times(Times(C1D4, Plus(Times(CN1D3, Power(C3, C1D2)), C1)), Power(integer(6L), C1D2))), Times( fraction(5L, 12L), Pi)), Set(ArcCos(Times(Times(C1D4, Power(C2, C1D2)), Power(Plus(Power(C5, C1D2), C5), C1D2))), Times( fraction(1L, 10L), Pi)), Set(ArcCos(Times(C1D2, Power(C2, C1D2))), Times(C1D4, Pi)), Set(ArcCos(Times(C1D2, Power(C3, C1D2))), Times(fraction(1L, 6L), Pi)), Set(ArcCos(Times(C1D2, Power(Plus(Power(C2, C1D2), C2), C1D2))), Times(fraction( 1L, 8L), Pi)), Set(ArcCos(Times(C1D2, Power(Plus(Times(CN1, Power(C2, C1D2)), C2), C1D2))), Times(fraction(3L, 8L), Pi)), Set(ArcCos(C1), C0), Set(ArcCos(C0), Times(C1D2, Pi)), Set(ArcCos(C1D2), Times(C1D3, Pi)), Set(ArcCos(CI), Plus(Times(Times( CN1, CI), ArcSinh(C1)), Times(C1D2, Pi))), Set(ArcCos(CInfinity), Times(CI, CInfinity)), Set(ArcCos(Cosh(C1)), CI), Set( ArcCos(Plus(Times(C1D4, Power(C5, C1D2)), C1D4)), Times(fraction(1L, 5L), Pi)), Set(ArcCos(Plus(Times(C1D4, Power(C5, C1D2)), Times(CN1, C1D4))), Times(fraction(2L, 5L), Pi))); @Override public IAST getRuleAST() { return RULES; } public ArcCos() { } @Override public IExpr evaluateArg1(final IExpr arg1) { if (AbstractFunctionEvaluator.isNegativeExpression(arg1)) { return Plus(Times(CN1, Pi), ArcCos(Times(CN1, arg1))); } return null; } @Override public IExpr numericEvalD1(final Num arg1) { return F.num(Math.acos(arg1.getRealPart())); } @Override public IExpr numericEvalDC1(final ComplexNum arg1) { return ComplexUtils.acos(arg1); } public double evalReal(final double[] stack, final int top, final int size) { if (size != 1) { throw new UnsupportedOperationException(); } return Math.acos(stack[top]); } @Override public void setUp(final ISymbol symbol) throws SyntaxError { symbol.setAttributes(ISymbol.LISTABLE | ISymbol.NUMERICFUNCTION); super.setUp(symbol); } }

Ads help maintain this website.