Diferencia entre revisiones de «Reacciones complejas - Grupo 1 A»

De MateWiki
Saltar a: navegación, buscar
(Página creada con «{{ TrabajoED | Reacciones complejas. Grupo 1A | Ecuaciones Diferenciales|Curso 2014-15 | María Ramírez Ig...»)
 
Línea 7: Línea 7:
  
  
 +
{{matlab|codigo=
  
 +
%PVI por Euler con una h = 0.1
 +
clear all
 +
% y'= cos(t) - y
 +
%condiciones dadas
 +
y0 = 0;
 +
t0 = 0;
 +
tN = 2;
 +
h = 0.1;
 +
%calculamos N
 +
N = (tN - t0)/h;
 +
%definimos la variable t
 +
t = linspace(t0,tN,N+1);
 +
%funcion y
 +
y = zeros(1,N+1);
 +
ya = zeros(1,N+1);
 +
yb = zeros(1,N+1);
 +
y(1) = y0;
 +
 +
for i=1:N
 +
    y(i+1) = y(i)+h*(cos(t(i)) - y(i));
 +
end
 +
 +
ya = 3.* t;
 +
yb = 1.* t;
 +
 +
%solucion grafica
 +
hold on
 +
plot(t,y)
 +
plot(t,ya)
 +
plot(t,yb)
 +
hold off
 +
}}
  
  

Revisión del 12:55 12 feb 2015

Trabajo realizado por estudiantes
Título Reacciones complejas. Grupo 1A
Asignatura Ecuaciones Diferenciales
Curso Curso 2014-15
Autores María Ramírez

Ignacio Posada

Antonio López-Mateos

Pablo Bueno

Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura


%PVI por Euler con una h = 0.1
clear all
% y'= cos(t) - y
%condiciones dadas
y0 = 0;
t0 = 0;
tN = 2;
h = 0.1;
%calculamos N
N = (tN - t0)/h;
%definimos la variable t
t = linspace(t0,tN,N+1);
%funcion y
y = zeros(1,N+1);
ya = zeros(1,N+1);
yb = zeros(1,N+1);
y(1) = y0;

for i=1:N
    y(i+1) = y(i)+h*(cos(t(i)) - y(i));
end

ya = 3.* t;
yb = 1.* t; 

%solucion grafica
hold on
plot(t,y)
plot(t,ya)
plot(t,yb)
hold off