/dports/math/giacxcas/giac-1.6.0/src/ |
H A D | intgab.cc | 209 f1=normal(_texpand(f1,contextptr),contextptr); in is_even_odd() 210 f2=normal(_texpand(f2,contextptr),contextptr); in is_even_odd() 493 dp=gcd(re(dp,contextptr),im(dp,contextptr),contextptr); in intgab_sincos() 761 gen reg=simplify(re(g,contextptr),contextptr); in intgab_r() 977 …po,contextptr)*pow(b-a,1+na*expo,contextptr)*Gamma(inv(na,contextptr),contextptr)*Gamma(expo+1,con… in intgab() 1009 …res=g0mult*pow(va1x,va2,contextptr)*pow(-vb1x,vb2,contextptr)*pow(b-a,va2+vb2+1,contextptr)*Beta(v… in intgab() 1015 …res=g0mult*pow(-va1x,va2,contextptr)*pow(vb1x,vb2,contextptr)*pow(b-a,va2+vb2+1,contextptr)*Beta(v… in intgab() 1310 g=_lin(trig2exp(g0,contextptr),contextptr); in intgab() 1348 radius=exp(re(beta,contextptr),contextptr); in intgab() 1477 tmp1=ratnormal(subst(tmp1,gx,pow(gx,inv(d,contextptr),contextptr),false,contextptr),contextptr); in in_sumab_int() [all …]
|
H A D | desolve.cc | 299 …it->num,l,contextptr),r2e(firstcoeff(a),lprime,contextptr),contextptr)*exp(r2e(rdiv(-b,a,contextpt… in pf_ilaplace() 312 sqrtdelta=normal(sqrt(r2e(delta,lprime,contextptr),contextptr),contextptr); in pf_ilaplace() 319 sqrtdelta=normal(sqrt(r2e(-delta,lprime,contextptr),contextptr),contextptr); in pf_ilaplace() 438 if (!is_linear_wrt(eval(ftmp,eval_level(contextptr),contextptr),t,a,b,contextptr)) in is_linear_diffeq() 922 t=simplify(exp(integrate_without_lnabs(u,x,contextptr)/2,contextptr),contextptr); in desolve_linn() 952 w=simplify(w/(u*derive(V,x,contextptr)-V*derive(u,x,contextptr)),contextptr); in desolve_linn() 1391 f=normal((derive(M,y,contextptr)-derive(N,x,contextptr))/N,contextptr); in desolve_f() 1404 f=normal((derive(N,x,contextptr)-derive(M,y,contextptr))/M,contextptr); in desolve_f() 1881 a=_eval(sqrt(1+4*a,contextptr),contextptr); in kovacic_iscase3() 1937 return ratnormal(_lin(_eval(e,contextptr),contextptr),contextptr); in explnsimp() [all …]
|
H A D | intg.cc | 383 return exp(ln(c,contextptr)/n,contextptr); in surd() 431 c=-pow(-c,inv(n,contextptr),contextptr); in makelnatan() 433 c=pow(c,inv(n,contextptr),contextptr); in makelnatan() 444 c=pow(c,inv(n,contextptr),contextptr); in makelnatan() 446 …,X,c,false,contextptr),contextptr)*lnabs2(X-c,X,contextptr)+normal_norootof(subst(N,X,-c,false,con… in makelnatan() 454 c=pow(-c,inv(n,contextptr),contextptr); in makelnatan() 645 c=pow(r2e(-*it/v.front(),lprime,contextptr),inv(d,contextptr),contextptr); in integrate_deno_length_2() 689 c=pow(c,inv(n,contextptr),contextptr); in integrate_deno_length_2() 1651 …lnpart=lnpart+rdiv(r2e(b,l,contextptr),r2e(a,l,contextptr),contextptr)*lnabs(r2e(it->fact,l,contex… in integrate_rational_end() 2028 return x*atan(x,contextptr)-rdiv(ln(abs(pow(x,2)-1,contextptr),contextptr),plus_two,contextptr); in int_atanh() [all …]
|
H A D | ti89.cc | 133 return rdiv(ln(e1,contextptr),ln(b,contextptr),contextptr)+n; in _logb() 265 …p(makevecteur(string2gen(eval(v[0],eval_level(contextptr),contextptr).print(contextptr)),v[1],1),c… in _InputStr() 1378 return _tran(_cumSum(_tran(v.front(),contextptr),contextptr),contextptr); in _cumSum() 1789 return evalf(_floor(g,contextptr),eval_level(contextptr),contextptr); in _iPart() 1791 return -evalf(_floor(-g,contextptr),eval_level(contextptr),contextptr); in _iPart() 1793 return _iPart(re(g,contextptr),contextptr)+cst_i*_iPart(im(g,contextptr),contextptr); in _iPart() 1834 return sto(_mRow(args,contextptr),v1,contextptr); in _mRow() 1986 gen args=evalf(eval(g,1,contextptr),eval_level(contextptr),contextptr); in _QR() 2741 return exact(re(g,contextptr),contextptr)+cst_i*exact(im(g,contextptr),contextptr); in exact() 2896 arc[g.print(contextptr)]=eval(g,eval_level(contextptr),contextptr); in _Archive() [all …]
|
H A D | subst.cc | 146 gen a=re(e,contextptr),b=im(e,contextptr); in exp2sincos() 147 …return exp(a,contextptr)*(cos(radtoangle(b,contextptr),contextptr)+cst_i*sin(radtoangle(b,contextp… in exp2sincos() 1446 tmp[s]=_floor(precision*evalf(re(ln(w[i],contextptr),contextptr),1,contextptr),contextptr); in ext_relation() 1447 tmp[s+1]=_floor(precision*evalf(im(ln(w[i],contextptr),contextptr),1,contextptr),contextptr); in ext_relation() 1607 …anch_evalf(rdiv(ln(r2e(a_orig,vars,contextptr),contextptr)-res,cst_pi_over_2*cst_i,contextptr),con… in expanded_ln() 1646 res=res+ln(r2sym(n,l,contextptr)/r2sym(d,l,contextptr),contextptr); in simplifylnarg() 1731 return exp(tmp1*_pow2exp(ln(v[0],contextptr),contextptr),contextptr); in _pow2exp() 2166 return _prod(expln2trig(inv(ff,contextptr),contextptr),contextptr); in expln2trig() 2561 …gen g1=normal(trigcos(reg,contextptr),contextptr)+cst_i*normal(trigcos(img,contextptr),contextptr); in simplify() 2562 …gen g2=normal(trigsin(reg,contextptr),contextptr)+cst_i*normal(trigsin(img,contextptr),contextptr); in simplify() [all …]
|
H A D | signalprocessing.cc | 867 return ln(_abs(x+b/a,contextptr),contextptr)+ln(_abs(a,contextptr),contextptr); in logabs_expand() 983 ret=ret*exp(cst_i*im(ex,contextptr),contextptr)*exp(re(ex,contextptr),contextptr); in lin_cplx() 990 gen f=factorise(lin_cplx(g,contextptr)+lin_cplx(h,contextptr),contextptr); in hcollect() 1061 … contextptr),contextptr),contextptr),contextptr),x,contextptr),contextptr),contextptr); in fourier_terms() 1632 gen g=factorise(_lin(expand(g_orig,contextptr),contextptr),contextptr); in is_expsum() 1806 is_zero(im(b,contextptr)) && is_positive(re(-a,contextptr),contextptr)) { in fourier_func_type() 2011 (_inv(_abs(s,contextptr),contextptr)+2*cst_euler_gamma*_Dirac(s,contextptr)); in fourier() 2073 ret+=_exp2pow(_lin(recursive_normal(t,contextptr),contextptr),contextptr); in fourier() 2135 ret+=_exp2pow(_lin(recursive_normal(t,contextptr),contextptr),contextptr); in fourier() 2411 gen expr=exp(alpha*ln(sin(k*_IDNT_pi()/(N-1),contextptr),contextptr),contextptr); in _cosine_window() [all …]
|
H A D | plot.cc | 206 return ratnormal(_lin(g*conj(g,contextptr),contextptr),contextptr); in abs_norm2() 502 …if (is_greater(eps*abs(pm,contextptr)*abs(wv,contextptr),scalar_product(pm,wv,contextptr),contextp… in find_best() 2815 return normal(im(g,contextptr)/re(g,contextptr),contextptr); in _slope() 4869 gen x=(b1-a1).re(contextptr),y=(b1-a1).im(contextptr),vx=v.re(contextptr),vy=v.im(contextptr); in inter2droites2() 5001 gen surface=re(AB,contextptr)*im(AC,contextptr)-im(AB,contextptr)*re(AC,contextptr); in inscrit() 5195 gen res1=re(a,contextptr)*re(b,contextptr)+im(a,contextptr)*im(b,contextptr); in normalized_scalar_product() 5239 …if ( is_zero(normal(subst(im(ee[0],contextptr),tt,*it,false,contextptr)-im(f,contextptr),contextpt… in on() 5930 res += im(cote2,contextptr)*re(cote1,contextptr)-re(cote2,contextptr)*im(cote1,contextptr); in _aire() 6622 …contextptr),contextptr) && is_zero(derive(fxy,x,contextptr),contextptr) && is_zero(derive(fyy,x,co… in equation2geo2d() 10035 …contextptr),false,contextptr)-im(M,contextptr),contextptr),test2=recursive_normal(subst(xt,var,re(… in tangent() [all …]
|
H A D | derive.cc | 108 return derive(*iti,i,contextptr)+derive(*(iti+1),i,contextptr); in derive_SYMB() 134 return (derive(u,i,contextptr)*v-u*derive(v,i,contextptr))/pow(v,2,contextptr); in derive_SYMB() 209 gen res=-derive(eq,x,contextptr)/derive(eq,y,contextptr); in derive_SYMB() 320 gen tmp=pow(v2,v0-1,contextptr)*pow(1-v2,v1-1,contextptr)*derive(v2,i,contextptr); in derive_SYMB() 345 gen tmp1=exp(-v1,contextptr)*pow(v1,v0-1,contextptr)*derive(v1,i,contextptr); in derive_SYMB() 498 return s.sommet(derive(s.feuille,i,contextptr),contextptr); in derive_SYMB() 607 eprime=ratnormal(derive(ecopie,vars,contextptr),contextptr); in derive() 684 res=ratnormal(derive(res,w[i],contextptr),contextptr); in _derive() 826 gen g=simplify(subst(d2,var,v[i],false,contextptr),contextptr); in critical() 948 gen dg(derive(g(_tmp,contextptr),_tmp,contextptr)); in _function_diff() [all …]
|
H A D | series.cc | 70 series_flags(contextptr)=series_flags(contextptr) | (1<<7) ; in taylor() 448 gen c_max = min(normal(a_min + b_max,contextptr),normal(b_min + a_max,contextptr),contextptr); in pmul() 1122 return pmul(res,normalize_sqrt(pow(first.coeff,e,contextptr),contextptr),res,contextptr); in ppow() 1441 tmp=normal(trig2exp(*lv1_it,contextptr),contextptr); in series__SPOL1() 1711 der=pow(x,a-1,contextptr)*exp(-x,contextptr); in series__SPOL1() 2528 gen e=Heavisidetosign(when2sign(piecewise2when(e0,contextptr),contextptr),contextptr); in in_limit() 2625 e_copy=_pow2exp(tan2sincos(exact(e,contextptr),contextptr),contextptr); in in_limit() 2954 faster_var_subst.push_back(pow(w,c,contextptr)*exp(normal(f-c*g,contextptr),contextptr)); in mrv_lead_term() 3005 p=subst(p,ln(w,contextptr),g,false,contextptr); in mrv_lead_term() 3007 p=subst(p,ln(w,contextptr),-g,false,contextptr); in mrv_lead_term() [all …]
|
H A D | global.cc | 278 if (contextptr && contextptr->globalptr ) in last_evaled_function_name() 286 if (contextptr && contextptr->globalptr ) in currently_scanned() 294 if (contextptr && contextptr->globalptr ) in last_evaled_argptr() 302 if (contextptr && contextptr->globalptr ) in language() 308 if (contextptr && contextptr->globalptr ) in language() 350 if (contextptr && contextptr->globalptr) in default_color() 380 if (contextptr && contextptr->globalptr) in spread_Row() 394 if (contextptr && contextptr->globalptr) in spread_Col() 408 if (contextptr && contextptr->globalptr) in printcell_current_row() 422 if (contextptr && contextptr->globalptr) in printcell_current_col() [all …]
|
H A D | mathml.cc | 473 +gen2mathml(eval(v[0],eval_level(contextptr),contextptr),contextptr); 476 +gen2mathml(eval(v[0],eval_level(contextptr),contextptr),contextptr); 1080 s=s+re(v[i],contextptr).print(contextptr)+" "+im(v[i],contextptr).print(contextptr)+" "; 1088 …contextptr).print(contextptr)+" "+im(v[i],contextptr).print(contextptr)+" , "+re(v[i],contextptr)… 1090 s=s+re(v[i],contextptr).print(contextptr)+" "+im(v[i],contextptr).print(contextptr)+" "; 1298 if (is_zero(eval(re(A,contextptr)-re(B,contextptr),eval_level(contextptr),contextptr))){ 1302 …gen a=eval((im(A,contextptr)-im(B,contextptr))/(re(A,contextptr)-re(B,contextptr)),eval_level(cont… 1303 gen b=eval(im(A,contextptr)-a*re(A,contextptr),eval_level(contextptr),contextptr); 1312 A=evalf(eval(A,1,contextptr),1,contextptr); B=evalf(eval(B,1,contextptr),1,contextptr); 1358 …s = s+re(evalf(v[i],1,contextptr),contextptr).print(contextptr)+" "+im(evalf(v[i],1,contextptr),co… [all …]
|
H A D | gauss.cc | 192 u1=recursive_normal(rdiv(derive(q,x[r],contextptr),plus_two,contextptr),contextptr); in gauss() 239 l1=rdiv(derive(q,x[r1],contextptr),2,contextptr); in gauss() 240 l2=rdiv(derive(q,x[r2],contextptr),2,contextptr); in gauss() 571 gen svp0(exact(sign(vp0,contextptr),contextptr)), in conique_reduite() 572 svp1(exact(sign(vp1,contextptr),contextptr)), in conique_reduite() 573 scoeffcst(exact(sign(coeffcst,contextptr),contextptr)); in conique_reduite() 722 P[0]=normal(_normalize(P[0],contextptr),contextptr); in quadrique_reduite() 723 P[1]=normal(_normalize(P[1],contextptr),contextptr); in quadrique_reduite() 814 gen a=sqrt(-d/s1,contextptr),b=sqrt(-d/s2,contextptr),c=sqrt(-d/s3,contextptr); in quadrique_reduite() 828 gen a=sqrt(d/s1,contextptr),b=sqrt(d/s2,contextptr),c=sqrt(-d/s3,contextptr); in quadrique_reduite() [all …]
|
H A D | moyal.cc | 155 return ratnormal(res*cst_pi/sin(cst_pi*a_,contextptr)/Gamma(n,contextptr),contextptr); in Beta() 158 return Gamma(a,contextptr)*Gamma(b,contextptr)/Gamma(a+b,contextptr); in Beta() 635 …return Gamma(nd+1,contextptr)/Gamma(kd+1,contextptr)/Gamma(nd-kd+1,contextptr)*pow(pd,kd,contextpt… in binomial() 653 return comb(n,k,contextptr)*pow(p,k,contextptr)*pow(1-p,n-k,contextptr); in binomial() 1068 return exp(-m,contextptr)*pow(m,k,contextptr)/_factorial(k,contextptr); in poisson() 1284 …+1,2,contextptr),contextptr)/Gamma(rdiv(n,2,contextptr),contextptr)/sqrt(n*cst_pi,contextptr)*pow(… in student() 1534 …return rdiv(pow(x,n2-1,contextptr)*exp(-x/2,contextptr),Gamma(n2,contextptr)*pow(2,n2,contextptr),… in chisquare() 1717 …return pow(a/b,a/2,contextptr)/Beta(a/2,b/2,contextptr) * pow(x,a/2-1,contextptr) * pow(1+a/b*x,-(… in snedecor() 2004 return k/lambda*pow(tmp,k-1,contextptr)*exp(-pow(tmp,k,contextptr),contextptr); in weibull() 2225 …return pow(x,alpha-1,contextptr)*exp(-beta*x,contextptr)*pow(beta,alpha,contextptr)/Gamma(alpha,co… in gammad() [all …]
|
H A D | usual.cc | 416 gen res=ln(abs(e,contextptr),contextptr)+cst_i*arg(e,contextptr); in ln() 516 return (ln(abs(e,contextptr),contextptr)+cst_i*arg(e,contextptr))/ln(10,contextptr); in log10() 541 return rdiv(ln(e1,contextptr),ln(10,contextptr),contextptr)+n; in log10() 1270 a=re(e,contextptr);b=ratnormal(im(e,contextptr),contextptr); in sqrt() 1479 e1= cos(e1,contextptr)*cosh(e2,contextptr)-cst_i*sinh(e2,contextptr)*sin(e1,contextptr); in cos() 1721 gen res=sin(e1,contextptr)*cosh(e2,contextptr)+cst_i*sinh(e2,contextptr)*cos(e1,contextptr); in sin() 7320 res=_floor(re(res,contextptr)+.5,contextptr)+cst_i*_floor(im(res,contextptr)+.5,contextptr); in _evalf() 9140 …erfc=evalf(2*inv(sqrt(cst_pi,contextptr),contextptr),1,contextptr)*exp(-inv(z*z,contextptr),contex… in erf0() 9668 …return _Ci(args,contextptr)*sin(args,contextptr)+(evalf(cst_pi/2,1,contextptr)-_Si(args,contextptr… in _SiCi_f() 9712 …return -_Ci(args,contextptr)*cos(args,contextptr)+(evalf(cst_pi/2,1,contextptr)-_Si(args,contextpt… in _SiCi_g() [all …]
|
H A D | lin.cc | 313 tadd(res,newcoeff,iscos?cos(tmp2,contextptr):sin(tmp1,contextptr),contextptr); in tconvolution() 314 tadd(res,newcoeff,iscos?cos(tmp1,contextptr):sin(tmp2,contextptr),contextptr); in tconvolution() 334 tadd(res,newcoeff,sin(tmp1,contextptr),contextptr); in tconvolution() 335 tadd(res,newcoeff,sin(tmp2,contextptr),contextptr); in tconvolution() 338 tadd(res,newcoeff,cos(tmp2,contextptr),contextptr); in tconvolution() 339 tadd(res,-newcoeff,cos(tmp1,contextptr),contextptr); in tconvolution() 594 …lifier(prod_expand(prod_expand(it,it+s/2,contextptr),prod_expand(it+s/2,itend,contextptr),contextp… in prod_expand() 657 *it=pow(expand(*it,contextptr),v[1],contextptr); in expand_pow_expand() 736 return pow(exp(arg,contextptr),coeff,contextptr); in exp_expand() 795 …return sin_expand(first,contextptr)*cos_expand(last,contextptr)+cos_expand(first,contextptr)*sin_e… in sin_expand() [all …]
|
H A D | rpn.cc | 620 …if (!is_zero(tmp.eval(eval_level(contextptr),contextptr).evalf(eval_level(contextptr),contextptr))) in _RPN_UNTIL() 959 if (contextptr->globalcontextptr!=contextptr){ in _purge() 1418 gen g=eval(args,eval_level(contextptr),contextptr); in _ADDROW() 1523 return sto(_SWAPCOL(g,contextptr),v0,contextptr); in _SWAPCOL() 1619 return _ranm(eval(args,eval_level(contextptr),contextptr),contextptr); in _RANDMAT() 1813 gen res=_egvl(evalf(args,1,contextptr),contextptr); in _EIGENVAL() 1970 v[i]=eval(v[i],eval_level(contextptr),contextptr); in hp38_eval() 2274 …return makevecteur(eval(a,eval_level(contextptr),contextptr),eval(b, eval_level(contextptr),contex… in _IS_LINEAR() 2471 gen gf=_floor(log10(abs(g,contextptr),contextptr),contextptr); in _XPON() 2768 value=evalf(subst(v0,v1,value,false,contextptr),eval_level(contextptr),contextptr); in _ITERATE() [all …]
|
H A D | optimization.cc | 193 …vecteur e_orig_simp=*expexpand(expand(_pow2exp(e_orig,contextptr),contextptr),contextptr)._VECTptr; in solve2() 208 vecteur e(*halftan(_texpand(hyp2exp(e_orig_simp,contextptr),contextptr),contextptr)._VECTptr); in solve2() 242 r[i]=simp(ln(c,contextptr),contextptr); in solve2() 2128 while (is_strictly_greater(_abs(c-d,contextptr),epsilon(contextptr),contextptr)) { in find_peak() 2335 eps=exact(epsilon(contextptr)/2,contextptr); in tprob() 2916 tp+=simp(c*ak,contextptr)*cos(simp((j+1)*twopi/T,contextptr)*x,contextptr); in triginterp() 2917 tp+=simp(c*bk,contextptr)*sin(simp((j+1)*twopi/T,contextptr)*x,contextptr); in triginterp() 3209 gen lambda=exp(_inv(k,contextptr)*ln(s0/n,contextptr),contextptr); in weibull_mle() 3293 gen e=ln(aidn,contextptr)-Psi(aidn,contextptr)-ln(mean,contextptr)+slog/N; in _fitdistr() 3641 parm,pvals,false,contextptr),contextptr),contextptr); in _conjugate_equation() [all …]
|
H A D | plot3d.cc | 259 …gen M=P+u__IDNT_e*(cos(theta,contextptr)*n+sin(theta,contextptr)*(cos(v__IDNT_e,contextptr)*n1+sin… in cone() 263 gen tmp=evalf_double(v[3]/cos(theta,contextptr),eval_level(contextptr),contextptr); in cone() 472 eqs=evalf_double(sign(eq,contextptr),1,contextptr); in polyedre() 735 gen v1(normal(b-O,contextptr)),v2(normal(c-O,contextptr)),v3(normal(d-O,contextptr)); in _tetraedre_centre() 737 gen B = normal(O - v1/3 - sqrt(2,contextptr)*v2/3 - sqrt(6,contextptr)*v3/3,contextptr); in _tetraedre_centre() 738 gen C = normal(O - v1/3 - sqrt(2,contextptr)*v2/3 + sqrt(6,contextptr)*v3/3,contextptr); in _tetraedre_centre() 771 gen D1 = normal(A+x-sqrt(2,contextptr)/3*(C-A)+sqrt(6,contextptr)/3*(D-A),contextptr); in _cube_centre() 1560 gen fx(derive(f,x,contextptr)),fy(derive(f,y,contextptr)),fz(derive(f,z,contextptr)); in equation2geo3d() 1564 …contextptr)),fxy(derive(fx,y,contextptr)),fyy(derive(fy,y,contextptr)),fxz(derive(fx,z,contextptr)… in equation2geo3d() 1578 …fx=normal(fx/d,contextptr); fy=normal(fy/d,contextptr); fz=normal(fz/d,contextptr); c=normal(c/d,c… in equation2geo3d() [all …]
|
H A D | solve.cc | 1336 testeq=abs(evalf(subst(e,x,m,false,contextptr),eval_level(contextptr),contextptr),contextptr); in in_solve_inequation() 1389 …lncollect(lnexpand(ln(simplify(a1,contextptr),contextptr)-ln(simplify(a2,contextptr),contextptr),c… in solve_inequation() 2800 arg1=ln(simplify(a1,contextptr),contextptr)-ln(simplify(a2,contextptr),contextptr); in _solve_uncompressed() 2812 …arg1=lnexpand(ln(simplify(a1,contextptr),contextptr)-ln(simplify(a2,contextptr),contextptr),contex… in _solve_uncompressed() 2907 …gen n=ln(abs(b-a,contextptr),contextptr)-ln(max(abs(b,contextptr),abs(a,contextptr),contextptr),co… in bisection_solver_sr() 4020 …if (is_greater(1e-8,abs(im(gguess,contextptr),contextptr)/abs(re(gguess,contextptr),contextptr),co… in in_fsolve() 4485 fa=evalf(eval(subst(f,x,a,false,contextptr),eval_level(contextptr),contextptr),1,contextptr); in newton() 4513 fa=evalf(eval(subst(f,x,a,false,contextptr),eval_level(contextptr),contextptr),1,contextptr); in newton() 4520 fa=evalf(eval(subst(f,x,a,false,contextptr),eval_level(contextptr),contextptr),1,contextptr); in newton() 4555 fb=evalf(eval(subst(f,x,b,false,contextptr),eval_level(contextptr),contextptr),1,contextptr); in newton() [all …]
|
H A D | misc.cc | 339 return inv(cos(args,contextptr),contextptr); in _sec() 347 return inv(sin(args,contextptr),contextptr); in _csc() 355 return rdiv(cos(args,contextptr),sin(args,contextptr),contextptr); in _cot() 363 return acos(inv(args,contextptr),contextptr); in _asec() 371 return asin(inv(args,contextptr),contextptr); in _acsc() 467 … return epsilon2zero(re(g,contextptr),contextptr)+cst_i*epsilon2zero(im(g,contextptr),contextptr); in epsilon2zero() 1484 …return _float2rational(re(g,contextptr),contextptr)+cst_i*_float2rational(im(g,contextptr),context… in _float2rational() 3286 …return f[1]*sqrt(Gamma(1+gen(2)/f[0],contextptr)-pow(Gamma(1+gen(1)/f[0],contextptr),2,contextptr)… in _stddev() 5413 gen A(re(ab,contextptr)),B(im(ab,contextptr)),C(re(ac,contextptr)),D(im(ac,contextptr)); in cross_prod() 7563 gen f1=_factor(derive(f,t,contextptr),contextptr),g1=_factor(derive(g,t,contextptr),contextptr); in step_param_() [all …]
|
H A D | prog.cc | 418 s += w[0].print(contextptr)+";"+w[1].print(contextptr)+";"+w[2].print(contextptr)+")"; in debug_print() 961 res += indent(contextptr)+it->print(contextptr); in printasprogram() 1743 if (!equaltosto(*it,contextptr).in_eval(eval_level(contextptr),res,contextptr)) in _bloc() 1939 …test=equaltosame(test.eval(eval_level(contextptr),contextptr)).eval(eval_level(contextptr),context… in ifte() 2936 contextptr=contextptr->previous; in leave() 3323 tmpsto=sto(gen(er,contextptr),id,contextptr); in _try_catch() 3350 tmpsto=sto(gen(er,contextptr),id,contextptr); in _try_catch() 7782 f=exact(gen(s,contextptr),contextptr); in _convert() 8227 return eval(quote_read(args,contextptr),eval_level(contextptr),contextptr); in _read() 13048 …test=equaltosame(test.eval(eval_level(contextptr),contextptr)).eval(eval_level(contextptr),context… in _piecewise() [all …]
|
H A D | sym2poly.cc | 803 gen q=evalf(b,1,contextptr)/evalf(num,1,contextptr); in sym2rxrootnum() 2040 return rdiv(r2sym(-f.num,l,contextptr),r2sym(-f.den,l,contextptr),contextptr); in r2sym() 2299 …(p._FRACptr->den,contextptr),lt,ltend,contextptr),r2sym(p._FRACptr->den*conj(p._FRACptr->den,conte… in r2sym() 2397 if (is_greater(abs(tmp3-tmp00,contextptr),abs(tmp2-tmp00,contextptr),contextptr)) in r2sym() 2531 return symb_rootof(r2sym(*(pp._EXTptr),l,contextptr),r2sym(f,l,contextptr),contextptr); in r2sym() 2664 res=sqrt(a,contextptr)+cs*sqrt(b,contextptr); in is_sqrtxy() 2805 …contextptr)*pow(doubl,expnum,contextptr)*pow(recursive_normal(r2e(s,lv,contextptr),contextptr),nov… in in_normalize_sqrt() 2962 if (contextptr && contextptr->previous) in normal() 3232 *it= pow(base,rdiv(num,den,contextptr),contextptr); in recursive_normal() 3265 res=simplifier(ratnormal(normalize_sqrt(e_copy,contextptr),contextptr),contextptr); in recursive_normal() [all …]
|
H A D | alg_ext.cc | 132 g=_factors(makesequence(g,rootof(g,contextptr)),contextptr); in galoisconj() 173 if (is_greater(borne*abs(R1,contextptr),abs(R1-R2,contextptr)*pow10,contextptr)){ in conj_in_nf() 252 …uble max_re(re(current,contextptr).evalf_double(1,contextptr)._DOUBLE_val),max_im(im(current,conte… in in_select_root() 255 …double cur_re(re(*it,contextptr).evalf_double(1,contextptr)._DOUBLE_val),cur_im(im(*it,contextptr)… in in_select_root() 280 gen r=in_select_root(a,is_real(v,contextptr),contextptr); in select_root() 286 …return rdiv(alg_evalf(a._FRACptr->num,b,contextptr),alg_evalf(a._FRACptr->den,b,contextptr),contex… in alg_evalf() 287 gen a1=a.evalf(1,contextptr),b1=b.evalf(1,contextptr); in alg_evalf() 549 K=_floor(re(K,contextptr),contextptr)+1; in common_minimal_POLY() 579 gen r2=_max(real_proot(v2,1e-12,contextptr),contextptr); in common_minimal_POLY() 1288 fraction fa(e2r(exact(g,contextptr),l,contextptr)); in sturm() [all …]
|
H A D | risch.cc | 59 g=tsimplify(pow2expln(g,x,contextptr),contextptr); in risch_tower() 252 prim=r2sym(R,lv,contextptr)*prim+r2sym(V,lv,contextptr)/r2sym(c,lv,contextptr); in SPDE() 261 prim=r2sym(N1,lv,contextptr)/r2sym(N2,lv,contextptr); in SPDE() 447 prim=tsimplify(exp(prim,contextptr),contextptr); in risch_desolve() 545 gen r=e2r(coeff,l,contextptr),ar=e2r(a,l,contextptr); in integrate_ei() 629 gen dX=ratnormal(derive(X,x,contextptr),contextptr); in risch_poly_part() 766 …al(allowed_lnarg*r2sym(it->num,l,contextptr)/(derive(allowed_lnarg,x,contextptr)*r2sym(it->den,l,c… in in_risch() 803 …d(makevecteur(recursive_normal(tmp1,contextptr),recursive_normal(r2sym(it->fact,l,contextptr),cont… in in_risch() 848 tmpcst=exp(tmpcst,contextptr)*exp(ratnormal(tmp-tmpcst,contextptr),contextptr); in risch_lin() 889 lin(trig2exp(e_orig,contextptr),vexp,contextptr); in risch() [all …]
|
/dports/math/giacxcas/giac-1.6.0/doc/en/ |
H A D | giac_bassin.cpp | 19 N_i = cpp_convert_2 (Ng, contextptr); in bassin() 21 P1 = giac::_diff (P, contextptr); in bassin() 22 P = giac::_symb2poly (P, contextptr); in bassin() 23 P1 = giac::_symb2poly (P1, contextptr); in bassin() 24 res_v = cpp_convert_7 (giac::_makelist ((N_i + 1) * (N_i + 1), contextptr), contextptr); in bassin() 25 r_v = cpp_convert_7 (giac::_proot (P, contextptr), contextptr); in bassin() 28 r[i]=cpp_convert_4(r_v[i],contextptr); in bassin() 30 …cpp_convert_1 (giac::_division (giac::makesequence (xmax + (-xmin), N_i), contextptr), contextptr); in bassin() 31 …cpp_convert_1 (giac::_division (giac::makesequence (ymax + (-ymin), N_i), contextptr), contextptr); in bassin() 34 …rt_4 (giac::_evalf (xmin + j_i * hx_d + gen (0, 1) * (ymin + k_i * hy_d), contextptr), contextptr); in bassin() [all …]
|