zjervoj97

40 Reputation

4 Badges

4 years, 95 days

MaplePrimes Activity


These are questions asked by zjervoj97

Hi,

I'd like to know, if it is possible to define any sort of range for parameters in NonlinearFit. E. g. I know that one of parameters should be somewhere between 0.2 - 0.4. I know there is a possibility of initalvalues, but using it doesn't lead into this range.

Thanks.

Hi,

I need to find out how to print the number of row or the numbers on right based only on the combination of the two numbers in first two columns. The external file looks like this:

I have prepared possible combinations of numbers and I need to read the numbers on the right which belongs to this combination. There is no such thing as the sum of two numbers = number of row, so it can't be done this easy.

Can somebody help please? Thanks a lot

Hello, 

is there any possibility how to write just a part of axis label e.g. in Italic? For example if I want want a label to be "ln(h)" and i want just the "h" to be written in Italic font.

Thanks a lot.

Hello, 

when trying to solve equation using fsolve, it returns just what i wrote and didn't solve. Can anyone help? I tried all i know. File is attached below.

Thanks a lot.

II.01.mw

 

Hi, 

when running a for do cycle the first two solutions are exact the same numbers, which is not correct. Can someone plesa help me fix this problem?

Thanks a lot.


 

restart

tmin:=10: #°C

tmax:=70: #°C

Tmin:=tmin+273.15:

Tmax:=tmax+273.15:

T:='T':

i:=0:

 

for T from Tmin by 2 to Tmax+1.25 do
i:=i+1:
a12:=2305.28444347652 - 9.14490843016421*T + 0.00680052257590234*T^2:
a21:=-6665.24838284836 + 46.0897018087247*T - 0.0694991633494123*T^2:
alfa:=0.3:
x2:=1-x1:
lng1:=x2^2*(tau21*(G21/(x1+x2*G21))^2+tau12*(G12/((x2+x1*G12)^2))):
lng2:=x1^2*(tau12*(G12/(x2+x1*G12))^2+tau21*(G21/((x1+x2*G21)^2))):
G12:=exp(-alfa*tau12):
G21:=exp(-alfa*tau21):
tau12:=a12/T:
tau21:=a21/T:
lnga1:=subs(x1=xa1,lng1):
lngb1:=subs(x1=xb1,lng1):
lnga2:=subs(x1=xa1,lng2):
lngb2:=subs(x1=xb1,lng2):
r1:=lnga1+ln(xa1)=lngb1+ln(xb1):
r2:=lnga2+ln(1-xa1)=lngb2+ln(1-xb1):
r:=fsolve({r1,r2},{xa1=0..0.22,xb1=0.22..1}):
assign(r):
xA1[i]:=xa1: xB1[i]:=xb1:
print(i,T,xA1[i],xB1[i]);
unassign('r1','r2','r','xa1','xb1'):
od:

1, 283.15, 0.4086195757e-1, .7650860222

 

2, 285.15, 0.4086195757e-1, .7650860222

 

3, 287.15, 0.4109304612e-1, .7564109743

 

4, 289.15, 0.4094687104e-1, .7455613500

 

5, 291.15, 0.4092946881e-1, .7345259571

 

6, 293.15, 0.4103899093e-1, .7233142804

 

7, 295.15, 0.4127502108e-1, .7119337099

 

8, 297.15, 0.4163857511e-1, .7003892047

 

9, 299.15, 0.4213210541e-1, .6886831535

 

10, 301.15, 0.4275956094e-1, .6768151230

 

11, 303.15, 0.4352648605e-1, .6647816285

 

12, 305.15, 0.4444014482e-1, .6525760020

 

13, 307.15, 0.4550971743e-1, .6401880864

 

14, 309.15, 0.4674655982e-1, .6276039487

 

15, 311.15, 0.4816452203e-1, .6148056135

 

16, 313.15, 0.4978038925e-1, .6017705621

 

17, 315.15, 0.5161445760e-1, .5884711603

 

18, 317.15, 0.5369126350e-1, .5748739966

 

19, 319.15, 0.5604059678e-1, .5609387793

 

20, 321.15, 0.5869882131e-1, .5466171048

 

21, 323.15, 0.6171072215e-1, .5318505354

 

22, 325.15, 0.6513208559e-1, .5165679944

 

23, 327.15, 0.6903339246e-1, .5006820935

 

24, 329.15, 0.7350533414e-1, .4840835134

 

25, 331.15, 0.7866732744e-1, .4666324128

 

26, 333.15, 0.8468127263e-1, .4481446382

 

27, 335.15, 0.9177527543e-1, .4283678555

 

28, 337.15, .1002878317, .4069372077

 

29, 339.15, .1107594115, .3832836608

 

30, 341.15, .1241550103, .3564111567

 

31, 343.15, .1425634130, .3241972432

(1)

 


 

Download lle.mw

lle.mw

1 2 Page 1 of 2