be brief

classic Classic list List threaded Threaded
7 messages Options
Reply | Threaded
Open this post in threaded view
|

be brief

‪Nastaran Rad‬ ‪
hello.
plz be brief and short my program that writed with MATLAB language.
i am waiting for you
my program:

 
     f1=inline('-3.9*m*t/(z^2)-1.2*m*a/(z^2)-14.5*m*b/(z^2)-1.34*m*s/(z^2)','m','t','a','b','s');
 
     z=1+3.3a+13.6b+33.3f+215c
    

     f2=inline('-3.9*m*t/(z^2)','m','t');
   
     z=1+3.3a+13.6b+33.3f+215c
    

     f3=inline('3.9*m*t/(z^2)-0.9*b/(z^2)-1.2*m*a/(z^2)','m','t','b','a')
     z=1+3.3a+13.6b+33.3f+215c
   

     f4=inline('2.3*b/(z^2)-1.34*m*s/(z^2)','b','m','s')
     z=1+3.3a+13.6b+33.3f+215c  


     f5=inline('0.9*a/(z^2)-14.5*m*b/(z^2)-2.3*b/(z^2)','a','m','b')
    
     z=1+3.3a+13.6b+33.3f+215c 

     f6=inline('1.2*a*m/(z^2)+14.5*m*b/(z^2)+1.34*m*f/(z^2)','a','m','b','s')

     z=1+3.3a+13.6b+33.3f+215c

 
       
     m(1)=3.8;
     t(1)=6.2;
     a(1)=0;
     s(1)=0;
     b(1)=0;
     c(1)=0;
     h=0.01;
     y=0:h:20;

