Code covered by the BSD License  

Highlights from
Constitutive Models of Linear Viscoelasticity with Laplace Transform

Constitutive Models of Linear Viscoelasticity with Laplace Transform

by

 

The program computes compliance and relaxation function.

[G]=gauss(A)
function [G]=gauss(A)
%gauss elimination
[m,n]=size(A);
j=0;	%column
for k=1:m-1;	%row diagonal
    %urceni prvniho clenu
    j=j+1;
    if j>n
        break
    end
    index=find(A(k+1:m,j)~=0);
    if (A(k,j)==0)&&isequal(index,[])
        continue
    elseif A(k,j)==0
        A_pom=A(k,:);        
        A(k,:)=A(k+index(1),:);
        A(k+index(1),:)=A_pom;
        for i=k+1:m %row
            A(i,:)=A(i,:)-A(i,j)/A(k,j)*A(k,:);
        end
    else
        for i=k+1:m %row
            A(i,:)=A(i,:)-A(i,j)/A(k,j)*A(k,:);
        end
    end
end
G=A;
end

Contact us