<![CDATA[
double mylog(double x) {
return log(x);
}
int main() {
array double x[3][3]={1,2,3,
3,4,5,
6,7,8};
array double y[3][3];
expm(y,x);
printf("x = \n%f",x);
printf("y = \n%f",y);
funm(x,mylog,y);
printf("x = \n%f",x);
}
]]>
<![CDATA[
x =
1.000000 2.000000 3.000000
3.000000 4.000000 5.000000
6.000000 7.000000 8.000000
y =
157372.953093 200034.605129 242697.257164
290995.910241 369883.552084 448769.193928
491431.845963 624654.472518 757878.099072
x =
1.000000 2.000000 3.000000
3.000000 4.000000 5.000000
6.000000 7.000000 8.000000
]]>