Thread Subject:

passing a polynomial to a function

Hello,

is there a way of passing a polynomial to a function. i have a polynomial as follows:

p= [A.^2 (2.*A.*C-1) ((2.*A.*B)+(C.^2)+D) 2.*B.*C B.^2];

which i pass into another function as:

[m]=newton_method(p,m0);

but when it is passed to the function it acquires some random value which is not the value of the polynomial that has been passed...

is there a way of doing that?

bye

Snow White wrote:
> Hello,
>
> is there a way of passing a polynomial to a function. i have a
> polynomial as follows:
>
> p= [A.^2 (2.*A.*C-1) ((2.*A.*B)+(C.^2)+D) 2.*B.*C B.^2];
...

What are the dimensions of A, B, ...?

If they're scalars, then p is a vector.

If they're (column) vectors, then p is an array of length(A) rows by 5.
  If they're (row) vectors, the p is a row vector of length(A)*5.
I guess you can extrapolate to matrices... :)

Short story is, evaluate

  p= [A.^2 (2.*A.*C-1) ((2.*A.*B)+(C.^2)+D) 2.*B.*C B.^2]

at the command line and all should be revealed (note no trailing ";")

--

A B C D are all column vectors.

dpb <none@non.net> wrote in message <i8vptr$kvn$1@news.eternal-september.org>...
> Snow White wrote:
> > Hello,
> >
> > is there a way of passing a polynomial to a function. i have a
> > polynomial as follows:
> >
> > p= [A.^2 (2.*A.*C-1) ((2.*A.*B)+(C.^2)+D) 2.*B.*C B.^2];
> ...
>
> What are the dimensions of A, B, ...?
>
> If they're scalars, then p is a vector.
>
> If they're (column) vectors, then p is an array of length(A) rows by 5.
> If they're (row) vectors, the p is a row vector of length(A)*5.
> I guess you can extrapolate to matrices... :)
>
> Short story is, evaluate
>
> p= [A.^2 (2.*A.*C-1) ((2.*A.*B)+(C.^2)+D) 2.*B.*C B.^2]
>
> at the command line and all should be revealed (note no trailing ";")
>
> --

Snow White wrote:
...[top posting repaired...don't do that; hard makes conversation follow]...

> dpb <none@non.net> wrote in message
> <i8vptr$kvn$1@news.eternal-september.org>...
>> Snow White wrote:
>> > Hello,
>> > > is there a way of passing a polynomial to a function. i have a >
>> polynomial as follows:
>> > > p= [A.^2 (2.*A.*C-1) ((2.*A.*B)+(C.^2)+D) 2.*B.*C B.^2];
>> ...
>>
>> What are the dimensions of A, B, ...?
>>
>> If they're scalars, then p is a vector.
>>
>> If they're (column) vectors, then p is an array of length(A) rows by
>> 5. If they're (row) vectors, the p is a row vector of length(A)*5.
>> I guess you can extrapolate to matrices... :)
...

> A B C D are all column vectors.

So, p is a 2D array.

What's the problem?

Are you searching for

doc polyvalm % , perhaps?

--

dpb <none@non.net> wrote in message <i91kgm$kc5$1@news.eternal-september.org>...
> Snow White wrote:
> ...[top posting repaired...don't do that; hard makes conversation follow]...
>
> > dpb <none@non.net> wrote in message
> > <i8vptr$kvn$1@news.eternal-september.org>...
> >> Snow White wrote:
> >> > Hello,
> >> > > is there a way of passing a polynomial to a function. i have a >
> >> polynomial as follows:
> >> > > p= [A.^2 (2.*A.*C-1) ((2.*A.*B)+(C.^2)+D) 2.*B.*C B.^2];
> >> ...
> >>
> >> What are the dimensions of A, B, ...?
> >>
> >> If they're scalars, then p is a vector.
> >>
> >> If they're (column) vectors, then p is an array of length(A) rows by
> >> 5. If they're (row) vectors, the p is a row vector of length(A)*5.
> >> I guess you can extrapolate to matrices... :)
> ...
>
> > A B C D are all column vectors.
>
> So, p is a 2D array.
>
> What's the problem?
>
> Are you searching for
>
> doc polyvalm % , perhaps?
>
> --
p is n x5

when i pass this to a function the values that i get are garbbage!

Snow White wrote:
> dpb <none@non.net> wrote in message
> <i91kgm$kc5$1@news.eternal-september.org>...
>> Snow White wrote:
>> ...[top posting repaired...don't do that; hard makes conversation
>> follow]...
>>
>> > dpb <none@non.net> wrote in message >
>> <i8vptr$kvn$1@news.eternal-september.org>...
>> >> Snow White wrote:
>> >> > Hello,
>> >> > > is there a way of passing a polynomial to a function. i have a
>> > >> polynomial as follows:
>> >> > > p= [A.^2 (2.*A.*C-1) ((2.*A.*B)+(C.^2)+D) 2.*B.*C B.^2];
>> >> ...
>> >>
>> >> What are the dimensions of A, B, ...?
>> >>
>> >> If they're scalars, then p is a vector.
>> >>
>> >> If they're (column) vectors, then p is an array of length(A) rows
>> by >> 5. If they're (row) vectors, the p is a row vector of
>> length(A)*5.
>> >> I guess you can extrapolate to matrices... :)
>> ...
>>
>> > A B C D are all column vectors.
>>
>> So, p is a 2D array.
>>
>> What's the problem?
>>
>> Are you searching for
>>
>> doc polyvalm % , perhaps?
>>
>> --
> p is n x5
>
> when i pass this to a function the values that i get are garbbage!

Can't see what you don't show...

--

"dpb" <none@non.net> wrote in message
news:i91rds$qu3$2@news.eternal-september.org...
> Snow White wrote:
>> dpb <none@non.net> wrote in message
>> <i91kgm$kc5$1@news.eternal-september.org>...
>>> Snow White wrote:
>>> ...[top posting repaired...don't do that; hard makes conversation
>>> follow]...
>>>
>>> > dpb <none@non.net> wrote in message >
>>> <i8vptr$kvn$1@news.eternal-september.org>...
>>> >> Snow White wrote:
>>> >> > Hello,
>>> >> > > is there a way of passing a polynomial to a function. i have a
>>> > >> polynomial as follows:
>>> >> > > p= [A.^2 (2.*A.*C-1) ((2.*A.*B)+(C.^2)+D) 2.*B.*C B.^2];
>>> >> ...
>>> >>
>>> >> What are the dimensions of A, B, ...?
>>> >>
>>> >> If they're scalars, then p is a vector.
>>> >>
>>> >> If they're (column) vectors, then p is an array of length(A) rows
>>> by >> 5. If they're (row) vectors, the p is a row vector of
>>> length(A)*5.
>>> >> I guess you can extrapolate to matrices... :)
>>> ...
>>>
>>> > A B C D are all column vectors.
>>>
>>> So, p is a 2D array.
>>>
>>> What's the problem?
>>>
>>> Are you searching for
>>>
>>> doc polyvalm % , perhaps?
>>>
>>> --
>> p is n x5
>>
>> when i pass this to a function the values that i get are garbbage!
>
From doc for polyvalm:

Y = polyvalm(p,X) evaluates a polynomial in a matrix sense. This is the same
as substituting matrix X in the polynomial p.

Polynomial p is a vector whose elements are the coefficients of a polynomial
in descending powers, and X must be a square matrix.

Your p is a matrix, not a vector.
--
Loren
http://blogs.mathworks.com/loren/
http://matlabwiki.mathworks.com/MATLAB_FAQ