Dobrý den, mám skript z matlabu, ale za boha nemohu udělat výsledný graf s logaritmickou osou y. Mohl by mi někdo poradit prosím?
n=60;
h=0.1;
f='(-log(11/2)/100)*(y-1)';
%f='[cos(x); x^2]'
%f='[y(2); y(1)+y(2)+x]'
t0=0;
y0=12
t=t0;
yp=y0;
hold on
grid on
fp=inline(f,'t','y')
for i=1:n
k1=feval(fp,t,yp);
ypom=yp+h/2*k1;
k2=feval(fp,t+h/2,ypom);
yn=yp+h*k2;
semilogy(yn,t,'.');
yp=yn;
t=t+h;
%disp(sprintf('%6.4f ',x,yn));
end;
xlabel('xIn [-]')
ylabel('t [s]')
Offline