diff options
Diffstat (limited to 'charf.c')
| -rw-r--r-- | charf.c | 27 |
1 files changed, 26 insertions, 1 deletions
@@ -220,15 +220,40 @@ void format_result(char* s, int index, int k, EXPTYPE p, VALUETYPE r, int format root_to_string(rts,r,p); l += sprintf(s+l,"|f^(%2d)|_%s: %s",k,e,rts); /*Compute and print also Mf.*/ + /*Compute and print also df and dMf.*/ /* VALUETYPE Mf[N]; compute_maximalfunction(f, Mf, d); - l += sprintf(s+l,"\nMf: "); + l += sprintf(s+l,"\nM: "); for(int i=0; i<d; i++){ char v[8]; valuetype_to_string(v,Mf[i]); l += sprintf(s+l,"%s ",v); } + VALUETYPE df[2][N]; + VALUETYPE dMf[2][N]; + for(int i=0; i<=1; i++){ + for(int n=0; n<d; n++){ + df[i][n] = f[n]; + dMf[i][n] = Mf[n]; + } + } + for(int m=0; m+1<=k; m++){ + differentiate(df[(m+1)%2],df[m%2],d,1); + differentiate(dMf[(m+1)%2],dMf[m%2],d,1); + l += sprintf(s+l,"\nf%d: ",m+1); + for(int i=0; i<d; i++){ + char v[8]; + valuetype_to_string(v,df[(m+1)%2][i]); + l += sprintf(s+l,"%s ",v); + } + l += sprintf(s+l,"\nM%d: ",m+1); + for(int i=0; i<d; i++){ + char v[8]; + valuetype_to_string(v,dMf[(m+1)%2][i]); + l += sprintf(s+l,"%s ",v); + } + } */ } else if(format == FORMAT_LATEX){ |