Trapetsregeln
Jag försöker använda mig av trapetsregeln för att beräkna integralen av sqrt(x-4)dx mellan -1 och 1 för h=2
Problemet är bara att det blir ett irrationellt tal vilket det inte ska bli. Siffrorna i det irrationella talet stämmer potentiellt över vad jag velat få från den reella delen. Är det någon som ser vad felet är? Jag gör inte :-(
Program:
clear all, close all, clc
f=inline('sqrt(x-4)');
a=-1; b=1; h=2;
N=(b-a)/h ;
x=a+[1:N-1]*h;
T0=(f(a)+f(b))/2;
T=h*(T0+sum(f(x)))
Ger:
T= 0.0000 + 3.9681i
Är du säker på att funktionen är ? Den funktionen är endast (reellt) definierad för x ≥ 4.
pepparkvarn skrev:Är du säker på att funktionen är ? Den funktionen är endast (reellt) definierad för x ≥ 4.
De har du ju rätt i!!! Ska kolla när jag kommer hem!
pepparkvarn skrev:Är du säker på att funktionen är ? Den funktionen är endast (reellt) definierad för x ≥ 4.
det var jag som skrivit fel! ska vara sqrt(x+4). Tack! såg inte det själv. Allt är rätt nu
Louiger skrev:pepparkvarn skrev:Är du säker på att funktionen är ? Den funktionen är endast (reellt) definierad för x ≥ 4.
det var jag som skrivit fel! ska vara sqrt(x+4). Tack! såg inte det själv. Allt är rätt nu
:)
