Ecuaciones Diferenciales de primer orden

·         Archivo M

function  dy=g1(t,y)

dy=funcion en términos t y Y;

·         Para cualquier tipo de ode:

[t y]=ode113(@g1,[valor inferior de t:espaciado h:valor inferior de t],valor inicial];

%En los cmondos anteriores puede ser cualquier tipo de ode

·         Para la solución exacta

Exacta=dsolve(‘Dy=funcion en términos de y y t’,’Y()=valor inicial’)

Y=subs(Exacta,t)

·         Error

Error=norm(Y-y,inf)

Ecuación de orden superior de 2 ecuaciones  

%En este caso se hará un proceso a través de método del anulador  para obtener la funcion DY Y DX

Archivo M

function p=f1(t,w)

p=zeros(size(w));

p=[numerox 0;0 numeroy]*w+[funcion restante de dx;function restante de dy];

·         Para cualquier tipo de Ode

t=[valor inferior:espaciado h;valor superior];

[t xy]=ode45(@f1,t,[valor X inicial;valor Y inicial])

·         Exacto

[x y]=dsolve(‘funcion DX','Funcion DY','x(0)=0','y(0)’)

ExactoX=subs(x,t);

ExactoY=subs(y,t);

Exacto=[ExactoX ExactoY]

errorx=norm(ExactoX-xy(:,1));

errory=norm(ExactoY-xy(:,2));

error=[errorx errory]