clear all, close all, clc
h=6;
B1=8;
B2=3;
V=h/3*(B1+sqrt(B1*B2)+B2)

baser=[h B1 B2];
fel=[0.3 0.2 0.1];
totalt_slutfel=0;
for i =1:3
    baser(i)=baser(i)+fel(i);
    Vfel=baser(1)/3*(baser(2)+sqrt(baser(2)...
        *baser(3))+baser(3));
    totalt_slutfel=totalt_slutfel+abs(Vfel-V);

    baser=[h B1 B2];
end
format short
disp(['volymen är ' num2str(V) ' plus/minus ' num2str(totalt_slutfel)]);