clear;
t=0;
t1=80; 
model=1;
if model==2
    x=[1,0];
end
if model==1
    x=[0,1.925,0.5];
end
h=0.005
i=1;
while t<t1
    loeul1(i)=x(1);
    loeul2(i)=x(2);
    if model==1
        loeul3(i)=x(3);
        x=x+h*lorenz(t,x);
    end
    if model==2
            x=x+h*harm_osc(t,x);
    end
    t=t+h;
    i=i+1;
end
if model==1
    plot3(loeul1,loeul2,loeul3)
else
    plot(loeul1,loeul2)
end