From: a a on
o 12 musim skoncit, nmeda viac casu
From: e on
"a a" <toth.aniko(a)azet.sk> wrote in message <i0uqog$49t$1(a)fred.mathworks.com>...
> o 12 musim skoncit, nmeda viac casu

to stihnem, prvy mam uz skoro hotovy
From: e on
"a a" <toth.aniko(a)azet.sk> wrote in message <i0uqog$49t$1(a)fred.mathworks.com>...
> o 12 musim skoncit, nmeda viac casu
zachvilu sem dam kod, aj ked nebude asi nic pocitat ale malo by to byt dobre
From: e on
"a a" <toth.aniko(a)azet.sk> wrote in message <i0uqog$49t$1(a)fred.mathworks.com>...
> o 12 musim skoncit, nmeda viac casu

zachvilu sem dam kod, aj ked nie je hotovy cely
From: e on
"a a" <toth.aniko(a)azet.sk> wrote in message <i0uqog$49t$1(a)fred.mathworks.com>...
> o 12 musim skoncit, nmeda viac casu
hlavny prog

clear all
close all
clc
global pol Y1 X1 P nI nL Kpa XL YL


pol=[5.437 0.9842 0.0028];

Y1=0.1;

X1=0.0496;

P=101325;

Kpa=4.6306e-004;

nL=250;
V=8;
nI=180;
YL=polyval(pol,XL);

[Y,Vn]=ode15s('ved2',[0 V],Y1);
Y2=Y(end)
plot(Y,Vn)

vedlajsi program ved 2

function dVdY = ved2(Y,Vn)
global pol Y1 X1 P nI nL Kpa XL YL
XL=X1-nI*(Y1-Y)/nL;
YL=polyval(pol,XL);
dVdY=-nI*(1+Y)*(1+YL)/Kpa/P/(Y-YL);