Integrering i Python
Lagt inn: 04/03-2013 13:30
Hei. Sitter med ei oppgave her eg ikkje får til. har fått funksjonen
f(t)=exp(-t)*sin(2*Pi*t)
Oppgava går ut på å lage eit plot av den integrerte funksjonen og den dobbel integrerte funksjonen.
Eg har laget dette programmet:
from scitools.all import *
def dydt(t):
return (exp(-t)*sin(2*pi*t))
a=0.0
b=4.0
N=50
dt=(b-a)/N
time=linspace(a,b,N+1)
yTrue=exp(-time)*sin(2*pi*time)
y=zeros(len(time))
y[0]=0
for i in range(0,len(time)-1,1):
y[i+1]=y+dt*dydt(time)
plot(time,yTrue,'k',time,y,'k--')
---------------------------------------------------------------------------------
Og når eg da plotter dette, så får eg opp plot for den funksjonen eg har gitt.(sjekket med maple)
Noen inteligente der ute som kan hjelpe meg videre? Står bom fast her, og kan på ingen måte skryte på meg noe kunnskap med Python.
på Forhånd takk
-Tor Olav-
f(t)=exp(-t)*sin(2*Pi*t)
Oppgava går ut på å lage eit plot av den integrerte funksjonen og den dobbel integrerte funksjonen.
Eg har laget dette programmet:
from scitools.all import *
def dydt(t):
return (exp(-t)*sin(2*pi*t))
a=0.0
b=4.0
N=50
dt=(b-a)/N
time=linspace(a,b,N+1)
yTrue=exp(-time)*sin(2*pi*time)
y=zeros(len(time))
y[0]=0
for i in range(0,len(time)-1,1):
y[i+1]=y+dt*dydt(time)
plot(time,yTrue,'k',time,y,'k--')
---------------------------------------------------------------------------------
Og når eg da plotter dette, så får eg opp plot for den funksjonen eg har gitt.(sjekket med maple)
Noen inteligente der ute som kan hjelpe meg videre? Står bom fast her, og kan på ingen måte skryte på meg noe kunnskap med Python.
på Forhånd takk
-Tor Olav-