help-octave
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Must all coefficients for fsolve be numeric???????


From: Doug Stewart
Subject: Re: Must all coefficients for fsolve be numeric???????
Date: Tue, 17 Mar 2020 19:50:37 -0400



On Tue, Mar 17, 2020 at 11:44 AM ttonon <address@hidden> wrote:
Thanks Kai, I should've known better.  Here's the relevant code, where the
... lines indicate missing lines.  Let me know if you need more detail, and
thanks again!  - Tom

...
...
>> Y_e = cosh(Alpha*X_bar) - cos(Alpha*X_bar) + Eta*(sinh(Alpha*X_bar) -
>> sin(Alpha*X_bar))
Y_e =  1.5875
...
...
>> function y = f (x)

Octave is case sensitive  and you have a small x above this line and a large X below this line!!!
Change to large X every where!

y = zeros (12, 1);
y(1) = Y_e*X(1)*sin(X(2)) +dStat*X(5)/2 + O;
y(2) = X(4) -X(3) + Og_o;
y(3) = X(7);
y(4) = X(1)*Tau^2*(Og_o + X(3))*X(4) +Cf*X(8);
y(5) = X(9) -X(3)*Tau*Y_e*X(1);
y(6) = X(10) -abs(X(3)*Tau*Y_e*X(1))*L_cm/Tau;
y(7) = X(11) -W_cm/X(10);
y(8) = X(12) -Thic_cm/X(10);
y(9) = X(6) -Ogto*X(3)*X(11);
y(10) = X(13) -Mpi*(-D1)*((X(3)*tau*Y_e*X(1)))*(1 - sin(X(2)));
...
...
y(51) = X(8) -X(39) -X(40) -X(41) -X(42) -X(43) -X(44) -X(45) -X(46) -X(47)
-X(48) -X(49) -X(50)
 -X(51);
endfunction
>> [x, fval, info] = fsolve (@f,
>> [0.1;0.01;1024;4;-1.0;0.001;0.0;0.5;71000;0.000007;0.000005])
error: 'Y_e' undefined near line 3 column 8
error: called from
    f at line 3 column 6
    fsolve at line 236 column 8
>> y(1) = 1.5875*X(1)*sin(X(2)) +dStat*X(5)/2 + O;
error: 'X' undefined near line 1 column 15
>> y(1) = Y_e*X(1)*sin(X(2)) +dStat*X(5)/2 + O;
error: 'X' undefined near line 1 column 12
>> [x, fval, info] = fsolve (@f,
>> [0.1;0.01;1024;4;-1.0;0.001;0.0;0.5;71000;0.000007;0.000005])
error: 'Y_e' undefined near line 3 column 8
error: called from
    f at line 3 column 6
    fsolve at line 236 column 8
>>
>>




--
Sent from: https://octave.1599824.n4.nabble.com/Octave-General-f1599825.html




--
DASCertificate for 206392


reply via email to

[Prev in Thread] Current Thread [Next in Thread]