## 220 Reputation

16 years, 286 days
Glasgow, United Kingdom

## Thx...

I really appreciate your help and patients. Do you use Tecplot or have any information about it?

## Re: some explanation of a few points...

`1. "By plotting it is 'obvious', that for r=0 this is always negative (so rmust decrease to have a chance to fullfill F=0."`
`The above point you mentioned should be understood as `
`when R3=r=0, this means that the gravity factor is ignored, so it is now`
` irrelavant that whether liquid is above vapour or vapour is above liquid.`
` Bcoz If gravity is zero then the problem, instead of a vertical through flow `
`becomes a  horizental flow problem. which is always stable thats why we have `
`sigma<0, for all l=lambda. `
` `
`2. "`
`For lambda = 10^3 and 10^4 compute solutions for F=0 with an initialguess of r = -2*lambda which then defines a line equationLin := lambda -> -1.591517842027199990210998791356366512040*lambda                 +2.229559088349405643936812337536203057778`
`"`
`This means that for large l=lambda the system is marignally stable, i.e., sigma=0`
` for large lambda. to me thats fine.`
` `

## re: 3d data plot...

this what you want?

 > restart:
 > with(plots):
 > Filelocation:=("c:/s.txt"); #ur data file location
 (1)
 > datalist := readdata(`Filelocation`,float,3): #where 3 is the number of columns

 > plots[pointplot3d](datalist, style=line, axes=boxed );
 >
 >
 >
 >

## re:...

bcoz u didnt assign the value to m

ur code is

m = 0.1083e-7;

which should be

m := 0.1083e-7;

its just a typo

## re:thanks...

Well thanks

I have sloved it thats why I put the correct form of the solution in the maple sheet.

(1) T1=1+fraction is bcoz I am trying to plot in the same way as given in  the paper. orginaly it is T1=fraction.

(2) in T2 there is no "E" it is typo in the paper.

I am looking to find a way to plot the results in a more sensiable way, which can be understand easly.

## Re...

This is not a problem of particular case, I am experiencing  again and again the same problem whenever I try to solve a differential equation.

## re...

restart:

fsolve(int(x^2+4,x=a..7)=40,a);

 Page 1 of 1
﻿