clc;

clear;

close all;

S=0;

k=input('k=please inter a number');

i=input('i=please inter a number');

j=input('j=please inter a number');

for s=1:.1:k

    for m=1:2:i

        for n=1:2:j

            W=(2.911650856*(.32*m^2+.32*n^2*s^2+11.2))*sin((1/2)*m*pi)*sin((1/2)*n*pi)/(m*n*(n^2*s^2+m^2)^2*((0.2e-1*(n^2*s^2+m^2))*(pi^4-96)+68.18636376));



        S=S+W;

      

        end

    end

    wbar=0.01*S;

    plot(s,wbar);

    hold on;
end