Forumet - Hjälp med Matlab-algoritm, funkar inte

Hjälp med Matlab-algoritm, funkar inte

194 0 2
Hej! Vi sitter med ett projekt i Beräkningsvetenskap där vi ska skriva en algoritm för att beräkna tryck i en damm som beskrivs genom en integral. Vi ska använda trapets-metoden och sedan plotta trycket (I1 här nedan).

Om någon har tid, ta gärna en titt på vår programmering och se om du hittar något fel! Tack så jättemycket på förhand!! [smile]

p=62.5;
V = [10:100];
for D = 10:100
n = 0;
h = 4; %steglängd
while n==0;
X1 = [0:h:D];
Q1 = [0:2*h:D];
X2 = p.*(D-X1).*(40-20.*exp((-0.01.*X1).^2));
Q2 = p.*(D-Q1).*(40-20.*exp((-0.01.*Q1).^2));
I1 = trapz(X1,X2);
I2 = trapz(Q1,Q2);
fel= abs((I1-I2)/3);
if (fel)>=(0.5e-4);
h = h/2; %Om felet är för stort så halveras steglängden
else n = 1; % När felet är tillräckligt litet anges trycket
end %if
end % whileloop
I(D-9)= I1;
I = [I1];
end % forloop
plot(V,I)