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]](/img/smilies/smile.gif)
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)
Testa att använda en zygohistomorfisk prepromorfism istället för en for-loop, det borde lösa problemet.
Gentlernen:
zygohistomorfisk prepromorfism
Haha shit, vad är det? 😀 vi la in en tom vektor som sedan fylldes i på slutet! det fixade det för oss 🙂 tack för hjälpen ändå <3