max=0;%max = Maximum for Selectivity of Para-xylene%
max2=0;%max2 = Maximum for Selectivity of Meta-xylene%
Obtimized_Length=1;
Obtimized_Length2=1;
for i=1:(length(y)-1)
    k11 = h*f1(m(i),t(i),a(i),b(i),s(i),c(i));
    k12 = h*f2(m(i),t(i),a(i),b(i),s(i),c(i));
    k13 = h*f3(m(i),t(i),a(i),b(i),s(i),c(i));
    k14 = h*f4(m(i),t(i),a(i),b(i),s(i),c(i));
    k15 = h*f5(m(i),t(i),a(i),b(i),s(i),c(i));
    k16 = h*f6(m(i),t(i),a(i),b(i),s(i),c(i));
   
    k21 = h*f1(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k22 = h*f2(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k23 = h*f3(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k24 = h*f4(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k25 = h*f5(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k26 = h*f6(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
   
    k31 = h*f1(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k32 = h*f2(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k33 = h*f3(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k34 = h*f4(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k35 = h*f5(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k36 = h*f6(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);

    k41 = h*f1(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k42 = h*f2(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k43 = h*f3(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k44 = h*f4(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k45 = h*f5(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k46 = h*f6(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);


    m(i+1) = m(i) + (1/6)*(k11+2*k21+2*k31+k41);
    t(i+1) = t(i) + (1/6)*(k12+2*k22+2*k32+k42);
    a(i+1) = a(i) + (1/6)*(k13+2*k23+2*k33+k43);
    b(i+1) = b(i) + (1/6)*(k14+2*k24+2*k34+k44);
    s(i+1) = s(i) + (1/6)*(k15+2*k25+2*k35+k45);
    c(i+1) = c(i) + (1/6)*(k16+2*k26+2*k36+k46);

   

    R1(i+1) = f1(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R2(i+1) = f2(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R3(i+1) = f3(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R4(i+1) = f4(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R5(i+1) = f5(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R6(i+1) = f6(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));

spx(i)=a(i+1)/c(i+1);
    if spx(i)>max
        max=spx(i);
        i=Obtimized_Length;
    end
    smx(i)=b(i+1)/c(i+1);
    if smx(i)>max2
        max2=smx(i);
        i=Obtimized_Length2;
       
    end
end
disp('max is = ');
disp(max);
disp('max2 is = ');
disp(max2);
disp('Obtimized Length is = ');
disp(Obtimized_Length);
disp('Obtimized Length2 is = ');
disp(Obtimized_Length2);
    

_______________________________________________
Help-octave mailing list
[hidden email]
https://lists.gnu.org/mailman/listinfo/help-octave
Reply | Threaded
Open this post in threaded view
|

Re: be brief

Francesco Potortì
>hello.plz be brief and short my program that writed with MATLAB
>language.i am waiting for youmy program:

Uh?
--
Francesco Potortì (ricercatore)        Voice:  +39.050.621.3058
ISTI - Area della ricerca CNR          Mobile: +39.348.8283.107
via G. Moruzzi 1, I-56124 Pisa         Skype:  wnlabisti
(entrance 20, 1st floor, room C71)     Web:    http://fly.isti.cnr.it


_______________________________________________
Help-octave mailing list
[hidden email]
https://lists.gnu.org/mailman/listinfo/help-octave
Reply | Threaded
Open this post in threaded view
|

Re: be brief

Andreas Weber-4
In reply to this post by ‪Nastaran Rad‬ ‪
Am 21.07.2015 um 23:37 schrieb ‪Nastaran Rad‬ ‪:
> plz be brief and short my program that writed with MATLAB language.
> i am waiting for you

How much do you want to pay for this?
-- Andy

_______________________________________________
Help-octave mailing list
[hidden email]
https://lists.gnu.org/mailman/listinfo/help-octave
Reply | Threaded
Open this post in threaded view
|

Re: be brief

Thomas D. Dean-2
In reply to this post by ‪Nastaran Rad‬ ‪
<snip>

Are you trying to use octave?  What version?  Have you read the first
part of the octave manual?

Tom Dean

_______________________________________________
Help-octave mailing list
[hidden email]
https://lists.gnu.org/mailman/listinfo/help-octave
Reply | Threaded
Open this post in threaded view
|

Re: be brief

nrjank
In reply to this post by ‪Nastaran Rad‬ ‪
On Tue, Jul 21, 2015 at 5:37 PM, <[hidden email]> wrote:
hello.
plz be brief and short my program that writed with MATLAB language.
i am waiting for you
my program:

 
     f1=inline('-3.9*m*t/(z^2)-1.2*m*a/(z^2)-14.5*m*b/(z^2)-1.34*m*s/(z^2)','m','t','a','b','s');
 
     z=1+3.3a+13.6b+33.3f+215c
    

     f2=inline('-3.9*m*t/(z^2)','m','t');
   
     z=1+3.3a+13.6b+33.3f+215c
    

     f3=inline('3.9*m*t/(z^2)-0.9*b/(z^2)-1.2*m*a/(z^2)','m','t','b','a')
     z=1+3.3a+13.6b+33.3f+215c
   

     f4=inline('2.3*b/(z^2)-1.34*m*s/(z^2)','b','m','s')
     z=1+3.3a+13.6b+33.3f+215c  


     f5=inline('0.9*a/(z^2)-14.5*m*b/(z^2)-2.3*b/(z^2)','a','m','b')
    
     z=1+3.3a+13.6b+33.3f+215c 

     f6=inline('1.2*a*m/(z^2)+14.5*m*b/(z^2)+1.34*m*f/(z^2)','a','m','b','s')

     z=1+3.3a+13.6b+33.3f+215c

 
       
     m(1)=3.8;
     t(1)=6.2;
     a(1)=0;
     s(1)=0;
     b(1)=0;
     c(1)=0;
     h=0.01;
     y=0:h:20;

max=0;%max = Maximum for Selectivity of Para-xylene%
max2=0;%max2 = Maximum for Selectivity of Meta-xylene%
Obtimized_Length=1;
Obtimized_Length2=1;
for i=1:(length(y)-1)
    k11 = h*f1(m(i),t(i),a(i),b(i),s(i),c(i));
    k12 = h*f2(m(i),t(i),a(i),b(i),s(i),c(i));
    k13 = h*f3(m(i),t(i),a(i),b(i),s(i),c(i));
    k14 = h*f4(m(i),t(i),a(i),b(i),s(i),c(i));
    k15 = h*f5(m(i),t(i),a(i),b(i),s(i),c(i));
    k16 = h*f6(m(i),t(i),a(i),b(i),s(i),c(i));
   
    k21 = h*f1(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k22 = h*f2(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k23 = h*f3(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k24 = h*f4(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k25 = h*f5(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k26 = h*f6(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
   
    k31 = h*f1(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k32 = h*f2(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k33 = h*f3(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k34 = h*f4(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k35 = h*f5(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k36 = h*f6(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);

    k41 = h*f1(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k42 = h*f2(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k43 = h*f3(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k44 = h*f4(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k45 = h*f5(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k46 = h*f6(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);


    m(i+1) = m(i) + (1/6)*(k11+2*k21+2*k31+k41);
    t(i+1) = t(i) + (1/6)*(k12+2*k22+2*k32+k42);
    a(i+1) = a(i) + (1/6)*(k13+2*k23+2*k33+k43);
    b(i+1) = b(i) + (1/6)*(k14+2*k24+2*k34+k44);
    s(i+1) = s(i) + (1/6)*(k15+2*k25+2*k35+k45);
    c(i+1) = c(i) + (1/6)*(k16+2*k26+2*k36+k46);

   

    R1(i+1) = f1(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R2(i+1) = f2(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R3(i+1) = f3(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R4(i+1) = f4(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R5(i+1) = f5(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R6(i+1) = f6(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));

spx(i)=a(i+1)/c(i+1);
    if spx(i)>max
        max=spx(i);
        i=Obtimized_Length;
    end
    smx(i)=b(i+1)/c(i+1);
    if smx(i)>max2
        max2=smx(i);
        i=Obtimized_Length2;
       
    end
end
disp('max is = ');
disp(max);
disp('max2 is = ');
disp(max2);
disp('Obtimized Length is = ');
disp(Obtimized_Length);
disp('Obtimized Length2 is = ');
disp(Obtimized_Length2);
    

_______________________________________________
Help-octave mailing list
[hidden email]
https://lists.gnu.org/mailman/listinfo/help-octave



That program would not run in matlab either.  (maybe with the Symbolic toolbox? but this is Octave, not Matlab)

you use variables before you assign values to them. this is not a symbolic algebra package. if you move the lines from m(1)=3.8 to y=0:h:20  up to the front of the program, at least all of your variable will be properly declared. Next, you cannot write :

3.3a

you have to write:

3.3*a

for all of the z equations.

next, why are you writing the same z = equation six times?  Last, the z equations use a variable f that you never declare.

Three may be more errors, but this is the most obvious initial problem that will prevent the program from even running.

Nick J.

_______________________________________________
Help-octave mailing list
[hidden email]
https://lists.gnu.org/mailman/listinfo/help-octave
Reply | Threaded
Open this post in threaded view
|

Re: be brief

nrjank
On Thu, Jul 23, 2015 at 6:38 AM, <[hidden email]> wrote:

From: Nicholas Jankowski <[hidden email]>
To: [hidden email]
Cc: "[hidden email]" <[hidden email]>
Sent: Wednesday, 22 July 2015, 19:15:56
Subject: Re: be brief

On Tue, Jul 21, 2015 at 5:37 PM, <[hidden email]> wrote:


hello.
plz be brief and short my program that writed with MATLAB language.
i am waiting for you
my program:

 
     f1=inline('-3.9*m*t/(z^2)-1.2*m*a/(z^2)-14.5*m*b/(z^2)-1.34*m*s/(z^2)','m','t','a','b','s');
 
     z=1+3.3a+13.6b+33.3f+215c
    

     f2=inline('-3.9*m*t/(z^2)','m','t');
   
     z=1+3.3a+13.6b+33.3f+215c
    

     f3=inline('3.9*m*t/(z^2)-0.9*b/(z^2)-1.2*m*a/(z^2)','m','t','b','a')
     z=1+3.3a+13.6b+33.3f+215c
   

     f4=inline('2.3*b/(z^2)-1.34*m*s/(z^2)','b','m','s')
     z=1+3.3a+13.6b+33.3f+215c  


     f5=inline('0.9*a/(z^2)-14.5*m*b/(z^2)-2.3*b/(z^2)','a','m','b')
    
     z=1+3.3a+13.6b+33.3f+215c 

     f6=inline('1.2*a*m/(z^2)+14.5*m*b/(z^2)+1.34*m*f/(z^2)','a','m','b','s')

     z=1+3.3a+13.6b+33.3f+215c

 
       
     m(1)=3.8;
     t(1)=6.2;
     a(1)=0;
     s(1)=0;
     b(1)=0;
     c(1)=0;
     h=0.01;
     y=0:h:20;

max=0;%max = Maximum for Selectivity of Para-xylene%
max2=0;%max2 = Maximum for Selectivity of Meta-xylene%
Obtimized_Length=1;
Obtimized_Length2=1;
for i=1:(length(y)-1)
    k11 = h*f1(m(i),t(i),a(i),b(i),s(i),c(i));
    k12 = h*f2(m(i),t(i),a(i),b(i),s(i),c(i));
    k13 = h*f3(m(i),t(i),a(i),b(i),s(i),c(i));
    k14 = h*f4(m(i),t(i),a(i),b(i),s(i),c(i));
    k15 = h*f5(m(i),t(i),a(i),b(i),s(i),c(i));
    k16 = h*f6(m(i),t(i),a(i),b(i),s(i),c(i));
   
    k21 = h*f1(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k22 = h*f2(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k23 = h*f3(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k24 = h*f4(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k25 = h*f5(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
    k26 = h*f6(m(i)+k11/2,t(i)+k12/2,a(i)+k13/2,b(i)+k14/2,s(i)+k15/2,c(i)+k16/2);
   
    k31 = h*f1(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k32 = h*f2(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k33 = h*f3(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k34 = h*f4(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k35 = h*f5(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);
    k36 = h*f6(m(i)+k21/2,t(i)+k22/2,a(i)+k23/2,b(i)+k24/2,s(i)+k25/2,c(i)+k26/2);

    k41 = h*f1(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k42 = h*f2(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k43 = h*f3(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k44 = h*f4(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k45 = h*f5(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);
    k46 = h*f6(m(i)+k31,t(i)+k32,a(i)+k33,b(i)+k34,s(i)+k35,c(i)+k36);


    m(i+1) = m(i) + (1/6)*(k11+2*k21+2*k31+k41);
    t(i+1) = t(i) + (1/6)*(k12+2*k22+2*k32+k42);
    a(i+1) = a(i) + (1/6)*(k13+2*k23+2*k33+k43);
    b(i+1) = b(i) + (1/6)*(k14+2*k24+2*k34+k44);
    s(i+1) = s(i) + (1/6)*(k15+2*k25+2*k35+k45);
    c(i+1) = c(i) + (1/6)*(k16+2*k26+2*k36+k46);

   

    R1(i+1) = f1(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R2(i+1) = f2(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R3(i+1) = f3(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R4(i+1) = f4(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R5(i+1) = f5(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));
    R6(i+1) = f6(m(i+1),t(i+1),a(i+1),b(i+1),s(i+1),c(i+1));

spx(i)=a(i+1)/c(i+1);
    if spx(i)>max
        max=spx(i);
        i=Obtimized_Length;
    end
    smx(i)=b(i+1)/c(i+1);
    if smx(i)>max2
        max2=smx(i);
        i=Obtimized_Length2;
       
    end
end
disp('max is = ');
disp(max);
disp('max2 is = ');
disp(max2);
disp('Obtimized Length is = ');
disp(Obtimized_Length);
disp('Obtimized Length2 is = ');
disp(Obtimized_Length2);
    

_______________________________________________
Help-octave mailing list
[hidden email]
https://lists.gnu.org/mailman/listinfo/help-octave



That program would not run in matlab either.  (maybe with the Symbolic toolbox? but this is Octave, not Matlab)

you use variables before you assign values to them. this is not a symbolic algebra package. if you move the lines from m(1)=3.8 to y=0:h:20  up to the front of the program, at least all of your variable will be properly declared. Next, you cannot write :

3.3a

you have to write:

3.3*a

for all of the z equations.

next, why are you writing the same z = equation six times?  Last, the z equations use a variable f that you never declare.

Three may be more errors, but this is the most obvious initial problem that will prevent the program from even running.

Nick J.



> hello.thank u for your answer. it is matLAB Language.f1 to f6   which is one function that other place is
> writed.but i write it whole of program(all   of program).i want that underestand you.i want brief all 
> program, delete Additional commands for example k41 to 46.(i want that write all of in 1
> command).thanks

Please include help-octave.gnu.org in any replies. And the list expects your replies to appear at the end of the other messages.

It seems that English is not your first language. There may be someone on the list able to speak to you directly in your native language, please let us know what that is, and maybe try to repeat your question in that language.

Here's the best I can do:

What you have written does resemble Matlab and  Octave m-code. However, it does not run in either.

After correcting the 3.3a to 3.3*a as I said above,

the first error MATLAB gives you is:

undefined function or variable 'a'. 

this is because you're trying to use a before you've defined it. same with other variables.

if you move the following lines to the top of the file:
     m(1)=3.8;
     t(1)=6.2;
     a(1)=0;
     s(1)=0;
     b(1)=0;
     c(1)=0;
     h=0.01;
     y=0:h:20;

then, you get the error:

undefined function or variable 'f'

if you want it to be f1, f2, f3, ...  then you need to use that name in the equations for z.

Assuming you meant to use f1, f2, f3, f4, f5 and f6 in the equations for z, such as:

 z=1+3.3*a+13.6*b+33.3*f1+215*c

then, you get the error: 
Undefined operator '*' for input arguments of type 'inline'.

It does not like when you multiply by f1.


Because I do not know what problem you are trying to solve, I have no idea how to correct your model. As I said, it currently doesn't even run in Matlab. 

If it is making use of the Symbolic Toolbox, then maybe it would run, but I believe you still have to change the setup and the way you declare variables

Nick J.


_______________________________________________
Help-octave mailing list
[hidden email]
https://lists.gnu.org/mailman/listinfo/help-octave
Reply | Threaded
Open this post in threaded view
|

Re: be brief

fbarbuto
On 23-07-2015 20:03, Nicholas Jankowski wrote:
> hello.thank u for your answer. it is matLAB Language.f1 to f6 which is one function that
> other place is writed.but i write it whole of program(all of program).i want that underestand you.i want
> brief all program, delete Additional commands for example k41 to 46.(i want that write all of in 1
> command).thanks
>

As I understand it, he wants to have his program shortened by eliminating
incorrect/unnecessary lines -- a tough thing to do when the objective of
the code is not known.

F.




_______________________________________________
Help-octave mailing list
[hidden email]
https://lists.gnu.org/mailman/listinfo/help-octave