## load values
load expm_data

## test matrix
disp ("Test matrix:");
disp (tmpa);

## calculation
disp ("Octave build in expm:");
disp (expm (tmpa));

## expected matrix
disp ("Expected matrix:");
disp (expa);

## 2nd norm
disp ("Error norm:");
disp (norm (expa - expm (tmpa), 2));
