- scl.jdg [c+]
[..] La fonction [jdg] betascl
Code
int beta_jdg(){
scltracefa(§, ƒ, ∅);
scltracefc("scl-%s (%s) : %s\n\n", sclver(), sclsec(), temps_char(0));
int an=1001, ai;
double *a=newtablin(an, -3, 3); //[-3..3]
int bn=1001, bi;
double *b=newtablin(bn, -3, 3); //[-3..3]
double **z=newmat<double>(an, bn);
int **zin=newmat<int>(an, bn);
double cb, za=∞rd, ze=εrd;
for(ai=0;ai<an;ai++){
for(bi=0;bi<bn;bi++){
cb = beta(a[ai], b[bi]);
if(cb==0) { z[ai][bi] = 0; }
else if(cb>0){ z[ai][bi] = log10(fabs(cb)); }
else { z[ai][bi] =-log10(fabs(cb)); }
zin[ai][bi] = 1;
if(estentier(a[ai])>0){ if(a[ai]<=0){ zin[ai][bi] = 0; } } //📝
if(estentier(b[bi])>0){ if(b[bi]<=0){ zin[ai][bi] = 0; } } //📝
if(zin[ai][bi]>0){
if(za>z[ai][bi]){ za = z[ai][bi]; }
if(ze<z[ai][bi]){ ze = z[ai][bi]; }
}
}
}
// Dynamique symetrique
ze = std::max(ze, -za); za=-ze;
scl::string_c imgname=scl::string_c(ƒ)+".pdf",
imgpath="./srt/"+imgname;
grafmat_dis(imgpath.c_str(), an, bn, z,
"x", a, "y", b, "instatus", zin,
/*"grid", "on",*/ "axis", "equal",
"title", "Beta function",
"subtitlex", imgname.c_str(), /*"subtitle", subtitle.c_str(),*/
/*"linespeca", "-a",*/ "colormap", "rhiver"/*"rqhult"*//*"rlcecd"*/,
"background-color", "w", /*"logaxis", "z",*/
"xlabel", "a", "ylabel", "b",
"azmax", ze, "azmin", za, ∅);
scltracefe(§, ƒ, ∅);
return 0;
}
/* ƒ décorée par
le 23-05-2026 19:11:34 */
Sortie
[>..\xpl\src\scl.jdg.cpp.beta_jdg]
scl-26.05 (gwin64) : 23-05-2026 19:11:34
<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
<< END OF DISLIN / VERSION 11.3.2 <<
<< Date : 23.05.2026 Time : 19:11:37 Pageformat: DA4L <<
<< Vectors : 84 Warnings: 0 Fileformat: PDF <<
<< Metafile: ./srt/beta_jdg.pdf <<
<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[<..\xpl\src\scl.jdg.cpp.beta_jdg]
