[..] 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]