ich habe eine Frage zum Rechteckverfahren (numerische Integration). Aufgabe ist es, eine R Funktion zu schreiben, getArea(lower,upper,n), die mit approximativ über Rechtecke die Fläche unter der Funktion y = x^2 im Intervall [lower,upper] schätzen soll. n ist dabei die Anzahl der Rechtecke zur Schätzung.
Meine Lösung dazu:
Code: Alles auswählen
foo=function(x) {
y=x^2
return(y)
}
getArea = function(lower,upper,n) {
for (i in 0:n-1) {
yi = foo(lower + i * (upper-lower)/n + (upper-lower)/(2*n))
return(yi)
}
z = ((upper-lower)/n) * sum(yi[1:length(yi)])
return(z)
}
Meine Frage: Wo liegt der Fehler?
Bin über jede Hilfe dankbar.
MfG
ein Anfänger